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The  infrared  intensity  data  for  the  methyl  halides  have 
been  interpreted  in  terms  of  the  atomic  polar  tensors  first 
introduced  by  Biarge,  Herranz,  and  Morcillo  (1961).  By  using 
the  principal  of  isotopic  invariance  and  the  approximate 
quantum  mechanical  CNDO/2  technique,  a unique  sign  choice  for 
the  3p/9Q's  for  CH^F  and  CH^Cl  was  determined.  By  analogy 
the  signs  for  the  9p/3Q's  for  CH  Br  and  CH  I were  inferred. 

The  errors  in  the  polar  tensor  elements  were  estimated  by 
comparing  the  CH^X  and  CD^X  tensors.  In  addition,  a normal 
coordinate  analysis  of  the  CH^X  and  CD^X  molecules  was 
carried  out  to  insure  consistency  in  coordinate  definitions 
using  the  most  recent  GF  force  fields  of  Duncan  et  al.  (1972) . 


(AP) 


The  magnitudes  of  the  dipole  moment  change  vectors, 
a 

, associated  with  the  atomic  polar  tensors  were  shown 
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to  have  smooth  trends  through  the  CH3X  series,  suggesting 
that  these  vectors  may  be  predictable  and  transferable  in- 
tensity parameters.  Similar  trends  were  observed  for  the 

H 

angle  between  the  (AP)  vectors  and  the  C-H  bond. 

The  effective  charges,  , first  defined  by  King,  Mast 
and  Blanchette (1972)  for  the  gth  atom  in  the  CH^X  series  were 
computed.  From  the  effective  charges  of  CH^F,  the  effective 
charges  for  CF^ , CF^H,  CF^H^ , and  CH^  were  predicted  and 
compared  to  the  experimentally  observed  intensity  sums  for 
these  molecules. 

The  CNDO/2  results  for  CH  F and  CH  Cl  were  critically 

J 

analyzed  and  the  absolute  intensities,  A.  , predicted  by 
CNDO , are  compared  to  the  experimental  data.  In  general, 
the  x-esults  for  CII^F  are  much  better  than  for  CH^Cl.  For 
example,  we  found  that  the  CNDO/2  calculated  bond  moment  for 
C-Cl  is  positive,  which  seems  contrary  to  chemical  intuition. 
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CHAPTER  1 


INTRODUCTION 

Infrared  spectroscopy  has  long  been  recognized  as  an 
aid  in  the  elucidation  of  the  structures  of  complex  molecu- 
lar systems  (1) . The  facts  that  different  spatial  arrange- 
ments of  the  atoms  give  rise  to  different  spectral  features 
and  that  in  many  cases  the  spectral  features  are  sensitive 
to  the  intramolecular  (as  well  as  intermolecular)  environ- 
ments, together  with  the  relative  ease  with  which  the 
spectra  are  obtained,  have  made  the  infrared  spectrometer 
a common  tool  in  most  chemical  laboratories.  In  the  past, 
the  emphasis  has  been  on  the  interpretation  of  the  spectra 
in  terms  of  frequencies  of  functional  groups  or  of  chemical 
bonds,  and  extensive  correlation  charts  have  been  derived 
that  are  quite  useful  for  qualitative  analysis.  Moreover, 
the  theoretical  basis  for  such  correlations  is  firmly 
established  (see  Reference  1,  for  example) . Much  less  atten- 
tion has  been  given  to  the  infrared  band  intensities  to 
obtain  the  additional  useful  information  about  the  molecular 
environment  which  they  may  contain  (2) . The  infrared  band 
intensity  is  a function  of  the  change  in  the  dipole  moment 
when  the  molecule  is  distorted  from  equilibrium;  it  is  a 
direct  reflection  of  the  changes  in  the  electron  distribu- 
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tion  in  the  molecule  as  a whole.  As  a consequence,  the  band 
intensity  is  expected  to  reflect  any  distortion  of  the  elec- 
tron cloud  produced  by  substituents  and/or  by  the  immediate 
environment  of  the  molecule.  Wexler  (2)  has  developed  an 
empirical  correlation  chart  for  band  intensities,  similar 
to  the  frequency  correlation  charts , that  demonstrates  quite 
well  the  potential  usefulness  of  the  band  intensities  as  an 
additional  probe  for  the  study  of  molecular  structure  and 
environments . 

In  general,  the  use  of  infrared  band  intensities  has 
been  limited  to  a semi-quantitative  tool  with  such  ambiguous 
terms  as  "strong  (s) , " "medium  (m) , ° or  "weak  (w) " used  to 
describe  the  band  intensities.  There  are  understandable 
reasons  for  the  lack  of  quantitative  infrared  band  intensity 
data  for  all  but  the  simplest  molecules.  Reliable  measure- 
ments of  the  absolute  (or  even  of  relative)  band  intensities 
is  not  a routine  procedure,  although  the  technique  has  been 
fairly  well  standardized  (3) . This  problem  may  well  be  elim- 
inated with  the  advent  of  the  computerized  spectrometer,  and 
reliable  band  intensity  measurements  may  soon  become  as  rou- 
tinely available  as  the  spectrum  itself.  There  is,  however, 
a more  fundamental  problem  associated  with  the  routine  use 
of  infrared  intensity  data;  that  is,  the  lack  of  a uniform, 
coherent  theoretical  basis  for  interpreting  them  (4). 

Hence,  while  it  seems  that  the  measurement  of  accurate 
band  intensities  may  be  fast  approaching  the  "routine,"  it 
now  seems  appropriate  that  the  theoretical  tools  should  be 
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developed  in  order  to  offer  a platform  for  interpretation. 

The  Morcillo  Theory 

One  of  the  primary  purposes  of  this  thesis  is  to  illus- 
trate a procedure  for  the  reduction  of  the  experimental  in- 
tensity data  to  a form  which  we  believe  is  susceptible  to 
theoretical  scrutiny  and  interpretation.  The  procedure 
was  first  introduced  by  Biarge,  Herranz , and  Morcillo  (5) 
and  independently  rediscovered  by  Person,  Peyton,  and  Jiang 
(6) , but  has  attracted  little  attention  since  its  conception 
in  1961.  With  the  exception  of  two  revealing  papers  by 
Morcillo  and  his  co-workers  on  the  application  of  the  technicrue 
(7  and  8) , little  work  has  been  done  employing  the  procedure. 

As  a consequence,  the  usefulness  of  the  procedure  as  a theo- 
retical basis  and  an  interpretative  tool  has  been  relatively 
unexplored. 

The  technique  is  discussed  in  Chapter  2,  where  we  refer 
to  it  as  the  "polar  tensor"  or  "Morcillo  technique."  It 
involves  the  reduction  of  the  experimentally  determined  com- 
plete infrared  band  intensity  data  to  a set  of  "atomic  polar 
tensors";  one  for  each  atom  in  the  molecule,  defined  in  a 
space-fixed  Cartesian  coordinate  system  for  the  molecule. 

The  polar  tensor  for  the  ath  atom  in  the  molecule  is  defined 
as  (5)  : 


x a 


a 


8p  /3y 
y a 


3p  /3z 
Y a 


(1-1) 


4 ' 

Here  a stands  for  the  ath  atom  in  the  molecule  and  the  ath 
tensor  is  composed  of  elements  measuring  the  change  in  the 
dipole  moment  components  with  respect  to  an  infinitesimal 
change  along  the  Cartesian  axes  (x,  y,  or  z)  on  the  ath  atom. 

The  atomic  polar  tensors  quantitatively  reflect  the 
contributions  of  each  atom  to  the  total  intensity  of  the 
infrared  bands  in  the  molecule.  In  the  past,  the  band 
intensities  have  been  interpreted  in  terms  of.  bond  dipole 
moments  and/or  of  dipole  moment  changes  with  respect  to  a 
set  of  defined  internal  coordinates  (3).  The  procedures, 
though  useful  in  their  own  way,  have  not  proved  to  be 
generally  satisfactory.  The  bond  moment  hypothesis,  for 
example,  which  assumes  that  the  dipole  moment  change  occurs 
parallel  to  the  bond  for  a bond-stretching  motion  and  per- 
pendicular to  the  bonds  for  an  angle-bending  motion,  and 
which  might  result  in  an  attractively  simple  theory,  has 
never  been  accepted  (9) . 

Since  the  polar  tensor  technique  resolves  the  band  in- 
tensities into  atomic  contributions,  we  are  able  to  examine 
atom  by  atom  the  contributions  to  the  intensities.  In  essence 
the  polar  tensors  offer  a more  detailed  breakdown  of  the  in- 
frared band  intensities  that  we  believe  is  one  of  its  major 
advantages. 

One  of  the  more  significant  advances  in  terms  of  the 
theoretical  interpretations  of  band  intensities  in  the  last 
decade  has  been  the  application  of  ab  initio  and  approximate 


quantum  mechanical  techniques  to  the  calculation  of  band  in- 
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tensities  (4).  As  we  shall  show,  the  polar  tensor  technique 
presents  the  experimental  data  in  a simple,  straightforward 
form  that  can  readily  be  compared  with  quantum  mechanical 
calculations  of  theoretical  polar  tensors. 

The  most  significant  and  potentially  useful  feature  of 
the  polar  tensor  technique  is  that  the  associated  vectors 
and  individual  elements  of  the  tensors  in  a series  of  struc- 
turally similar  molecules  exhibit  simple  and  useful  trends. 
The  fact  that  certain  polar  tensor  elements  are  quite  in- 
sensitive to  changes  in  substituents  and  others  vary  in  a 
predictable  manner  suggests  that  it  may  be  possible  to  pre- 
dict quantitatively  the  infrared  spectra  of  more  complex 
systems  from  a knowledge  of  the  polar  tensors  of  simpler 
systems . 

King,  Mast,  and  Blanchette  (KMB)  have  defined  a quantity 

which  they  called  the  effective  charge,  £ , for  the  ath  atom 

a — 

in  a molecule  (10) , and  which  bears  a simple  relationship  to 

the  intensity  sum  rule  (11).  They  found  that  the  effective 

charge  of  the  hydrogen  atoms  in  a series  of  hydrocarbons  was 

constant.  It  has  been  shown  that  the  squared  effective 

charge  derived  by  KMB  is  equal  to  the  trace  of  the  product 

of  the  ath  atomic  polar  tensor  with  its  transpose  (12).  This 

tensor  invarient,  E,  , appears  to  be  another  transferable  or 

a 

predictable  intensity  parameter,  as  we  shall  show  later  in 
Chapter  4. 
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In frarea  Gas  Phase  Band  Intensities 

As  discussed  by  Overend  (3) , the  experimentally  deter- 
mined intensity  (lb)  of  the  ith  fundamental  absorption  band 

* i ^ r 

is  related  to  the  transition  moment  < 'H  „ p ¥ , > by: 

v v 

r.  = (1/n  SL)J  {In  (I  /I)  } d(lnv) 

1 J band  0 

= 87r^N/3hc  < Y „ | p | Y . (1-2) 

v v 1 

Here  n is  the  concentration,  % is  the  pathlength,  I is  the 
intensity  of  the  incident  radiation  of  wavenumber,  v,  and  I 
is  the  intensity  of  the  radiation  after  passage  through  the 

•k 

sample;  ¥ and  ¥ are  the  wave  functions  for  the  ground 
v"  v ' 

and  excited  vibrational  states  respectively  (and  are  usually 

A 

taken  to  be  harmonic  oscillator  functions) , and  p is  the 
dipole  moment  operator.  The  concentration  term  N and  h,  and 
c have  their  usual  meanings.  Assuming  both  mechanical  and 
electric  harmonicity.  Equation  1-2  reduces  to  the  simple 
form  for  the  gas  phase  infrared  band  intensities: 

r±  = (Ntt/3c2V  ) Op/8Qi)2  , (1-3) 

where  v.  is  the  harmonic  wavenumber  (in  cm  "S  of  the  ith 

i — 

fundamental  band.  The  quantity  3p/3Q.  is  the  dipole  moment 
derivative  relative  to  the  ith  normal  coordinate,  Q..  As 

— l 

Overend  has  pointed  out,  even  though  Equation  1-3  is  approxi- 
mate, the  errors  involved  in  the  approximations  are  generally 
regarded  as  negligible  relative  to  the  error  involved  in  the 
actual  integration  of  the  spectrum  (3) , and  Equation  1-3  is 
quite  sufficient  for  the  interpretation  of  molecular  systems 
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of  general  interest. 

The  dimensions  of  V are  (concentration  x length) 

Another  vibrational  intensity  measure,  A,  is  often  used, 
and  it  is  related  to  F by: 

A.  - T.V.  = (Ntt/3c2)  (3p/3Q.  ) 2 , (1-4) 

X 1 1 i 

and  has  dimensions  of  (wavenumber/ concentration  x length) . 

2 

The  units  of  Y are  generally  cm  /mol,  or  in  SI  units  (4) 

-3 

where  the  concentration  is  in  mol  m , the  pathlength  in 

-1  2 

m,  and  the  wavenumber  in  m ; F has  the  units  m /mol.  The 

-1 

SI  units  of  A are  m mol  where  the  pathlength  is  in  m and 

-3 

the  concentration  in  mol  m (Equation  1-2)  . In  this  work 

2 

we  shall  use  the  intensity  measurement  F in  units  of  cm  /mol. 

The  dipole  moment  derivatives  have  dimensions  of 

charge  x length/  length,  and  we  shall  use  the  units  of 
o-l  9-1 

DA  , or  electrons  (e)  (1.0  D A =0.208  e ).  In  terns  of 

SI  units  the  dipole  moment  is  expressed  in  Coulombs  meters 

(Cm),  with  1 esu  • cm  = 1 D = 3.335  x 10  ^ Cm.  Thus  DA 

0-1  -20 

expressed  m SI  units  becomes  1 DA  = 3.335  x 10  C = 

-2  -10  o-l  -2 

3.335  x 10  aC,  and  1 e =4.8  x 10  DA  = 16.02  x 10  aC. 

In  Equation  1-4  we  see  that  either  measure  of  the  ex- 
perimental intensity  is  directly  proportional  to  the  square 
of  the  dipole  moment  derivative.  Hence  Equation  1-4  gives 
the  relationship  between  these  experimentally  observed 
quantities  and  a set  of  molecular  parameters  subject  to 
theoretical  interpretation  and  computation.  Only  the  magni- 
tude  of  3p/3Cb  can  be  determined  from  such  experimental  in- 
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tensity  measurements;  if  the  band  intensity  is  to  be  in- 
terpreted in  terms  of  changes  in  the  electron  density  (or 
dipole  moment  derivatives  with  respect  to  internal  coordi- 
nates) , the  sign  of  9p/3Q.  must  be  determined  as  well.  As 
we  have  pointed  out,  the  data  are  usually  reduced  to  a set 
of  dipole  moment  derivatives  with  respect  to  a set  of  de- 
fined internal  displacement  coordinates  (3)  or,  as  in  the 
case  of  the  polar  tensors,  to'  the  set  of  Cartesian  displace- 
ment coordinates.  We  shall  not  pursue  this  subject  further 
since  the  procedures  are  described  in  detail  in  Chapter  2. 

It  will  suffice  here  to  point  out  that  even  when  the  data 
are  reduced  in  polar  tensor  forav  the  sign  ambiguity  is 
still  present  and  a complete  description  of  the  normal  coordi- 
nates is  still  required.  For  the  reduction  of  the  data,  or 
for  a description  of  if  3p/3Q_^  is  to  form  the  basis  of 
interpretation,  the  normal  coordinate  transformation  matrix, 

L,  is  required.  The  point  to  be  made  here  is  that  for  a 
meaningful  interpretation  of  infrared  band  intensities,  we 
must  not  only  have  accurate  intensity  measurements,  but  also 
some  rules  to  determine  the  signs  of  the  3p/3Qbs  and  a 
rather  complete  knowledge  of  the  force  field  with  a normal 
coordinate  analysis  of  the  molecular  system  to  give  L (13) . 
Sign  of  Dipole  Moment  Derivatives 

In  recent  years  the  application  of  ab  initio  and 
approximate  quantum  mechanical  techniques  has  proven  quite 
helpful  in  determining  the  signs  of  the  dipole  moment  deriva- 
tives. The  CNDO/2  approximate  quantum  technique  (14)  has  been 
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highly  successful,  and  a fairly  high  degree  of  confidence  has 
been  expressed  in  its  ability  to  predict  signs  for  molecules 
containing  first  row  elements.  This  subject  was  recently 
reviewed  in  some  detail  (4) . 

There  are  also  two  experimental  procedures  which  aid 
in  the  determination  of  the  relative  signs  of  the  dipole 
derivatives.  The  first  involves  a detailed  analysis  of  the 
Coriolis  perturbation  between  two  fundamental  vibrational 
modes  of  appropriate  symmetry;  from  these  studies  the  rela- 
tive signs  of  the  two  3p/3Q^'s  for  the  coupled  modes  can  be 
determined  (15) . The  second  method  involves  the  comparison 
of  the  integrated  band  intensities  of  a molecule  with  its 
isotopically  substituted  species  and  has  been  referred  to  as 
the  "isotopic  invariance  principle"  (11) . (We  shall  discuss 
these  two  techniques  in  more  detail  in  Chapter  3.)  Thus,  if 
the  pertinent  Coriolis  data  and  if  the  intensities  of  iso- 
topically substituted  molecules  are  both  available,  one  may 
eliminate  a large  number  of  possible  sign  combinations  of 
the  3p/3Q^'s.  These  two  procedures,  coupled  with  the  quantum 
mechanical  calculations,  make  the  choice  of  the  signs  of  the 
3p/3Q^'s  more  objective,  as  we  shall  demonstrate  in  Chapter  3, 
and  requires  less  reliance  on  chemical  intuition,  but  usually 
they  are  insufficient  to  determine  only  one  set  of  signs. 

The  Methyl  Halides 

The  methyl  halides  (CH3F,  CH3Cl,  CH-jBr,  and  CH  I)  are 
perhaps  the  best  series  of  structurally  similar  molecules 
for  use  as  a model  system  to  demonstrate  the  polar  tensor 
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technique  and  its  potential  usefulness. 

During  the  last  decade,  the  methyl  halides  have  been 
the  subject  of  several  extensive  studies,  especially  calcu- 
lations to  establish  their  force  fields.  Notable  among  these 
are  the  classic  papers  by  Aldous  and  Mills  that  appeared  in 
1962  (16)  and  1963  (17) , where  an  attempt  was  made  to  deter- 
mine the  general  harmonic  force  fields  for  the  molecules  in 
the  series.  In  their  force  constant  refinement  calculations, 
they  utilized  frequency  data  for  the  CH^X  and  CD^X  molecules, 
as  well  as  Coriolis  constants  and  centrifugal  distortion  con- 
stants for  both  CH^X  and  CD^X.  In  all,  twenty-one  pieces  of 
independent  experimental  data  were  used  to  calculate  twelve 
harmonic  force  constants  for  each  CH^X  molecule.  In  spite 
of  the  large  number  of  experimental  data  used,  they  were  not 
completely  successful  in  fixing  the  force  fields.  {The 
general  problem  involved  here  is  discussed  thoroughly  by 
Aldous  and  Mills  (16) , by  Mills  (13) , and  also  briefly  again 
in  Chapter  2.}  For  CH^F  they  found  that  two  off-diagonal 
force  constants  in  the  A^  class  were  essentially  undetermined 
and  that  the  E class  force  constants  converged  to  two  alter- 
native sets.  For  CH^Cl  the  force  field  did  not  converge,  but 
wandered  slowly  through  a family  of  solutions,  each  of  which 
gave  a good  fit  to  the  experimental  data.  For  CH^Br  and  for 
CH^I  one  set  of  force  constants  was  reported  for  each  mole- 
cule, but  these  sets  were  not  significant  since  they  were  just 
one  arbitrary  choice  from  a family  of  possible  solutions. 

In  1966,  Russell,  Needham,  and  Overend  (RNO)  (18)  carried 
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out  a force  constant  refinement  calculation  on  the  methyl 
halides,  using,  in  addition  to  the  data  used  by  Aldous  and 
Mills,  the  frequencies  of  the  CHD2X  and  CH^DX  molecules. 
Although  these  data  are  not  independent  (19),  the  calcula- 
tions converged  to  a unique  set  of  force  constants  for  each 
molecule.  However,  several  of  the  off-diagonal  force  con- 
stants had  such  large  dispersions  that  they  were  essentially 
undetermined. 

Duncan,  Allen,  and  McKean  in  1970  (20)  and  Duncan , McKean 

13 

and  Speirs  in  1972  (21)  used  isotope  frequencies  from  C 
substituted  molecules,  together  with  all  the  data  previously 
used  by  Overend,  finally  to  fix  the  general  harmonic  force 
constants  within  a fairly  narrow  uncertainty  range  for  all 
four  methyl  halides.  The  force  constants  found  by  Duncan, 
et  crl.  (21)  (hereafter  cited  as  Duncan)  were  quite  similar 
to  those  of  Overend,  but  with  significantly  smaller  disper- 
sions for  many  of  the  force  constants.  In  general,  the  dis- 
persions in  the  off-diagonal  force  constants  in  Duncan's 
force  field  were  an  order  of  magnitude  smaller  than  those 
for  the  force  field  of  RNO.  In  particular  the  off-diagonal 
force  constants  (f-,  0,  f ) that  were  undetermined  in  Over- 
end's  work  have  been  determined  fairly  well  by  Duncan's  work. 
(In  RNO's  work,  the  dispersions  for  these  force  constants 
ranged  from  100  to  200%  of  their  values,  whereas  in  Duncan's 
work  the  uncertainty  is  reduced  to  less  than  20  to  50%  of 
the  force  constant  for  the  series.) 


The  infrared  band  intensities  of  the  fundamentals  of 
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CH  X and  CD  X (where  X is  Cl,  Br,  or  I) , have  been  measured 
3 3 

fairly  accurately  by  Dickson,  Mills,  and  Crawford  (22).  The 
band  intensities  of  CH^F  an<^  CD^F  have  been  measured  by 
Barnett  and  Eggers  (as  reported  by  RNO) . With  the  exception 
of  CD^I,  the  intensity  measurements  were  complicated  by  over- 
lapping bands,  the  most  serious  of  which  occurred  in  CH^F 
and  CD^F . With  the  exception  of  the  fluoride,  however,  the 
errors  in  the  band  intensities,  including  the  estimated  errors 
involved  in  band  separation,  ranged  from  about  5 to  20%  for 

the  series.  For  CH-,F  and  CD  F each  fundamental  band  is  over- 

J 3 

lapped  with  one  other  fundamental,  and  the  intensities  of 

two  bands,  v_  of  CH_F  and  v_  of  CD  F,  have  estimated  errors 
2 3 5 3 

of  about  100%  and  200%  respectively.  The  other  bands  in  the 
fluoride  molecules  have  estimated  errors  of  about  10  to  20% 
(18) . Although  the  errors  for  the  methyl  halide  intensities 
may  appear  excessive,  they  are  of  the  magnitude  one  must 
expect  for  intensity  measurements  of  this  kind.  The  problem 
of  overlapping  band  separation  has  been  very  aptly  described: 
"...separation  is  a major  source  of  error,  demanding,  as  it 
usually  does,  a keen  eye,  a sharp  pencil,  a certain  artistic 
talent,  and  simple  faith  in  the  symmetry  of  vibrational- 
rotational  bands  and  the  cooperation  of  nature"  (18). 

Russell,  Needham,  and  Overend's  work  (18)  on  the 
methyl  halide  series  gives  the  best  previous  estimate  of  the 
3p/3Qi's  for  CH^X  and  CD^X  molecules.  Although  many  sign 
combinations  of  the  3p/3Q^'s  were  eliminated,  they  were  not 
able  to  select  a unique  set  of  3p/3Cb  's  for  any  molecule  in 
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the  series.  The  primary  purpose  of  their  work  was  to 
illustrate  the  procedure  for  calculating  the  errors  in 
dipole  moment  derivatives  and  normal  coordinates , and  this  pro- 
vides another  reason  for  our  choice  of  this  series  as  a 
model  system:  we  can  use  the  errors  reported  by  them  as  a 
starting  point,  and  propagate  them  through  to  the  polar 
tensors  to  give  a measure  of  their  reliability. 

In  the  remaining  chapters  of  this  thesis,  we  shall 
attempt  to  show  the  potential  of  the  polar  tensor  technique 
using  the  methyl  halide  series  as  a model  system.  Our  pri- 
mary purpose  will  be  to  illustrate  the  technique  and  show 
how  it  can  be  used  for  interpretive  purposes.  In  addition, 
we  shall  point  out  trends  in  the  tensors  that  we  believe  are 
indicative  of  the  predictive  nature  of  the  polar  tensors. 

We  shall  use  the  approximate  quantum  mechanical  technique, 
CNDO/2,  to  aid  in  our  choice  of  signs  for  the  8p/3Q^'s  and 
critically  analyze  the  results  of  the  calculations  for  CH^F 
and  CH  Cl.  In  doing  this  we  also  hope  to  illustrate  the 
potential  of  the  polar  tensor  technique  for  quantum  mechanical 
calculations  and  theoretical  interpretations. 


CHAPTER  2 


POLAR  TENSORS 

As  we  pointed  out  in  the  last  chapter,  the  polar  tensor 
concept,  though  fairly  old,  is  not  well  known  among  spectros- 
copists.  In  light  of  this  fact  we  shall  briefly  develop  in 
this  chapter  the  concept  of  an  atomic  polar  tensor  as  first 
introduced  by  Biarge,  Herranz , and  Morcillo  (5),  but  using 
the  revised  notation  of  Person  and  Newton  (12) . In  order 
to  develop  the  mathematical  expressions  for  the  polar  tensors 
and  to  familiarize  the  reader  with  the  notation  to  be  employ- 
ed, we  must  digress  momentarily  and  examine  the  mathematical 
transformations  used  in  a normal  coordinate  analysis  of 
molecular  systems. 

Normal  Coordinate  Transformations 

The  normal  coordinate  analysis  procedures  have  been  well 
documented  (13,  23,  24)  and  we  shall  outline  only  the  steps 
pertinent  to  the  development  of  the  polar  tensors.  For  a 
molecular  system  containing  N atoms  there  are  3N  - 6 (or 
3N  - 5 if  the  molecule  is  linear)  normal  modes  of  vibration. 
In  order  to  determine  the  nature  of  these  normal  modes,  one 
generally  begins  with  a space-fixed  Cartesian  coordinate 
system  on  which  the  instantaneous  location  of  individual 
atoms  are  given  by  the  Cartesian  displacement  coordinates. 
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These  3N  Cartesian  displacement  coordinates  can  be  represented 
by  a column  vector  X of  dimensions  3N  x 1 as  shown  below: 


Here  , y^,  and  z ^ refer  to  the  Cartesian  displacement 
coordinates  of  atom  i and  the  i subscript  runs  from  1 to  N. 
The  normal  coordinates  can  be  expressed  in  terms  of  this 
coordinate  system,  but  it  has  been  found  more  convenient  to 
give  the  normal  coordinates  in  terms  of  molecule-fixed  in- 
ternal coordinates  or  symmetry  coordinates. 

The  transformation  from  Cartesian  displacement  coordi- 
nates to  internal  displacement  coordinates  (r ^ , where  j 
identifies  the  jth  internal  coordinate)  is  achieved  by  the 
transformation  matrix  B. 

R = B X.  (2-1) 

Here  R is  a (3N  - 6)  xl  column  vector  of  the  independent 
internal  displacement  coordinates  (bond  stretches  or  bends) , 


16 


expressed  as  linear  combinations  of  the  Cartesian  displace- 
ment coordinates.  B is  a (3N  - 6)  x 3N  transformation  matrix 

with  elements  of  the  form  b. . = 3r./9x.  so  that  an  r element 

13  1 3 1 

of  R has  the  form 


3N  3N 

r,  = l (3r./3x.)  x.  = £ b.x.  . (2-2) 

j=l  1 3 3 j=l  13  3 


The  transformation  matrix  B (or  its  element,  b. .)  is  obtained 

!3 

by  inspection  of  the  geometry  of  the  molecule  and  a convenient 
prescription  for  the  evaluation  of  the  elements  has  been  given 
by  Wilson,  Decius,  and  Cross  (23).  It  should  be  pointed  out 
that  the  transformation  just  described  transforms  a set  of 
3N  space-fixed  Cartesian  displacement  coordinates  into  a set 
of  3N  - 6 molecule- fixed  internal  displacement  coordinates. 

The  space-fixed  coordinate  system  includes , in  addition  to 
the  coordinates  of  the  vibrational  degrees  of  freedom,  coordi- 
nates for  the  three  rotational  and  three  translational  degrees 
of  freedom,  for  in  the  space-fixed  axis  system,  the  displace- 
ments cause  the  molecules  to  translate  and  rotate  relative 
to  a fixed  axis  system.  In  the  molecule-fixed  axis  system, 
the  coordinate  axis  rotates  and  translates  with  the  molecule 
as  the  atoms  are  displaced,  and  only  the  3N  - 6 vibrational 
degrees  of  freedom  are  evident  in  this  system.  Although 
this  may  appear  to  be  trivial  now,  it  will  be  important  in 
our  development  of  the  polar  tensors  and  we  shall  return  to 
this  fact  later  in  this  chapter. 

In  order  to  take  advantage  of  any  symmetry  of  the 
molecule  we  can  transform  the  internal  coordinates  to  symmetry 
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coordinates  by  the  following  transformation: 


S = U R , ' (2-3) 

where  S is  a (3N  - 6)  x 1 column  vector  of  elements,  , 
composed  of  a linear  combination  of  the  internal  displacement 
coordinates,  r . U is  a (3N  - 6)  x (3N  - 6)  orthogonal 
transformation  matrix  obtained  by  considering  the  symmetry 
point  group  to  which  the  molecule  belongs  (25) . 

The  normal  coordinates  (Q  ) are  related  to  the  symmetry 
displacement  coordinates  by  : 

S = L Q , or  Q = L-1S  . (2-4) 


Here  Q is  a (3N  - 6)  x 1 column  vector  composed  of  elements, 

CK  , which  are  linear  combinations  of  the  symmetry  displace- 
ment coordinates.  L is  the  (3N  - 6)  x (3N  - 6)  normal  coordi- 
nate transformation  matrix  and  L 1 is  its  inverse.  The 

general  form  of  an  element  of  L is  L. . = 3S./8Q.  in  terms  of 

- 13  i 3 

symmetry  displacement  coordinates  or  1„  = 9r^/3Q^  and 

l^j  = 3x^/8Qj  in  terms  of  internal  displacement  coordinates 
and  Cartesian  displacement  coordinates,  respectively.  These 
three  L matrices  are  not  identical  and  the  matrix  is  generally 
denoted  (as  here)  by  different  symbols  for  each  set  of  dis- 
placement coordinates.  We  shall  not,  however,  make  an  explicit 
distinction  in  L for  the  different  coordinates  since  we  shall 
be  speaking  in  general,  or  the  particular  set  of  coordinates 
will  be  obvious  from  the  equations. 


The  evaluation  of  the  B matrix  and  the  U matrix  is 
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straightf orward , but  the  evaluation  of  L is  not,  in  general, 
a trivial  task,  since  it  depends  on  the  values'  of  the  force 
constants,  as  discussed  by  Mills  (13).  In  general,  these 
values  are  not  known  (see  Reference  13) . A great  deal  of 
effort  has  been  expended  in  trying  to  evaluate  good  normal 
coordinates  (that  is,  good  L matrices)  for  polyatomic  systems 
(see  Reference  16  or  17  for  examples) . Although  we  are  not 
explicitly  interested  in  the  details  of  the  calculations, 
the  results  of  such  calculations,  as  we  shall  show,  are  of 
great  importance  in  evaluating  and  interpreting  the  polar 
tensors.  Once  a set  of  internal  displacement  coordinates 
has  been  chosen,  the  L matrix  can  be  evaluated  from  the 
solution  of  the  eigenvalue  equation 


G F L = L A 


(2-5) 


Here  G is  the  symmetric  inverse  kinetic  energy  matrix  and 
is  related  to  the  B matrix  by 


, -1 

G = B M B 


(2-6) 


• . -1 

where  B is  the  transpose  of  B and  M is  a 3N  x 3N  diagonal 

matrix  whose  elements  are  the  reciprocal  masses  of  the  atoms 
of  the  molecule.  F is  the  (3N  - 6)  x (3N  - 6)  symmetric 
force  constant  matrix  comprised  of  elements,  f^j,  which 
constitute  the  "force  field"  (2V  = potential  energy  = the 

I 

quadratic  form,  R F R)  of  the  molecule.  The  A is  the  (3N  - 6) 
x (3N  - 6)  diagonal  eigenvalue  matrix  whose  elements,  , are 
the  frequency  parameters  which  are  related  to  the  (3N  - 6) 
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frequencies  of  the  normal  modes  of  vibration.  The  matrix  L, 
the  normal  coordinate  transformation  matrix,  is  the  eigen- 
vector matrix  of  the  GF  product  matrix.  Although  there  are 
several  well-known  mathematical  techniques  (26)  available 
for  solving  the  eigenvalue  problem  expressed  by  Equation  2-5, 
they  all  involve  the  direct  or  indirect  solution  of  a set  of 
vibrational  secular  equations  given  by  the  determinant 

I G F - X±E  | =0  , (2-7) 

where  E is  a (3N  - 6)  x (3N  - 6)  identity  matrix.  Once  the 
eigenvalues  have  been  determined,  then  the  eigenvectors  can 
be  computed. 

The  formulation  of  the  solution  to  the  eigenvalue  pro- 
blem is  straightforward,  but  in  practice  we  may  generally 
encounter  a major  difficulty.  As  explained  by  Mills  (13,  27), 
the  problem  is  that  for  most  molecules  we  have  more  unknown 
force  constants  {approximately  £n^  (ng  + l)/2,  where  n is  the 
number  of  coordinates  of  symmetry  type  S}  than  we  have  known 
A^  's.  This  situation  renders  the  secular  equation  unsolvable 
unless  we  (1)  constrain  the  force  field,  (2)  assume  values  for 
the  unknown  force  constants,  or  (3)  obtain  additional  indepen- 
dent experimental  data  such  as  frequencies  (A^'s)  from  iso- 
topically  substituted  molecules,  Coriolis  coupling  constants, 
centrifugal  distortion  constants,  or  other  data  of  this  type. 
The  point  to  be  made  is  that  if  we  constrain  the  force  field 
or  make  a bad  assumption  for  any  unknown  force  constants,  then 
these  constraints  and  assumptions  will  be  reflected  in  L.  In 
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short,  an  ill-defined  force  field  yields  an  ill-defined  set 
of  normal  coordinates.  Thus,  we  cannot  expect'  any  subsequent 
calculations  involving  L to  be  any  better  than  the  set  of 
force  constants  used  to  obtain  L.  Up  to  this  point,  we  have 
assumed  that  the  geometry  and  the  vibrational  frequencies  of 
the  molecule  are  known.  If  there  are  errors  in  either  of 
these  experimental  quantities,  then  as  RNO  have  demonstrated 
(18) , the  normal  coordinates  will  reflect  them  also.  Gen- 
erally this  kind  of  lack  of  knowledge  is  not  as  serious  a 
problem  as  the  one  to  obtain  the  complete  set  of  force  con- 
stants . 

We  have  outlined  above  the  relation  between  the  normal 
coordinate  and  a set  of  Cartesian  displacement  coordinates. 
The  inverse  transformations  from  the  noimal  coordinates  to 
the  Cartesian  displacement  coordinates  are  outlined  below. 

The  inverse  L matrix  is  related  to  L,  F,  and  A ^ by 
the  relationship 


A-V 


(2-8) 


i -1  “J- 

where  L is  the  transpose  of  L and  Q = L S_.  The  A is 
easily  evaluated  since  A is  a diagonal  matrix  and  F is  known 
(after  the  normal  coordinate  calculations) . The  internal 
coordinates  are  obtained  from  the  symmetry  coordinates  by 


R = U S , (2-9) 

• • -1  ' 
since  U is  an  orthogonal  matrix  (U  = U ) . Since  B is  not 

a square  matrix,  we  cannot  evaluate  B 1 by  direct  inversion. 
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Instead,  we  calculate  the  inverse  transformation  matrix  A to 
carry  the  internal  displacement  coordinates  back  to  a set  of 
Cartesian  displacement  coordinates  (X  = A R)  by  the  following: 

A = M-1  b'  G-1  . (2-10) 

-1  . • 

Here  M is  easily  obtained  since  M is  diagonal;  B is  just 

the  transpose  of  B,  and  the  inverse  of  G can  be  computed  by: 
G-1  = (L_1) ' (L-1)  . (2-11) 

Thus  we  can  write 

X = A R = A U'L  Q (2-12) 

for  the  transformation  from  internal  displacement  coordinates 
to  the  Cartesian  displacement  coordinates. 

There  is  a subtlety  in  this  last  coordinate  transforma- 
tion, however,  that  should  be  pointed  out.  A carries  R into 
a set  of  molecule- fixed  Cartesian  displacement  coordinates, 
whereas  by  using  B in  Equation  2-1  we  transformed  a set  of 
space-fixed  Cartesian  displacement  coordinates  into  a set  of 
molecule-fixed  internal  displacement  coordinates.  Thus,  the 
X vector  of  Equation  2-12  is  not  the  same  as  the  X vector  we 
started  with  in  Equation  2-1.  This  brings  us  back  to  the 
point  mentioned  earlier  in  this  section  about  transforming 
from  space-fixed  to  molecule-fixed  coordinate  systems,  and 
we  shall  now  explore  this  in  a little  more  detail. 

In  order  to  transform  a set  of  space-fixed  Cartesian 
coordinates  into  a set  of  space-fixed  internal  coordinates. 


22 


the  augmented  B matrix  is  used: 


(2-13) 


where  R is  the  (3N  - 6)  x 1 column  vector  of  internal  coordi- 
nates previously  defined  and  £ is  a 6 x 1 column  vector  of 
the  translations  and  rotations  of  the  whole  molecule.  B is 
the  transformation  matrix  defined  previously  and  Ms  a 
6 x 3N  rotation  matrix  relating  £ to  X.  The  inclusion  of 
the  £ (and  3_)  is  an  explicit  expression  of  the  Eckart  condi- 
tions (28)  for  the  space-fixed  molecular  coordinate  frame. 

A typical  8 matrix  for  the  ath  atom  is  given  below. 


Here  m is  the  mass  of  the  a atom;  M is  the  total  mass  of  the 
a 


molecule;  I , , I , 
x y 


and  I are  the  moments  of  inertia  about 

. . . t , (e)  (e)  (e) 

the  equilibrium  principal  axes;  x^  , y^  and  are  the 


coordinates  of  the  a atom  at  equilibrium,  measured  from  the 
origin  at  the  center  of  mass,  along  the  principal  axes  of 
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inertia  (12) . There  are  N such  6x3  matrices  juxtaposed  to 
form  3_.  The  p_  vector  is  just  X.  To  transform  a set  of 
molecule-fixed  internal  displacement  coordinates  back  to  a 
set  of  space-fixed  Cartesian  displacement  coordinates,  the  p 
vector  and  the  augmented  A matrix  must  be  employed  as  shown 
in  Equation  2-14: 


where  X is  now  in  terms  of  space-fixed  Cartesian  displacement 
coordinates  equivalent  to  those  in  Equation  2-1.  The  product 
ap  forms  the  rotational  and  translational  corrections  necessary 
for  transforming  from  the  molecule-fixed  reference  frame  to 
the  space-fixed  reference  frame.  A summary  of  the  normal 
coordinate  transformation  notation  used  in  this  section  is 
given  in  Table  2-1. 

Atomic  Polar  Tensors 

As  we  indicated  in  Chapter  1 (Equation  1-2) , the  inte- 
grated infrared  band  intensity  of  a normal  vibrational  mode 
is  related  to  the  square  of  the  dipole  moment  derivative 
with  respect  to  the  normal  coordinate  (3p/3Q.)  of  the  particu- 
lar mode.  In  the  following  discussion,  we  shall  derive  the 
atomic  polar  tensors  for  a molecular  system  in  terms  of 
space-fixed  Cartesian  displacement  coordinates  starting  with 
the  experimentally  determined  3p/3Q^'s.  For  each  normal 
coordinate,  Q.,  we  can  write  a 3 x 1 column  vector  in  terms 


(2-14) 


A,  R,  and  p have  been  defined  and  a = M 6.  Thus 


X = A R + ap 


(2-15) 


l 
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TABLE  2-1 


SUMMARY  OF  NOTATION  FOR  A NORMAL  COORDINATE  TRANSFORMATION 


Matrix 

Typical  Element 

Dimensions 

(Unweighted) 

X 

x.  , y.  , z. 

X X 1 

(length) 

R,  S 

a 

(s) , si  (s) 

(length) 

b 

r . (b)  . s . (b) 

1 1 

— 

B 

b.  . (s) 

il 

— 

b (b) 

kl 

(length)  ^ 

-1 

M 

1/nr 

-1 

(mass) 

L 

L.  . (s) 

il 

(mass)  2 

(b) 

_P  - 

(mass)  2 (length) 

Q 

Q. 

i 

% 

(mass)  (length) 

GC 

g. . (s,s) 

il 

-1 

(mass) 

g±j  (sfb) 


g . . (b,b) 
11 


-1  -1 

(mass)  (length) 


-1  v -2 
(mass)  (length) 


Continued  on  next  page 
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(Table  1-1  continued) 


Matrix  Typical  Element  Dimensions 

(Unweighted) 


f . . 
11 

(s,  s) 

(energy) 

(length) 

f. 

Ij 

(s,b) 

(energy) 

(length) 

f 

j j 

(b,b) 

(energy) 

A, 

l 

-2 

(time) 

» 


a.  s indicates  a bond  stretching  motion. 

b.  b indicates  an  angle  bending  motion. 

c.  We  shall  use  a weighted  G matrix  so  that  the  s-b  ^nd  s-s 
elements  are  weighted  by  d = 1.0  A and  d = (1.0A) 
respectively,  and  the  force  constants  have  units  of 
Mdynes/S  for  all  f elements. 

ID 
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of  the  x,  y,  and  z components  of  the  dipole  moment  derivative 
as  described  by  Morcillo  and  co-workers  (5) . 


a.  9p  /9Q. 
ixi 


^i 


ck  9p  /9Q^  (charge)  (mass) 


-% 


(2-16) 


°i  3V3Qi 


Here  is  the  sign  (+  or  -)  of  the  dipole  moment  derivative 
with  respect  to  the  ith  normal  coordinate  and  the  subscripts 
x,  y,  and  z refer  to  the  respective  components  of  the  dipole 
moment  change  in  the  molecule-fixed  Cartesian  coordinate 
frame . 

For  molecules  with  reasonable  symmetry,  the  dipole 
moment  change  will  occur  parallel  to  only  one  of  the  coordi- 
nate axes  so  that  only  one  component,  say  9p  /9Q  , will  be 

x i 

non-zero.  For  some  molecular  systems  this  may  not  be  true 

and  in  such  cases  one  would  have  to  resolve  the  9p/9Q^ 

vector  into  its  respective  components.  Unless  the  particular 

vibrational  mode  in  question  happens  to  be  a degenerate  mode, 

it  is  not  at  once  obvious  how  the  experimentally  determined 

9p/9Q.  should  be  partitioned  into  its  components  (3p  /3Q., 
i ' x 1 

3p„/3Q. , 9p  /9Q.)  and  in  such  cases  this  may  negate  the 
y i z i 

advantages  of  the  interpretation  of  infrared  intensities  by 
the  polar  tensor  technique.  For  the  molecules  we  shall 
discuss,  and  for  a large  number  of  other  molecules,  this  is 
not  a problem  and  in  the  following  development  of  the  polar 
tensors  we  shall  assume  that  this  is  always  true. 
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We  can  form  a 3 x (3N  - 6)  matrix,  P^,  by  the  juxta- 
position of  the  3N  - 6 column  vectors,  P^  , so  that 

0 -i 


P = fP  P P 

-O  -Qi  ' -02  * -Q3 

• • • 


• En  } 

~^3N  - 6 


(2-17) 


We  can  express  the  dipole  moment  derivatives  in  terms 
of  symmetry  coordinates  by 

-1 


P = P^  L 

-S  ~Q  “ 


(charge) 


(2-18) 


where  P is  a 3 x (3N  - 6)  matrix  with  elements  of  the  form 
-S 

9p^_/9S..  A typical  element  of  P^  is  given  by 

3N  - 6 

*S(T,k)  = 8V8Sk= 


3N  - 6 -i 

1 <*_.  (3p_/3Q..)Lik 


(2-19) 


i=l  1 T 

where  T represents  x,  y,  or  z,  and  k = 1,  2,  3 . . . . 3N  - 6. 

Transformation  to  dipole  derivatives  with  respect  to 
internal  displacement  coordinates  is  accomplished  by  the 
use  of  the  transformation  matrix  U. 


P = P U 
~ R ~S  “ 


(2-20) 


P is  a 3 x (3N  - 6)  matrix  containing  elements  of  the  form 
R 

9p  /9r  and  given  by 
t k 

3N  - 6 

*R(T,k)  = 3V»rk  = <3PT/3^(3S./3rk) 


3N  - 6 

jll  ^S(T#j)  Ujk 


(2-21) 


We  should  note  that  Equation  2-19  and  2-21  represent  the 
reduction  of  the  intensity  data  from  the  9p/9Q. 's  that  we 
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referred  to  in  Chapter  1.  The  reduction  of  the  data  to 

~y 

3p/3S_.'s  (or  3p/3R^'s)  and  their  use  for  subsequent  inter- 
pretation of  the  intensities  has  been  discussed  by  Overend 
(3)  and  in  the  past,  interpretation  of  infrared  intensities 
has  been  based  on  P and/or  on  P 


~S 


-R 


We  can  transform  the  P matrix  to  a 3 x 3N  matrix  con- 

— R 

taining  the  dipole  moment  derivatives  with  respect  to  mole- 
cule-fixed Cartesian  displacement  coordinates  using  the  B 
matrix.  This  transformation  is  given  by  Equation  2-22: 


V 

“X 


(2-22) 


In  general,  we  are  more  interested  in  the  dipole  moment 
change  with  respect  to  the  space-fixed  Cartesian  coordinates 
and  in  order  to  transform  to  that  coordinate  system  we  must 
use  the  augmented  B and  P matrices  as  given  by 

- -p 


= P B + P 8 . (2-23) 

R “p  “ 

Here_P  is  a 3 x 6 matrix  which  gives  the  dipole  moment  change 
P 

with  respect  to  translations  and  rotations  of  the  molecule. 

Since  a translation  of  the  molecule  does  not  change  the 

dipole  moment,  the  translational  terms  in  P are  zero.  Using 

-P 

the  relationship  described  by  Biarge,  Herranz,  and  Morcillo 

(5)  for  computing  3p/3p,  we  find  that  the  P matrix  has  the 

“P 


form 
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0 % 

P /I 
z y 


0 ^ 

-P  /I 
*z  x 


0 

-p  /I 

y 2 

o 

p„/i 


% 


o 

p /i 
y x 


o 

-p  /i 
x y 


0 / 


(2-24) 


The  quantities  1,1,  and  I are  the  three  principal  moments 

x y z 

of  inertia  assuming  that  x,  y,  and  z are  drawn  to  coincide 

with  these  principal  axes,  and  p°,  p°,  and  p^  are  the  x,  y, 

X y z 

and  z components  of  the  permanent  dipole  moment.  P 3 con- 

-p- 

stitutes  the  rotational  correction  to  the  dipole  moment 

derivatives  for  the  space-fixed  coordinate  frame  (11) . 

Notice  that  if  the  molecule  contains  no  permanent  dipole 

moment,  P becomes  the  null  matrix  and  no  rotational  correc- 
P 

tion  is  made. 

A typical  element  of  P^  , Sp^/Sx^, is  given  by 
3N 

-X(T,k)  = 


3N 

= I p„ , . . b ; k = 1,2,3... 3N 

j=l  R(t,3)  jk 


(2-25) 


The  P^  matrix  can  be  partitioned  into  a set  of  N.  3 x 3 
matrices  associated  with  each  atom  (designated  by  the  super- 
script) so  that 


r (1) 

P = {P 
-x  -X 


(2) 

P 

-X 


(3) 

P 

“X 


(N) 

P } 
“X 


(2-26) 


(a) 


The  submatrix,  P , is  the  desired  atomic  polar  tensor  for 

X. 
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atom  a and  has  the  form 


/ 3p  /3x 
x a 

3P  /3y 
x a 

3p  /3z 

X 

X P 
II 

3p  /3x 
y a 

3P  /3Ya 

y a 

3p  /3z 
Y 

3p  /3x 
\ z a 

3p  /3y 
z a 

3p  /3z 
z 

\ 

(2-27) 

a 


a / 


The  rotational  correction  (P  ' 8_)  is  not  explicit  in  Equation 

p a 

2-27  but  is  implied  by  our  notation.  In  writing  P or  P 

-x  -x 

it  is  understood  that  we  are  referring  to  a space-fixed 
coordinate  frame  and  the  rotational  corrections  are  implicitly 


included.  Notice  that  the  use  of  polar  tensors  carries  the 

reduction  of  the  data  in  the  form  of  3p/3Q^  one  step  further 

than  the  traditional  interpretation  in  terms  of  3p/3R. 's. 

The  inverse  transformations  from  P to  P can  be  carried 

-X  -Q 

out  by  employing  the  inverse  coordinate  transformations  given 

in  the  last  section.  Thus,  P can  be  computed  from  P by 

-R  -X 

right  multiplying  Equation  2-23  by  the  A matrix  so  that 
(PR  B + Pp  3_)A  = P^A.  Notice  that  since  3_  A = 0 (12)  the 
rotational  corrections  drop  out  so  that 


and 


P„  = P A , 
-R  -X  - 


P„  = P U . 
~S  -R  - 


Finally,  P is  given  by 
~Q 


P = P L 
~Q  S ~ 


P A U'L 


(2-28) 


(2-29) 


(2-30) 
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In  order  to  obtain  meaningful  atomic  polar  tensors 

we  need  accurate  intensity  measurements  and  a well  defined 

force  field;  moreover,  we  need  to  know  the  signs  (a's)  of 

the  3p/9Cb  's.  If  the  symmetry  properties  of  the  polar  tensors 

are  neglected  and  if  rotational  corrections  are  applied,  there 
„ (3N  - 6)  r (3N  - 5)  . 

are  2 {or  2 if  the  direction  of  the  permanent 

dipole  moment  is  known}  possible  sign  combinations  possible 
for  a set  of  (3N  - 6)  3p/9Q^'s.  In  reality,  the  symmetry 
properties  of  the  tensors  reduce  this  number  considerably,  as 
we  shall  demonstrate  in  the  next  chapter.  In  addition,  other 
sets  may  be  eliminated  by  use  of  Coriolis  coupling  data,  in- 
tensities from  isotopically  substituted  molecules,  or  chemical 
intuition  (using  the  latter  with  a great  deal  of  caution) . 

We  alluded  to  the  effective  charge,  £ , for  a given  atom, 

a 

a,  in  Chapter  1.  In  terms  of  the  P tensor,  the  effective 

— A 

charge  is 

K = (Tr  p“  (Pa)  , (2-31) 

a XX 

where  Tr  is  the  trace  of  the  product  matrix  pa(pa)  (12). 

X X 

The  associated  vectors  of  the  polar  tensors,  referred 
to  earlier,  are  the  dipole  moment  change  vectors  produced 
by  a unit  displacement  of  a given  atom  parallel  to  one  of  the 
Cartesian  coordinate  axes  (x,  y,  or  z) . That  is,  they  give 
the  magnitude  and  direction  of  the  dipole  moment  change  pro- 
duced by  the  motion  of  one  atom  in  a single  given  direction. 

We  shall  denote  these  vectors  by  (Ap)  where  the  superscript 
a indicates  the  atom  and  t indicates  the  coordinate  axis  (x. 
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->  C i 

y,  or  z) , along  which  atom  a is  moved.  The  (Ap)  ' vectors  are 

T 

defined  by 


+ a 
(Ap)x 


Op  /3t)  i , i 


Opy/3T)  j,j 


(8p  /3x)  k,k 
z 


t = x,  y,  or  z 


(6t)  . 


(2-32) 


The  quantity  in  brackets  is  just  the  xth  column  of  the  ath 

^ /s 

tensor,  and  i,  j,  k are  the  unit  vectors  parallel  to  the  x, 
y,  or  z coordinate  axes,  respectively.  The  quantity  x takes 
on  the  values  of  i,  j,  or  k depending  on  the  axis  under  con- 
sideration and  (6t)  is  a unit  displacement  parallel  to  the 
xth  axis.  These  vectors  are  shown  in  plots  given  in  Chapter 
4 for  the  methyl  halides,  and  their  direction  and  magnitude 
form  a useful  pictorial  representation  of  the  atomic  polar 
tensors . 

Another  useful  property  of  the  polar  tensors  derived 
by  Biarge,  Herranz , and  Morcillo  (5)  is 

„ a 

l P = 0.  (2-33) 

a x 

a 

Here  P is  the  polar  tensor  for  the  ath  atom  in  the  space- 
fixed  coordinate  frame  (rotationally  corrected  polar  tensor) . 
This  relationship  is  particularly  useful  since  it  shows  that 
there  are  only  N - 1 independent  tensors  (N  is  the  total 
number  of  atoms  in  the  molecule) . For  additional  properties 
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and  relationships  of  the  polar  tensors,  the  reader  is  referred 
to  the  papers  by  Biarge,  Herranz,  and  Morcillo'  (5)  , or  to 
Person  and  Newton  (12). 

Morcillo  (5)  has  discussed  the  fact  that  if  the  bond 
moment  hypothesis  were  true  one  would  expect  the  polar  ten- 
sors to  be  diagonal.  That  is 

- / . \ 

I 3p/3xa  0 0 

o 3p /3ya  0 

0 0 9p/3z  , 

/ 

and  the  dipole  moment  changes  would  only  be  induced  parallel 
to  the  direction  of  motion.  As  a consequence,  if  off-diagonal 
non-zero  elements  appear,  they  should  be  a reflection  of  the 
error  involved  in  the  bond  moment  hypothesis  and,  as  we  shall 
show  in  Chapter  4,  this  "error"  seems  to  be  quite  sensitive 
to  the  polarizability  of  adjacent  substituents  and  the  result- 
ing induced  dipole  moments . 

Signs  of  3p/9Q^'s 

At  this  point  we  must  interject  a word  of  warning  con- 
cerning the  arbitrary  nature  of  the  signs  of  the  3p/3Q.'s. 

To  illustrate  this  point,  consider  a triatomic  molecule  XY 

2 

°f  C symmetry  with  internal  displacement  coordinates,  r , 

, . ->0  1 
r^ , and  a,  and  a permanent  dipole  moment  p as  shown  in 

Figure  2-1 (a) . We  can  construct  two  symmetry  displacement 

coordinates,  S_^  and  S2,  corresponding  to  the  totally  symmetric 

vibrational  modes.  Two  ways  of  constructing  these  symmetry 

coordinates  (differing  only  in  relative  phase)  are  shown  in 


a 

P 

-X 
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Figures  2-1 (b)  and  2-1  (c) . Either  way  of  defining  the  symmetry 

coordinates  is  acceptable  since  the  only  criterion  for  these 

coordinates  is  that  they  transform  according  to  the  totally 

symmetric  representation  of  the  C point  group.  Now  assume 

2v 

that  a normal  coordinate,  Cb  , exists  that  is  composed  of  a 

certain  percentage  of  and  (the  percentage  contributions 

are  unimportant  here) . Associated  with  the  normal  coordinate 

is  a dipole  moment  change  Ap.  The  Q.  may  be  a linear  combi- 
b b 1 

nation  of  S_^  and  S^,  as  shown  in  Figure  2-1  (d)  . However,  an 
equivalent  statement  of  Q is  shown  in  Figure  2-1  (e)  as  a 
linear  combination  of  and  S Assuming  that  we  know  the 
direction  of  the  equilibrium  dipole  moment  as  indicated,  then 
the  sense  of  3p/3Q.  is  negative  in  Figure  2-1 (e)  if  Q.  is 

l l 

positive  in  that  figure;  and  the  sense  of  3p/3Q  is  positive 

i 

in  Figure  2-1 (d)  if  Q.  is  positive  in  that  figure.  The 


l 

problem  is,  how  do  we  know  whether  Q.  is  "positive"  or  "nega- 
tive" since  both  representations  of  Q.  are  equivalent  except 
for  a phase  difference?  There  is  no  a priori  reason  to 
choose  Figure  2-1 (b)  or  Figure  2-1 (c)  as  the  "positive"  dis- 
placement, but  it  is  obvious  that  the  resulting  sign  of  3p/3Q 
will  depend  on  just  that  choice.  In  reality,  we  are  not  at 
liberty  to  choose  2-1 (b)  or  2-1 (c)  as  positive;  the  choice 
is  made  for  us  by  the  coordinate  definitions  and  the  result- 
ing L matrix  from  the  normal  coordinate  analysis.  If  indeed 
there  is  any  meaning  to  the  sign  of  3p/3Q.,  it  is  only  within 
the  context  of  the  particular  coordinate  definitions  and 
that  particular  L matrix  from  the  normal  coordinate  calcula- 


Figure  2-1:  Schematic  diagram  of  two  ways  to  define 

totally  symmetric  symmetry  coordinates  of 
an  XY^  molecule. 

(a)  Equilibrium  configuration. 

(b)  and  (c)  Two  definitions  of  the 

totally  symmetric  vibrations 
differing  only  in  phase. 


(d)  and  (e)  The  resulting  two  definitions 
of  the  normal  coordinate. 
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(a) 


Si 


37 


tion.  As  a consequence  of  all  this,  if  one  speaks  of  the 
sign  of  a dipole  moment  derivative  with  respect  to  any  set 
of  coordinates  (X,  R,  £,  or  Q) , one  cannot  interpret  the 
results  unless  the  following  definitions  are  known : 

(1)  the  normal  coordinate  transformation  matrix  L or 


(2)  unambiguous  definitions  of  the  internal  coordinates 
R and 

and  (3)  an  unambiguous  definition  of  the  Cartesian  coordi- 
nate axis  and  the  direction  of  the  dipole  moment. 

It  is  unfortunate  that  in  many  papers  published  on 
normal  coordinate  analyses,  insufficient  information  is  given 
on  the  actual  coordinate  definitions  used.  This  means  that 
one  should  be  very  cautious  in  using  another's  published 
normal  coordinates  for  calculating  dipole  moment  derivatives. 
In  order  to  insure  consistency  in  coordinate  definitions,  it 
is  advisable  for  one  person  to  do  both  the  normal  coordinate 
calculation  and  the  dipole  moment  derivative  calculation. 

In  addition,  one  must  also  be  aware  of  the  difficulties 
that  may  arise  due  to  the  arbitrary  signs  of  the  eigenvectors 
of  GF_  (the  columns  of  L)  , particularly  when  comparisons  are 
made  between  independent  workers.  These  signs  are  quite  arbi- 
trary and  depend  on  the  diagonalization  and  normalization 
routines  used  in  the  normal  coordinate  analysis.  As  a conse- 
quence, even  if  two  independent  workers  use  identical  coordi- 
nate definitions  in  the  normal  coordinate  analysis,  they  may 
obtain  an  apparent  inconsistency  in  the  signs  of  the  GF  eigen- 
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vector  when  their  results  are  compared.  For  example,  we 
have  found  in. this  laboratory  that  we  obtain  different  signs 
for  the  columns  of  L for  the  and  symmetry  blocks  of 
F2CO  when  we  use  two  different  normal  coordinate  analysis 
computer  programs  {program  CHARLY  from  the  University  of 
Minnesota  (19),  and  Schachtschneider 1 s program  GVIB  (29)}, 
even  though  the  input  for  each  program  was  identical.  This 
fact  probably  contributed  to  some  of  the  difficulties  en- 
countered by  McKean,  Bruns,  Person,  and  Segal  (30)  in  their 
paper  on  F2C0* 

In  order  for  the  reader  not  to  be  misled,  it  is  important 
to  realize  that  the  sign  of  Q or  9p/9Cb  in  itself  is  not  im- 
portant; what  is  important  is  the  actual  direction  of  the 
displacements  associated  with  a given  Q . This  can  be  de- 
termined only  in  light  of  the  previous  history  of  the  coordi- 
nate definitions  used  to  define  Q..  In  the  next  chapter,  we 
shall  dwell  at  length  on  the  signs  of  dipole  moment  deriva- 
tives but  keep  in  mind  that  the  signs  have  meaning  only  with 
reference  to  the  defined  coordinate  axis,  the  internal  coordi- 
nate definitions,  and  L. 


CHAPTER  3 


POLAR  TENSORS  FOR  THE  METHYL  HALIDES 
Introduction 

In  this  chapter  we  present  the  atomic  polar  tensors  for 
the  hydrogen,  carbon,  and  halogen  atoms  derived  from  the  ex- 
perimental intensity  data  for  CH^X  and  CD^X  molecules  (X  = 

F,  Br,  Cl,  I)  for  all  the  possible  choices  for  the  signs  of 
the  3p/3(L  's,  together  with  the  polar  tensors  for  CH^F  and 
CH^Cl  calculated  by  the  CNDO/2  method.  We  then  extend  the 
error  propagation  treatment  given  by  RNO  (18)  to  the  polar 
tensors  and  follow  by  a discussion  of  our  best  estimate  of 
the  signs  of  the  3p/3Cb  's. 

We  shall  restrict  our  attention  exclusively  to  the 

a 

rotationally  corrected  tensors  P since,  as  Crawford  has 

— X 

explained  (11) , rotational  corrections  are  required  before 

comparisons  can  be  made  between  such  intensity  parameters 
a 

(P  , or  3p/3R.'s)  from  different  isotopically  substituted 
—x  3 

molecules.  Our  neglect  of  the  non-rotationally  corrected 

a a 

tensors  (V  ) is  not  meant  to  imply  that  the  V 's  are  not 

useful  for  interpretative  purposes.  It  is  just  that  we  wish 

to  use  the  isotopic  invariance  to  help  in  the  determination 

of  the  correct  signs  of  the  3p/3Q  's.  We  also  neglect  the 

i 

->■ 

Pg  tensors  (3p/3S  ' s)  since  RNO  have  dealt  with  them  in 
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detail  (18) . However,  both  the  Va  and  P0  tensors  for  the 

— — O 

methyl  halide  series  have  been  computed  but  are  not  reported 
in  this  thesis. 

In  order  to  insure  consistency  in  the  definition  of 
coordinates,  we  repeated  the  normal  coordinate  analyses  for 
the  entire  CH^X  and  CD^X  series  using  the  more  recent,  and 
apparently  more  well  defined,  force  field  of  Duncan  (20) . 

A complete  summary  of  the  results  from  the  normal  coordinate 
analyses  is  given  in  Appendix  I. 

Calculation  of  the  Polar  Tensors  from  Experimental  Data 

The  integrated  infrared  band  intensities  for  all  the 
CH^X  and  CD^X  molecules  were  taken  from  the  summary  of  RNO 
(18)  and  the  corresponding  values  of  | 3p/3Q . | ' s , calculated 
by  Equation  1-4,  are  given  in  Table  3-1.  The  P tensors  are 

O 

calculated  from  Table  3-1  and  Equation  2-18,  using  the  L ^ 
matrix  given  in  Table  AI-1  in  Appendix  I.  From  these,  V 

— A 

is  calculated  using  Equation  2-12.  Applying  the  rotational 
corrections  to  V as  shown  in  Equation  2-23,  we  calculated 

X. 

OL 

the  P tensors  as  shown  in  Equation  2-26.  The  data  required 
for  the  calculation  of  the  rotational  corrections  are  given 
in  Tables  3-2  and  3-3,  and  a diagram  of  the  Cartesian  coordi- 
nate axis  used  in  the  calculation  of  Pv  is  given  in  Figure 

2V 

3-1.  The  permanent  dipole  moment  is  shown  as  a vector 

defined  to  be  positive  when  it  is  pointing  from  - to  +, 

parallel  to  the  C-X  bond  axis,  as  shown.  The  calculation 

of  P from  EA  was  carried  out  in  one  calculation  on  the 
X Q 

computer,  using  the  program  PVDTEN  described  in  Appendix  II. 
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TABLE  3-2 

CENTER  OF  MASS  CARTESIAN  COORDINATES 
OF  THE  METHYL  HALIDES  (UNITS  ARE  ®)a 


CH  F 
3 


Atom  # 
1 
2 

3 

4 

5 


(e) 

x 

0.000000  0 

0.000000  0 

-1.038888  0 

0.519444  -0 

0.519444  0 


(e) 

(e) 

y 

z 

000000 

-0.740982 

000000 

0.641018 

000000 

-1.087010 

899703 

-1.087010 

899703 

-1.087010 

1 

2 

CD  F 3 

3 

4 

5 


0.000000 

0.000000 

-1.038886 

0.519444 

0.519444 


0.000000 

0.000000 

0.000000 

-0.899703 

0.899703 


-0.652389 

0.729611 

-0.998420 

-0.998420 

-0.998420 


1 

2 

CH3C1  3 

4 

5 


0.000000 

0.000000 

-1.030521 

0.515262 

0.515262 


0.000000 

0.000000 

0.000000 

-0.892459 

0.892459 


-1.223346 

0.554654 

-1.559617 

-1.559617 

-1.559617 


Continued  to  next  page 
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(Table  3-2  continued) 


Atom 

#b  x(e> 

(e) 

y 

(e) 

z 

1 

0.000000 

0.000000 

-1.134529 

2 

0.000000 

0.000000 

0.643471 

"D3C1 

3 

-1.030521 

0.000000 

-1.470800 

4 

0.515262 

-0.892459 

-1.470800 

5 

0.515262 

0.892459 

-1.470800 

1 

0.000000 

0.000000 

-1.615035 

2 

0.000000 

0.000000 

0.319946 

:K3Br 

3 

-1.033611 

0.000000 

-1.941685 

4 

0.516807 

-0.895135 

-1.941685 

5 

0.516807 

0.895135 

-1.941685 

1 

0.000000 

0.000000 

-1.554581 

2 

0.000000 

0.000000 

0.380419 

3D3Br 

3 

-1.033611 

0.000000 

-1.881230 

4 

0.516807 

-0.895135 

-1.881230 

5 

0.516807 

0.895135 

<r 

-1.881230 

1 

0.000000 

0.000000 

-1.903073 

2 

0.000000 

0.000000 

0.232926 

3H3I 

3 

-1.034705 

0.000000 

-2.222871 

4 

0.517353 

-0.896082 

-2.222871 

5 

0.517353 

0.896082 

-2.222871 

Continued  to  next  page 
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(Table  3-2  continued) 


b 

Atom  # 

(e) 

X 

(e) 

y 

(e) 

z 

1 

0.000000 

0.000000 

-1.856776 

2 

0.000000 

0.000000 

0.279223 

3 

-1.034705 

0.000000 

-2.176575 

4 

0.517353 

-0.896082 

-2.176575 

5 

0.517353 

0.896082 

-2.176575 

a.  Calculated  by  program  CART  (J.H.  Schachtschneider , 
"Vibrational  Analysis  of  Polyatomic  Molecules."  V. 
Technical  Report  No.  231-64,  Shell  Development  Co., 
Emeryville,  California  (1964),  pp.  3-11).  The  center 
of  mass  for  the  isotopically  substituted  CH3X  molecules 
will  shift  on  the  2 axis  due  to  the  difference  in  the 
masses  of  H and  D. 

b.  Atom  # refers  to  the  atom  sequence  in  Figure  3-1,  as 
defined  in  Figure  AI-1. 
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Figure  3-1:  Definition  of  the  coordinate  axis 

used  for  the  calculation  of  the 
polar  tensors. 

Atom  1 is  carbon,  atom  2 is  halogen 
and  atoms  3 , 4 , and  5 are  hydrogens 
(See  Figure  AI-1.) 
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9 

In  general  there  are  2 possible  sign  combinations  for 

the  8p/9Q^'s  for  each  methyl  halide  (if  there  were  no 

constraints).  However,  since  the  (3p/3Q.)  's  and  (3p/3Q.)  's 

1 y 

belong  to  the  degenerate  E class,  the  sign  combination  for 
these  two  rows  must  be  identical.  This  reduced  the  total 

g 

possible  number  of  sign  combinations  to  2 . In  addition, 

for  CH^F  and  CD^Cl , DiLauro  and  Mills  (15)  have  shown  that 

the  signs  of  the  3p/9Qi's  for  and  , as  well  as  and 

vc  , must  be  the  same.  This  reduced  the  total  number  of 
6 4 

possible  sign  combinations  to  2 . Although  there  are  no 

experimental  Coriolis  data  to  indicate  the  relative  signs 

for  CH^Br  or  CH^I,  we  shall  show  later  in  this  chapter  that 

it  is  reasonable  that  the  same  relative  signs  are  correct 

for  them. 

The  experimental  polar  tensors  for  the  CH^X  and  CD^X 
molecules  are  given  in  Tables  3-4  to  3-11  for  all  possible 
choices  for  the  signs  of  the  Sp/SQ^'s  consistent  with  the 
Coriolis  data.  The  signs  in  the  first  column  of  the  tables 
are  arranged  in  order  of  the  symmetry  blocks.  The  first  row 
of  signs  represent  the  chosen  signs  of  the  three  3p/3Q^'s  for 
the  E degenerate  vibrations;  the  second  row  of  signs  repre- 
sent  the  chosen  signs  of  the  three  3p/3Q^'s  for  the  E^  degen- 
erate vibrations;  and  the  third  row  of  signs  is  the  sign  choice 
for  the  three  3p/3Q.'s  of  the  A,  vibrations.  The  sets  of 
signs  listed  in  the  tables  and  identified  by  the  letters 

(A-H)  correspond  to  the  signs  in  the  order  given  below: 

(turn  to  page  74) 
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3P/3Q4 

3P/3Q4 

3P/3Q-L 


3p/3Qc 

D 

3p/3Qc 

D 

3p/3Q2 


3p/3Q 

3p/3Q 

3p/3Q 


We  reported  only  eight  possible  sign  combinations  in 

the  tables,  and  by  considering  the  Coriolis  data  we  can 

construct  eight  additional  acceptable  sign  sets,  all  subject 

to  the  two  constraints  just  mentioned.  These  additional 

sets  will  later  be  designated  by  a double  letter  notation 

given  by:  AB,  BA,  CII,  HC,  DG,  GD,  EF , FG.  By  this  notation 

we  mean  that  the  E and  E rows  of  the  set  indicated  by  the 

x y 

first  letter  is  combined  with  the  A,  row  of  the  set  indi- 

1,  z 

cated  by  the  second  letter.  Thus,  for  the  tensor  indicated 
by  set  AB,  we  mean  the  tensor  formed  by  combining  the  E and 

X 

E rows  of  set  A with  the  An  row  of  set  B. 
y l#z 

Calculation  of  the  Polar  Tensor  for  CH  F and  CH  Cl  by  CNDO/2 
Methods  ^ ^ 

The  CNDO/2  approximate  quantum  mechanical  technique  has 

been  described  in  detail  by  Pople  and  Beveridge  (31) . The 

calculations  of  the  polar  tensors  were  made  using  the  IBM 

370/165  computer  with  Segal's  CNDO  program  from  QCPE  (32). 
a 

The  P tensors  for  CH_F  and  CH_C1  were  calculated  by 
-X  3 3 

Op/3x)a  = (Ap/Ax)  a = (P+  - P~)/2.0(x+  - x")  . (3-1) 

T T 


t = x,  y,  or  z. 

+ _ 

Here  P^  and  P^_  are  the  calculated  permanent  dipole  moment 
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components  for  the  displacement  of  the  ath  atom  along  the 
positive  and  negative  Tth  axis  on  either  side  of  the  equili- 
brium position  of  the  atom.  The  quantities  x+  and  t are 
the  actual  displacements  (0.01A  and  -0.01A,  respectively) 
along  the  Tth  axes.  Thus,  we  approximated  (3p/3x)  by 
averaging  the  permanent  dipole  moments  calculated  from  dis- 
placement on  either  side  of  the  equilibrium  position  of  the 
ath  atom.  The  atomic  displacements  were  made  relative  to 
the  experimental  geometries  as  suggested  by  Schwendeman,  who 
found  that  from  the  point  of  view  of  the  spectroscopist , 
the  experimental  equilibrium  geometry  is  acceptable  since 
the  theoretically  calculated  minimum  geometry  "...may  lead 
to  errors  in  the  spectroscopic  parameters  which  are  larger 

than  those  due  to  all  other  causes  combined"  (33) . 

a 

The  calculated  polar  tensors  (P  ) for  CH  F and  CH  Cl 

X 3 3 

are  given  at  the  bottom  of  Tables  3-4  and  3-6,  respectively. 
Error  Propagation 

In  RNO's  analysis  of  the  propagation  of  the  experimental 
errors  for  the  3p/3S  's,  the  contributions  from  the  random 

j 

errors  in  the  experimental  intensities,  force  constants, 

molecular  geometries,  and  in  the  harmonic  frequencies  were 

taken  into  account.  We  used  RNO's  estimated  errors  in  the 

3p/3Sj's  of  the  methyl  halides  as  a starting  point  for  the 

estimation  of  the  uncertainties  in  the  atomic  polar  tensors. 

a 

An  element  of  P can  be  expressed  in  terms  of  the  3p/3S.'s 

k 3 

(the  P tensor  elements)  as 
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9p/3x.  = Y (9p/9S  ) b..  , 

1 j j 


(3-2) 


where  b..  is  the  jith  element  of  the  symmetrized  B matrix 
3 x 

and  the  sum  is  over  all  symmetry  coordinates  in  a given 

symmetry  block.  Since  RNO  have  shown  that  the  errors  in 

the  molecular  geometries  are  negligible,  then  the  only  errors 

contributing  to  the  9p/3x.  are  those  associated  with  the 

1 

3p/3Si's.  We  used  the  dispersions  listed  by  RNO  as  the 

estimated  errors  in  the  9p/3S.  elements  and  calculated  the 

a 

standard  deviation,  for  each  polar  tensor  element  by 


e?  = i l 6(3p/3S.)2(b. .)2}H 
1 j 1 ix 


(3-3) 


a 

This  value,  e , is  the  quantity  we  used  as  a measure  of  the 
errors  in  the  ath  atomic  polar  tensor  element  (i) . 

Since  the  propagated  errors  are  independent  of  the  sign 
choices  for  the  3p/3Q^'s  (18),  only  one  propagation  calcula- 
tion of  the  standard  deviations  for  each  polar  tensor  is 
required  and  the  results  are  reported  in  Tables  3-4  to  3-11 

for  each  CH  X and  CD  X atomic  polar  tensor. 

3 3^ 

We  did  not  use  the  force  fields  given  by  RNO  for  the 
calculation  (see  Appendix  I)  of  normal  coordinates  of  the 
methyl  halides,  but  rather  the  force  fields  of  Duncan. 

Since  the  latter  force  field  has  smaller  dispersions  (Chapter 
1) , the  propagated  errors  in  our  L matrix  are  expected  to  be 
smaller  so  that  errors  given  in  the  polar  tensors  (and  derived 
from  the  error  in  L given  by  RNO)  may  be  slightly  overesti- 
mated. It  was  not  our  intent  to  carry  out  an  exact  error 
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analysis,  but  just  to  obtain  an  estimate  of  the  reliability 
of  the  elements  listed  for  the  atomic  polar  tensors. 

Sign  Choice  of  the  3p/3Q.'s 

Although  eight  of  the  sixteen  possible  sign  combinations 
of  the  3p/3Ch's  were  eliminated  by  RNO,  they  were  not  able 
to  make  a unique  choice  as  to  the  correct  sign  combination 
preferred  for  the  methyl  halides.  Their  criterion  for 
selecting  the  most  probable  sign  combination  was  based  on 
the  isotopic  invariance  principle  (11) . According  to  that 
criterion,  if  the  differences  in  the  3p/3S's  for  CH  X and 
for  CD^X  exceeded  the  sums  of  the  errors  in  the  3p/3S's  for 
CH^X  and  for  CD^X  for  a given  sign  combination  of  the  3p/3Q^'s, 
then  that  sign  choice  was  unacceptable.  However,  the  eight 
remaining  sign  sets  were  consistent  with  the  Coriolis  data 
(15)  discussed  earlier  in  this  chapter.  They  were  unable  to 
choose  the  best  set  from  these  eight. 

We  employed  an  additional  criterion  to  allow  rejection 
of  seven  other  sign  sets;  namely,  agreement  of  the  experi- 
mental polar  tensors  with  the  CNDO-calculated  tensors. 

For  the  criterion  of  isotopic  invariance,  we  calculated 

a 

two  statistical  parameters  a (a)  and  d for  each  atomic 

D r ri 

tensor.  The  quantity  a (a)  is  an  average  standard  deviation 
of  the  propagated  errors  in  the  ath  tensor  elements  for  the 
whole  tensor  and  is  defined  as 

a (a)  = 1/n  {?(|e“  |2  + | e“  |2)  }h  , 

j J / u J ' n 


(3-4) 


78 


where  n is  the  number  of  non-zero  elements  in  the  ath  tensor 

and  the  sum  is  over  all  non-zero  elements.  The  quantity 
a a 

£j  ^ and  ^ are  the  propagated  errors  computed  by  Equation 

3-3  for  the  jth  elements  of  the  ath  tensor  for  CD^X  and 

CH^X , respectively.  Thus,  a (a)  gives  an  estimate  of  the 

total  error  in  the  ath  atomic  tensor  for  CH  X and  CD„X  due 

— 3 3 

to  errors  in  the  experimental  intensities,  force  constants, 
geometries  and  harmonic  frequencies  based  on  RNO ' s error 
analysis . 


a 


The  second  statistical  parameter,  d is  the  root- 

H , D 

mean-squared  difference  of  the  ath  tensor  of  CD-.X  and  CH  X 

— 3 3 

defined  as 


a _ , a a 2 % 

d = 1/n  Ulpi  (D)  “ P.(H)|  } , (3-5) 

D,H  j 3 j 

a a 

where  P. (D)  and  P. (H)  are  the  jth  element  of  the  ath  tensor 
3 3 — — 

for  CD^X  and  CH^X,  respectively,  for  a positive  sign  choice, 

and  n and  j are  the  same  as  in  Equation  3-4.  The  quantity 

d gives  a measure  of  the  difference  between  the  ath  tensor 
D / ri  

from  CD  X and  that  from  CH  X.  These  two  quantities,  a (a) 
a 3 3 

and  d , are  characteristic  of  the  tensor  as  a whole  so  that 
D r H 

the  isotopic  invariance  criterion  is  based  on  a tensor-by- 
tensor comparison  rather  than  on  an  element-by-element  com- 
parison. We  picked  the  tensor  comparison  because  it  seemed 
to  be  straightforward  and  effective  for  the  intensity  param- 
eters involving  the  entire  ath  atoms.  We  have  restricted 

, . HD  X 

our  analysis  to  the  P^,  (or  P ) and  P^  tensors,  neglecting 

the  carbon  tensors,  since  the  P tensor  can  be  found  from 

-X 
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a 


the  others.  Notice  that  d should  be  zero  within  experi- 

D,H 

mental  error  measured  by  a (a),  if  the  correct  signs  of  the 


3p/3Q.'s  were  used  in  computing  the  tensors. 

For  a given  methyl  halide,  we  computed  the  sum  of  a (H) 


+ a(X),  designated  as  a , which  is  a measure  of  the  total 

experimental  error  for  that  methyl  halide  and  its  trideuter- 

H X 

ated  isotope.  We  then  computed  d + d , designated  as 

D , H D,H 

dT  for  the  same  methyl  halide  for  each  sign  set  of  3p/3Q^'s; 

if  d^  was  greater  than  twice  the  sum  of  the  average  standard 

deviation  (2a  ) for  any  sign  set,  we  considered  that  the 
T 

isotopic  invariance  principle  was  not  satisfied  and  eliminated 

that  sign  set  from  further  consideration. 

The  second  criterion  is  based  on  the  quality  of  agree- 

a a 

ment  between  the  experimental  P tensors  and  the  P tensors 

“X  -x 

calculated  by  the  CNDO  technique.  For  comparison  of  the 

experimental  tensors  to  the  CNDO  tensors,  we  calculated  a 

_H 

root-mean-squared  difference  (d  ) between  the  two,  similar 

T a 

to  Equation  3-5,  for  the  CH^X  tensors  by  replacing  the  P. (D) 

tensor  in  Equation  3-1  with  the  CNDO  calculated  tensor  element 

a a 

P.(CNDO).  In  like  manner,  we  calculated  d for  each 

“i  D ,CNDO  a 

tensor  element  of  the  CH_^X  molecule  by  replacing  the  (H) 

tensor  element  by  the  CNDO  calculated  tensor  element.  The 

H X -HD 

sums  d + d designated  as  d , and  d_  „„„  + 

x H , CNDO  H ,CNDO  _D  T D,CNDO 

d , designated  as  d , were  calculated  for  the  eight 

H , CNDO  T 

sign  sets  remaining  after  applying  the  first  criterion.  We 

chose  as  the  most  probable  sign  set  that  set  for  which  the 

-D  -H 

d and  d were  smallest. 

T T 
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Selection  of  the  Signs  of  the  3p/3Q  's  for  CH  F 
i 3— 

The  P tensors  for  CH  F and  CD  F are  given  in  Table  3-4 
-X  3 3^ 

and  3-5  for  all  possible  sign  choices  of  the  3p/3Q^'s  con- 
sistent with  the  Coriolis  data.  The  statistical  data, 

described  in  the  previous  section,  used  for  the  selection 

->■ 

of  the  sign  set  for  the  3p/3Ch's  are  given  in  Table  3-12. 

By  using  criterion  1 (isotopic  invariance) , we  eliminated 
sets  C,  D,  G,  H,  CH,  HC,  DG,  and  GD.  ENO  arrived  at  the 
same  conclusions  with  respect  to  non-acceptable  sign  combi- 
nations with,  however,  one  exception.  They  found  that  they 

could  reject  the  combination  for  the  E modes.  As  shown 

in  Table  3-12,  we  cannot  reject  this  sign  choice  (which 
appears  in  sets  A and  AB)  based  solely  on  isotopic  invariance. 

Of  the  remaining  sign  sets  (A,  B,  E,  AB,  BA,  EF , and  FE) , 

-H  -D 

set  B gives  a minimum  for  both  d and  dm  , and  therefore  is 
our  preferred  sign  set  for  CH^F. 

Restricting  our  attention  to  the  sets  not  ruled  out  by 
the  first  criterion  (sets  A,  B,  AB , BA,  EF,  and  FE) , we 
should  point  out  a fact  that  is  implicit  in  the  numbers  in 
Table  3-12,  but  perhaps  not  obvious.  If  we  examine  the  A 

1 , z 

rows  of  the  E and  F sign  sets  in  Table  3-4  (or  Table  3-5)  for 
F 

the  P tensor,  we  see  that  the  experimental  C-F  bond  moment 
~^X 

derivative  is  a relatively  large  positive  number.  This 

result  is  intuitively  unreasonable,  since  it  implies  that 

F 

the  fluorine  atom  "effective  charge"  3p  /3z  is  positive. 

z 

Previous  studies  of  fluorine  containing  compounds  (such  as 
CF4 , F2CO,  and  CgFg)  support  the  intuitive  idea  that  this  C-F 
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bond  moment  parameter  is  negative  (4).  As  a consequence,  we 
feel  that  sign  sets  E and  F should  be  eliminated  from  con- 
sideration. This  conclusion  requires  that  the  sign  sets  EF , 
and  FE,  also,  be  discarded,  since  any  sign  set  from  the  E 

x 

and  E rows  of  set  E or  F must  be  restricted  with  a set  of 


signs  from  the  A rows  consistent  with  the  Coriolis  data. 

± rz 

By  means  of  this  pragmatic  argument,  only  sets  A,  B,  AB,  and 

-H 

BA  must  still  be  considered  possible.  Examination  of  d 

-D  . T 

and  d in  Table  3-12  for  these  sets  shows  little  difference 
T 

between  the  four  sets.  If,  however,  we  examine  only  the  P 
D 

(or  P ) tensor  of  set  A (in  Table  3-4)  , we  see  that  the  first 
element  of  the  E row,  3p  /3x,  has  a relatively  large  posi- 

X X 

tive  value  although  CNDO/2  predicts  an  equally  large  negative 

value;  also  the  first  element  of  the  A row,  3p  /3x  , for 

1,  z z 

set  A is  a relatively  large  positive  number  and  CNDO/2  pre- 
dicts a negative  value  for  that  element.  However,  the  signs 

of  these  elements  from  set  B are  consistent  with  the  CND0/2- 

F 

predicted  signs.  (Note  that  the  P tensors  for  both  set  A 

~ X 

and  set  B are  very  similar  and  we  do  not  gain  any  additional 
information  by  examining  these  tensors.) 


Determination  of  the  Siqns  of  3p/3Q  for  CH  Cl 
— — — — 3 — 

The  polar  tensors  for  CH  Cl  and  CD  Cl  for  all  possible 

33. 

sign  choices  of  the  3p/3Q.'s  are  given  in  Tables  3-6  and 

3-7.  As  in  the  case  of  CH^F,  DiLauro  and  Mills  (15)  have 

shown  that  the  signs  of  the  3p/3Q. 's  for  v and  v , as  well 

l 2 3 

as  for  v and  v , must  be  identical,  and  as  a consequence, 

4 6 

CH^Cl  has  the  same  number  of  possible  sign  sets  as  CH^F. 
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Following  a similar  procedure  as  for  CH^F  (as  described 
below) , set  B is  found  to  be  the  preferred  set  of  signs 
for  CH^Cl  and  CD^Cl.  The  data  used  in  making  this  choice 
are  given  in  Table  3-13. 

We  have  used  the  same  criteria  (isotopic  invariance 

and  CNDO-calculated  tensors)  for  selecting  a preferred 

sign  set  for  CH^Cl  as  was  used  for  CH^F  with,  however,  a 

modification  of  the  CNDO  criterion.  From  Table  3-6  we  see 

that  the  CND0/2-calculated  bond  moment  parameter  3p  /3z 

z 

for  the  C-Cl  bond  is  positive.  We  feel  that  this  is  clearly 
an  error  on  the  part  of  the  CNDO/2  calculation  for  chlorine, 
based  on  the  same  reasoning  given  for  rejection  of  corre- 
sponding positive  C-F  bond  moment  parameters.  Hence,  we 

, , H , D H 

used  only  the  P and  P tensors  to  compute  d and 

D ~X  “X  H , CNDO 

dD  CNDO  ' even  thou9h  the  inability  of  CNDO/2  to  predict 
the  C-Cl  bond  moment  sign  introduces  a larger  degree  of 
uncertainty  in  the  entire  calculation. 

Using  the  isotopic  invariance  criterion,  we  eliminated 
sets  C,  D,  G,  H,  CH,  HC,  GD,  and  DG.  Of  the  remaining  sets 

(A,  B,  E,  F,  AB,  BA,  EF , and  FE) , set  B gives  the  minimum 

H n 

value  for  d and  dD 

D,CNDO  H , CNDO 

Determination  of  the  Signs  of  3p/3Qi's  for  CH^Br  and  CH3 I 

The  polar  tensors  for  CH^Br , CD^Br , CH^I,  and  CD^I  are 
given  in  Tables  3-8  to  3-11  for  all  possible  sign  combina- 
tions of  the  3p/3Q^'s.  We  will  treat  these  molecules 
collectively  since  the  problem  of  determination  of  the 

—y 

signs  of  the  3p/3Q.'s  are  common  for  both  CH_Br  and  CH  I. 

i 3 3 
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We  were  unable  to  eliminate  any  of  the  possible  sign 

sets  for  CH  Br  using  the  isotopic  criterion.  For  each  sign 
H 3 Br 

set,  d + d was  smaller  than  twice  the  sum  of  the  stan- 


H ,D  H , D 

dard  deviation  as  shown  in  Table  3-14;  however,,  a less 

H Br 

rigorous  criterion  (that  d + d be  less  than  a)  would 
eliminate  all  sets  except  E,  F,  EF,  and  FE. 

Since  we  do  not  have  a CNDO/2  calculation  of  the  polar 
tensors  for  molecules  containing  atoms  outside  the  first  and 
second  row  elements,  it  appears  that  all  sign  sets  for  CH^Br 
might  have  to  be  considered.  In  addition,  there  is  no 
direct  experimental  evidence  (for  example,  Coriolis  data) 
that  requires  9p/3Q  and  3p/3Q  , or  3p/3Q  and  3p/3Q  ,to 

2 o 3 6 

have  the  same  signs,  as  was  the  case  for  CH  Cl  and  CH  F. 

3 3 

If  this  requirement  is  removed,  then  we  should  consider  an 

additional  forty-eight  sign  sets. 

In  the  case  of  CH^I,  we  were  able  to  reduce  the  number 

of  possible  sign  sets  by  the  isotopic  invariance  criterion. 

The  same  sets  as  were  eliminated  for  CH„F  and  CH  Cl  are 

3 3 

eliminated  for  CH^I,  using  the  same  criterion.  The  statis- 
tical data  are  given  in  Table  3-15.  As  in  the  case  of  CH  Br 

3 

there  are  no  Coriolis  data  available  to  determine  whether 
the  relative  signs  of  3p/3Q  and  3p/3Qj-  or  3p/3Q  and  3p/3Q 

2 5 3 6 

are  the  same.  Without  this  information,  there  are  eight 

more  possible  sign  combinations  that  should  be  considered, 

in  addition  to  those  listed  in  Table  3-15  for  CH  I. 


It  is  apparent  that  without  the  Coriolis  data  and  CNDO/2 


data  for  CH  Br  and  CH  I we  have  a rather  large  number  of 
3 3 
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TABLE  3-14 

CRITERIA  FOR  REJECTION  OF  SIGN  SETS  FOR  CH^Br 


a 

H 

Br 

Set 

Reason  for 

d 

d 

Rejection 

D ,H 

D,H 

A 

- 

.033 

. 055 

B 

.027 

.058 

C 

.034 

.083 

D 

.027 

.068 

E 

.024 

.030 

F 

.010 

.035 

G 

.031 

.075 

H 

.036 

.077 

AB 

.028 

.058 

BA 

.032 

.054 

CH 

. 034 

.077 

HC 

.035 

.083 

DG 

.027 

.075 

d 

T 

.088 

.085 

.117 

.095 

.054 

.045 

.106 

.113 

.086 

.086 

.111 

.118 

.102 


Continued  on  next  page 
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(Table  3-14  continued) 


Set 

a 

Reason  for 

H 

d 

dBr 

d 

Rejection 

D,H 

D,H 

T 

GH 

.030 

.068 

.098 

EF 

- 

. 014 

.033 

.047 

FE 

.020 

.029 

.049 

cr(H)  = .0085 
o(Br)  = .0530 

2 x{a  (H)  + CT(Br)}  = 2a  = .123 

T 


a.  For  CH^Br  we  were  unable  to  reject  any  sets  of  signs 
for  the  3p/3Q^'s  by  the  isotopic  invariance  criteria 
(see  text) , and  the  CNDO/2  technique  has  not  been 
used  here  for  elements  outside  the  second  row  of  the 
periodic  chart. 


TABLE  3-15 


CRITERIA  FOR 

REJECTION  OF 

SIGN  SETS  FOR  CH  I 

3 

a 

H 

I 

Set  - 

Reason  for 

d 

d 

d 

Rejection 

H,D 

H,D 

T 

A 

.034 

.042 

. 076 

B 

. 034 

.041 

.075 

C 

1 

.126 

.030 

.156 

D 

1 

.119 

.028 

.147 

E 

.011 

.039 

.060 

F 

. 012 

.034 

.046 

G 

1 

. 124 

.025 

.149 

H 

1 

.121 

.036 

.157 

AB 

.035 

.040 

.075 

BA 

.033 

.042 

.075 

CH 

1 

.122 

.031 

.153 

HC 

1 

.126 

.035 

.161 

DG 

1 

.124 

.027 

.151 

Continued  on  next  page 


(Table  3-15  continued) 


Set 

Reason  for 
Rejection 

dH,D 

dH,D 

dT 

GH 

1 

.119 

.027 

.126 

EF 

.013 

.037 

.050 

FE 

- 

.010 

.036 

.046 

cr(H)  = 

. 0075 

o (I)  = 

.0511 

2{a (H) 

+ a (I) } = 2ctt  = 0.117 

a.  1 means  the  isotopic  invariance  principle. 
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possible  sign  sets  that  should  be  considered.  It  appears 
to  be  impossible  at  this  point  to  make  an  objective  choice. 
However,  it  is  quite  reasonable  to  assume  for  CH^Br  and 
CH^I  that  we  can  proceed  by  analogy  with  the  CH-jF  and  CH3C1 
sign  choice.  This  argument  is  based  on  the  fact  that  the 
methyl  halides  are  structurally  similar  and  that  intuitively 
one  would  expect  to  find  a smooth  trent  in  the  physical 
properties.  We  observed  no  anomalous  behavior  in  the 
normal  coordinates  or  the  force  fields  (Appendix  I)  for 
the  series. 

The  requirement  that  the  relative  signs  of  the  dipole 

moment  derivatives  of  v.  and  vc  , and  v~  and  v,.  be  the  same 

Id  do 

for  CH^Cl  and  CH^F  is  based  on  the  knowledge  (from  a normal 
coordinate  analysis)  of  the  absolute  sign  of  the  Coriolis 
coupling  constant,  g , and  the  sign  (determined  from  the 
Coriolis  effect)  of  the  product  _ Op/3Q0  x Zp/dQ^)  (see 
Reference  34).  If  the  product  is  positive,  then  there  is 
a positive  perturbation  between  the  v and  v modes  (v^  > v_) 
which  produces  an  enhancement  of  the  R branch  and  depletion 
of  the  P branch  of  the  low  frequency  mode  (v  ) and  the 

O 

enhancement  of  the  P branch  and  depletion  of  the  R branch 
of  the  high  frequency  mode  (v  ) . If  the  product  is  negative, 
then  there  is  a negative  perturbation  between  the  coupled 
modes  and  the  intensity  perturbation  in  the  P and  R branches 
are  just  reversed.  In  principle,  by  simple  inspection  of 
the  gas  phase  infrared  spectrum,  one  can  determine  whether 
the  perturbation  is  positive  or  negative  and  a knowledge  of 
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the  sign  of  ^ allows  us  to  infer  the  relative  signs  of  the 

derivatives.  In  practice,  however,  a complete  band  shape 

analysis  (34,  35)  as  a function  of  would  have  to  be 

R , S 

carried  out. 

Y y 

We  calculated  £ „ and  £ for  the  CH  X series  (see 

2,5  36  3 

Appendix  III)  in  order  to  see  if  the  £ 's  changed  signs  for 
any  member  of  the  series.  We  found  the  calculated  values 
of  C 's  do  not  change  signs;  moreover,  their  magnitudes  were 
surprisingly  constant  or  varied  smoothly  through  the  series. 
This  fact  is  encouraging  in  that  it  suggests  that  there  is 
no  anomalous  behavior  through  the  series  with  respect  to 
the  Coriolis  constants  for  the  modes  in  question.  Since  we 
found  no  anomalies  in  the  calculated  £ constants , we  predict 

that  there  will  be  no  anomalies  in  the  signs  of  the  dipole 

moment  derivatives  for  CH  Br  and  CH  I.  We  do  not  mean  to 

3 3 

imply  that  the  calculated  signs  for  the  £ 1 s are  proof  that 

the  relative  signs  of  the  3p/3Q,' s of  and  , and 

and  v are  really  the  same  in  CH_Br  and  CHI,  but  rather 
6 3 3 

offer  the  calculation  as  evidence  to  support  our  prediction 
of  similar  behavior  through  the  CH^X  series. 

We  believe  that  it  is  reasonable  to  assume  that  the 
same  type  of  Coriolis  perturbations  found  in  CH^F  and  CH^Cl 
will  also  occur  for  CH^Br  and  CH^I.  Based  on  this  expected 
similarity,  we  have  selected  sign  set  B,  by  analogy  with 
CH^F  and  CH^Cl,  as  our  preferred  set  of  signs  for  CH^Br  and 
CH.,1,  also. 
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Summary 

We  have  assumed  that  the  total  number  of  possible  sign 

sets  is  the  same  for  each  of  the  CH  X molecules.  As  we  have 

3 

explained,  this  assumption  is  based  on  the  idea  that  there 
should  be  a close  similarity  in  the  Coriolis  perturbations 
through  the  series;  an  argument  was  presented  to  support 
this  assumption',  based  on  calculated  £ constants.  We  were 
able  to  eliminate  several  of  the  possible  sign  sets  for 
CH^F , CH^Cl,  and  CH^I  based  on  the  isotopic  invariance  prin- 
ciple. For  the  remaining  sets  (A,  B,  E,  F,  AB,  BA,  EF , and 
FE)  we  ruled  out  sets  E,  F,  EF , FE  for  CH^F  and  CH^Cl  because 
of  their  lack  of  agreement  with  the  CNDO/2  results  (and  be- 
cause of  the  pragmatic  argument  based  on  the  sign  of  the 
C-X  bond  moments) . The  bond  moment  argument  may  also  be 
used  to  eliminate  the  same  sign  sets  for  CH^Br  and  CH^I. 

We  feel  that  one  of  the  remaining  sets  (A,  B,  AB,  and  BA)  is 
the  most  probable  sign  set  for  the  series,  and  we  chose  set 
B for  CH^F  and  CH^Cl  based  on  its  agreement  with  the  CNDO/2 
calculations.  By  analogy,  we  picked  the  same  set  to  be  the 
most  probable  sign  set  for  CH^Br  and  CH^I. 

We  believe  that  the  isotopic  invariance  criterion  is 
fairly  sound,  although  it  is  applied  subject  to  the  errors 
in  the  experimental  intensity  measurements;  for  and 
of  CH^F t these  errors  are  quite  large  (18).  The  argument 
concerning  the  relative  signs  of  3p/3Q^  and  3p/3Q^  , and 
3p/3Q^  and  3p/3Q^  for  CH^Br  and  CH  I is  intuitive,  but  we 
believe  it  is  justified.  What  is  needed  experimentally  in 
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this  regard  is  a thorough  examination  of  the  Coriolis  pertur- 
bations in  CH  Br  and  CPI  I (or  the  deuterated  isotopes)  . We 
3 3 

should  point  out  that  there  is  also  Coriolis  coupling  between 

the  v and  v modes  of  the  CH_X  series  and  we  believe  these 
14  3 

intensity  perturbations  have  not  been  investigated.  This 
information  could  provide  an  additional  check  on  our  choice 
of  sign  sets.  In  regard  to  the  CNDO/2  criterion,  we  recog- 
nize the  approximate  nature  of  the  calculations,  but  a high 
degree  of  confidence  in  its  ability  to  predict  the  signs  of 
dipole  moment  derivatives  for  molecules  made  up  of  first  row 
elements  has  been  expressed  (4).  As  a consequence,  we  have 
a.  fairly  high  degree  of  confidence  in  our  preferred  sign 

choices  for  CH_F  and  CH  Cl. 

3 3 

To  summarize,  we  collect  our  preferred  polar  tensors  for 
the  CH^X  series  in  Table  3-16.  The  tensor  elements  reported 
in  the  table  are  the  average  of  the  CH_X  and  CD  X tensor 
elements  for  sign  set  B in  Tables  3-4  to  3-11,  with  the 
difference  between  the  CH^X  and  CD^X  tensor  elements  as  a 
measure  of  the  experimental  error  in  a given  element,  report- 
ed in  the  table.  (We  should  point  out  that  the  errors 
reported  in  Table  3-16  are  generally  smaller  than  the  errors 
propagated  through  from  the  8p/3S  values.  This  might  suggest 
that  we  did  overestimate  the  errors  in  the  original  error 
analysis . ) 
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CHAPTER  4 


DISCUSSION  OF  RESULTS 


Introduction 

In  this  chapter  we  shall  examine  the  atomic  polar  ten- 
sors for  the  halogen  and  hydrogen  atoms  for  the  preferred 
sign  choice  given  in  Table  3-16.  For  the  following  dis- 
cussion, we  have  changed  the  units  of  the  tensor  elements 
from  D/A  to  the  number  of  electronic  charges  (e)  (1  D/&  = 

0.208  e_)  in  order  to  conform  more  closely  to  recommendations 
(4)  based  on  SI  units.  The  preferred  polar  tensors  given 
in  Table  3-16  are  given  again  in  Table  4-1  in  these  units. 

We  shall  be  primarily  concerned  with  the  total  dipole 
moment  change  induced  in  the  molecule  due  to  a unit  dis- 
placement, 6x,  etc. , for  the  ath  atom  parallel  to  the 
Cartesian  coordinate  axis  (x,  y,  or  z)  as  given  by  Equation 
2-32.  For  a unit  displacement  (in  A)  of  the  ath  atom  parallel 
to  the  xth  Cartesian  axis,  the  dipole  moment  change  vectors, 
(AP) a , have  the  units  of  e A. 

We  shall  first  discuss  the  P tensors  and  the  trends 

X 

X 

observed  in  the  (AP)  vectors  as  the  X atom  is  changed, 

x 

followed  by  a discussion  of  the  P tensors  and  their  associ- 

X 

H 

ated  dipole  moment  change  vectors,  (AP) 
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X 


X 

The  P„  Tensors 

X 

In  this  section  we  shall  devote  our  attention  to  the 

, , . X 

halogen  atom  polar  tensors  given  in  Table  4-1.  The  P 

-X 

• X 

tensors  are  diagonal  so  the  magnitudes  of  the  (AP)  vectors 

x 

are  just  the  magnitudes  of  the  diagonal  elements. 

For  motion  of  the  X atoms  parallel  to  the  C-X  bond  (the 

z axis),  there  is  a steady  decrease  in  the  magnitude  of  ( AP ) ^ 

- * 

as  the  halogen  atom  changes  from  F to  I . Each  displacement 
, X 

vector  (AP)^  is  directed  along  the  C-X  bond  pointing  from 
the  negative  X atom  to  the  positive  carbon  atom.  The  (AP) 

z 

vectors  are  shown  plotted  to  scale  (decreasing  in  length 
^ ...  X 

from  F to  I)  m Figure  4-1.  The  (AP)  displacement  vectors 

x 

for  motion  of  the  X atom  perpendicular  to  the  C-X  bond  along 

the  x axis  are  shown  in  Figure  4-2.  The  most  striking  fea- 

x 

ture  of  the  (AP)^  vectors  is  that  there  is  so  little  change 

in  magnitude  as  X is  changed  from  F to  I . 

Because  of  the  cylindrical  symmetry,  the  dipole  vector 
X 

displacements  (AP)  for  motion  along  the  y axis  have  the 

Y X 

same  magnitudes  as  the  (AP)  vectors,  but  are  oriented  along 

x 

y.  Note  that  the  direction  of  all  the  dipole  displacement 
vectors  is  consistent  with  movement  of  a negatively  charged 
halogen  atom.  Movement  of  the  F atom,  which  is  expected  by 
electronegativity  conditions  to  have  the  highest  negative 
charge  in  a C-X  bond  is  consistent  with  this  idea  of  motion 

along  z.  (However,  the  relatively  constant  values  of  (AP) 

X X 

and  (AP)^  for  all  the  CH^X  compounds  suggests  that  this  idea 

is  too  naive . ) 
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Tensor  Rotation 

For  the  analysis  of  the  hydrogen  atom  polar  tensors, 

H 

/ it  is  more  convenient  and  informative  to  rotate  the  Car- 
tesian coordinate  system  from  that  (x,  y,  z)  shown  in  Figure 
4-1  to  one  (x!,  y z')  with  the  z'  axis  parallel  to  a C-H 
bond,  shown  in  Figure  4-3.  This  tensor  is  more  directly 
related  to  the  intensities  of  the  C-H  bond  stretching  and  C-H 
bending  motions.  The  corresponding  T'  tensor  in  the  new 
coordinate  system  has  elements  related  to  the  old  tensor  by: 


3 3 

I I a . a . T . 

i = l j = l k,l  &/J 


k,l  = 1,  2,  3. 


(4-1) 


Here  1,  2,  and  3 correspond  to  i ' , and  k'  respectively  in 

the  new  coordinate  system,  and  a . is  the  k,ith  dot  product 

k , i 

of  the  unit  vectors;  for  example,  a k i1  • j , etc.  Thus 

1,2 

the  • coefficients  are  the  cosines  of  the  angles  between 

the  original  Cartesian  axes  and  the  rotated  axes.  T is 

1 1 3 

the  i,  jth  element  of  the  original  tensors,  and  T'  is  the 

JC  f Aj 

k,  &th  element  of  the  rotated  tensor  (36). 

CL 

We  shall  designate  the  rotated  tensors  by  P , where  the 

_X ' a 

subscript  x'  indicated  the  polar  tensors  in  the  rotated  P 

X' 

coordinate  system.  We  chose  to  rotate  about  the  y axis,  so 

that  the  z'  axis  lies  along  the  C-H  bond  in  the  zx  plane, 

H 

as  shown  in  Figure  4-3.  The  resulting  P tensors  given  in 

H ~x' 

Table  4-2  were  calculated  from  the  P tensors  in  Table  4-1 

“X 

by  the  computer  program  TRANSC  described  in  Appendix  II. 
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x 


(a)  Original  Coordinate  System  ( b)  Coordinate  System  Rotated 

Through  Angle  /^About  Y 
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TABLE  4-2 

AVERAGE  OF  THE  P*? , TENSORS3  OF  CH.X  AND  CD.X 
—A  3 3 

FOR  THE  PREFERRED  SIGN  CHOICE  OF  8p/3Q'S 
(Units  are  e) 


0.012 

±.009 

0.000 

0.021 

±.005 

CH  F 
3 

0.000 

0.073 

±.009 

0.000 

-0.007 

±.015 

0.000 

-0.143 

±.005 

-0.036 

±.0004 

0.000 

0.056 

±.0006 

CH3C1 

0.000 

0.073 

±.009 

0.000 

-0.024 

±.011 

0.000 

-0.061 

±.0008 

-0.057 

±.005 

0.000 

0.057 

±.002 

CH^Br 

0.000 

0.068 

±.010 

0.000 

-0.031 

±.008 

0.000 

-0.042 

±.002 

Continued  on  next  page 


(Table  4-2  continued) 


-0.066 

0.000 

0.055 

+ .019 

+ .002 

0.000 

■ 0.060 

0.000 

+ .010 

-0.038  ' 

0.000 

-0.028 

+ .001 

±.004 

a.  As  defined  in  Figure  4-3 

b.  1 D/&  = 0.208  e. 
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The  P Tensors 
x 1 

We  shall  first  examine  the  dipole  displacement  vectors, 

H 

(AP)  , for  motion  of  the  hydrogen  atom  parallel  to  z ' . We 

z ' 

see  from  the  last  column  of  each  tensor  in  Table  4-2  that 

these  vectors  have  both  x'  and  z'  components.  This  means, 

^ H 

of  course,  that  (AP)  is  not  parallel  to  the  C-H  bond  as 

z ' 

was  required  by  symmetry  for  the  halogen  atom  tensors  in  the 
other  coordinate  system.  If  the  bond  moment  hypothesis  (as 

described  in  Chapter  1)  were  true,  we  should  expect  these 

H H 

P tensors  to  be  diagonal.  The  fact  that  P is  not 

x'  x ' 

diagonal  is  indicative  of  the  "error"  in  the  bond  moment 
hypothesis;  in  particular,  in  the  assumption  that  the  dipole 
moment  change  is  generated  in  the  direction  of  motion  of  the 
atoms . 

The  magnitudes  of  the  (AP)  vectors  are  given  in  Table 

4-3  along  with  the  angles  0 that  the  vectors  make  with  the 

z 

C-H  bond  axis.  The  positive  sense  of  9 is  defined  as  a clock 

z 

wise  rotation  from  the  C-H  bond  when  looking  down  towards  the 
ori9in  along  the  positive  y'  axis.  A pictorial  representation 
is  given  in  Figure  4-4. 

The  errors  given  in  Table  4-3  were  computed  by 

H 2 2 2 2 % 

(6P)z,=  { (3r/3P^  f + (3r/3P  ()  ap  } (4-2) 

for  the  errors  in  (AP)  and 

z ' 

6(6)  = (O0  /3P  )2  a2  + (30  /3P  ,)2  a } (4-3) 

z x'  p z z P , 

x'  ^ 

for  the  errors  in  0 . In  Equation  4-2,  r is  the  absolute 
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TABLE  4-3 


THE  (AP)  VECTORS  AND  ANGLE  0 
z 1 z 

FOR  THE  MOTION  OF  H PARALLEL  TO  THE  Z ' AXIS 


CH  F 
3 

ch3ci 

CH  Br 
3 

CH  I 
3 


(AP)Z,  (e  A) 
0.145  ± .005 

0.083  ± .001 

0.071  ± .002 

0.062  ± .003 


a 

0 (degrees) 
z 

8.2  ± 1.9 


42.9  ± 0.5 


53.5  +1.8 


62.7  ± 3.2 


a.  0 is  measured  as  a clockwise  rotation  of  (AP)  , from 
tne  C-H  bond  axis  as  seen  by  looking  down  the  positive 
y'  axis  centered  on  the  H atom. 
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X 2 2 

magnitude  of  the  (AP)  given  by  r = (P  + P ) , where  P and 

Z v i 1 v 1 


X1 

P , are  the  x''  and  z'  components  of  (AP)H,  given  in  Table  4-2; 
z z' 

c and  a are  the  experimental  errors  in  P , and  P , as  derived 
^x  pz  x z' 

from  the  difference  in  tensor  elements  of  CH3X  and  CD3X  as 

given  in  Table  4-2.  The  partial  derivatives  of  in  Equation 

4-3  were  computed  from  0 = cos  1(r/P  . 

z z 

There  are  two  trends  in  the  (AP)H,  vectors  shown  in 

z 1 

X 

Figure  4-4.  As  found  in  the  case  of  the  (AP)  vectors,  there 

z 

is  a smooth  decrease  in  the  magnitude  of  the  (AP)H,  vectors 

z 

when  the  halogen  atoms  are  changed  through  the  series  from  F 

to  I.  In  addition,  the  angle  0 increases  through  the  series. 

z 

An  examination  of  the  value  of  0^  demonstrates  quite  dramat- 
ically the  error  involved  in  the  bond  moment  hypothesis 

assumption  that  (AP)  , should  be  along  the  C-H  bond.  Figure 

z 

4-4  implies  that  the  "error"  increases  as  the  polarizability 
of  the  adjacent  substituent (in  this  case,  the  halogen  atom) 
increases,  or  as  its  electronegativity  decreases.  Note  in 

TJ 

Figure  4-4  that  the  direction  of  (AP)^,  implies  a negative 
effective  charge  on  H for  this  6z'  displacement. 

The  (AP)  , vectors  in  Table  4-2  for  the  motion  of  H 

y 

parallel  to  y'  are  allowed  by  symmetry  to  have  only  one  com- 
ponent, P^ , , parallel  to  the  y'  axis  and  pointing  in  the 

positive  y'  direction,  suggesting  a positive  effective  charge 

on  H for  this  6y'  displacement.  The  (AP)  , are  pictured  in 
y ' 

Figure  4-5,  and  are  the  same  within  experimental  error  for 
this  "out-of-plane"  motion  of  the  hydrogen  in  all  four  CH^X 
compounds. 
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H 

The  (AP)  vectors  for  the  motion  of  the  hydrogen  atom 

x' 

perpendicular  to  the  C-H  bond  and  parallel  to  the  x'  axis 

pp 

also  have  non-zero  x'  and  z’  components.  Like  the  (AP) 

z ' 

pp 

vectors,  the  (AP)  vectors  are  not  parallel  to  the  coordi- 
nate axis  and  indicate  a negative  effective  charge  on  H 

for  this  motion.  The  vectors  and  the  angle  0 {defined  as 

x 

pj 

in  the  (AP)^(  case)  are  given  in  Table  4-4.  A diagram  of 


the  vectors  is  shown  in  Figure  4-6 

H 
x 


The  large  error  in  (AP)“(  for  CH^F  and  in  the  associated 


angle  0 makes  this  vector  essentially  undetermined;  never- 
theless,  there  is  a discernible  trend  in  the  magnitudes  of 
the  vectors  as  the  halogen  atom  is  changed.  It  is  interest- 
ing that  this  trend  is  in  the  opposite  direction  from  the 

H 

trend  found  for  the  (AP)  vectors.  As  the  polarizability 

z 1 

(or  electronegativity)  of  the  halogen  decreases,  the  magni- 

H _ H 

tude  of  (AP)  t increases;  whereas,  for  (AP)  , , the  opposite 

occurred.  Note  that  (except  for  CH^F  when  it  is  undetermined) 

the  value  of  0 is  essentially  constant  for  the  CH  X molecules, 
x 3 

Correlations  with  Electronegativity 

Since  the  electronic  character  of  the  halogen  atoms  and 
substituents  is  represented  qualitatively  by  their  electro- 
negativities, this  parameter  is  a logical  choice  for  correlat- 

Oi 

ing  with  the  (AP)  vectors.  The  trends  in  the  (AP)  and 

x t 

pp 

(AP)^_  vectors  with  changing  X shown  in  the  last  section 

correlate  quite  well  with  the  electronegativities  (EN)  of 

x 

the  halogen  atoms.  Figure  4-7  shows  a plot  of  the  (AP) 

x 

vectors  vs.  EN  for  the  halogen  atoms.  As  the  EN  of  the  X atom 
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H 

(AP) 

x 1 

CH3F  0.014  ± 

CH3C1  0.043  ± 

CH3Br  0.065  ± 

CH  I 0.076  ± 

3 


AXIS3 


(degrees) 
x 

61.5  ± 55.8 


303.3  ± 11.9 


6.6 

016  300.2  ± 7.1 


TABLE  4-4 


, % H 

THE  (AP)  , VECTORS  AND  ANGLE  6 
FOR  THE  MOTION  OF  H PARALLEL  TO  THE  X' 


(e  A) 


.011 


. 006 


006 


298.7  ± 


a.  The  errors  were  computed  as  described  previously  for 
the  (AP)  , vectors. 

b.  0 is  mealureg  as  a clockwise  rotation  from  the  CH  bond 
axis  to  (AP)  , as  seen  by  looking  down  the  positive  y' 
axis  centered  on  the  H atom. 
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X 

Figure  4-7:  A plot  of  the  magnitude  of  the  (AP)T 

vectors  against  the  electronegativity 
of  the  halogen  atoms . 

The  electronegativity  values  of  the 
halogen  atoms  were  taken  from  M.C.  Day 
and  J.  Selbin,  Theoretical  Inorganic 
Chemistry  (Rienhold  Publishing  Company, 
New  York,  1962),  p.  114. 
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increases,  the  magnitude  | (AP)  | increases,  with  quite  a 

steep  slope  (nearly  0.7).  On  the  other  hand,  the  magnitudes 
X X 

| (AP) x | and  | ( AP ) | , are  almost  independent  of  the  EN  changes 

(see  Figure  4-7) . 

By  applying  the  simple  model  for  the  dipole  moment 

change  developed  by  Bruns  and  Person  (37)  for  the  diatomic 

hydrides,  HY,  (Y'=  B,  C,  N,  O,  and  F) , we  can  attempt  to 

separate  different  contributions  to  the  (AP)X  vectors  and  to 

z 

rationalize  the  observed  trends  in  the  (AP)  vectors  with  EN. 

z 

In  order  to  do  so,  we  make  two  additional  assumptions.  The 
first  assumption  is  that  the  "out-of -plane"  (AP)X  vectors 

y 

reflect  only  the  dipole  moment  change  generated  by  moving  the 
fixed  charge  on  the  halogen  atoms  out  of  the  plane,  to  give 
(AP)^.  (Bruns  and  Person  refer  to  this  as  q^.)  This  assump- 
tion is  reasonable  since  this  vector  is  essentially  constant, 
indicating  that  the  fixed  (equilibrium)  charge  on  X is  nearly 
the  same  for  each  halogen  for  the  series;  furthermore , one 
would  not  expect  the  rehybridization  moment  to  contribute  an 
appreciable  amount  to  this  motion.  The  second  assumption  in- 
volves the  rehybridization  term  (AP)  (called  q by  Bruns 

sp  sp 

and  Person) . We  shall  assume  that  this  term  is  constant  for 

the  CH^X  series,  by  analogy  with  the  HY  series;  moreover,  we 

shall  assume  the  same  constant  value  for  (AP)  as  Bruns  and 

sp 

Person  calculated  for  the  HY  series  { (AP)  = +0.20  e}.  The 

sp 

model,  based  on  the  CNDO/2  breakdown  of  contributions  to  the 
total  dipole  moment  (14),  assumes  that  the  contributions  to 
the  dipole  moment  change  vectors,  (AP)  , come  from  three 

T 
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different  terms: 

(AP)  = (AP)“  + (AP)a  + (AP) a , (4-4) 

1 2 hyb 

CL 

where  (Ap)^  is  the  dipole  moment  change  due  to  the  equilibrium 

a 

charge  on  each  atom,  and  (AP)^  , called  by  Bruns  and  Per- 
son, is  the  dipole  moment  change  due  to  the  change  in  the 

charge  when  the  C-X  bond  is  stretched.  The  (AP)a  term  con- 

hyb 

the  contribution  to  the  dipole  moment  change  due  to  the 

rehybridization  (or  polarization)  during  bond  stretching.  For 

a 

the  first  row  elements  (AP)  consists  only  of  sp  contribu- 

hyb  — 

tions, (AP)  , but  for  second  row  elements,  both  S£  and  pd 

contributions  are  present  { (AP)  = (AP)  + (AP)  } . In 

hyb  sp  pd 

Table  4-5,  the  resulting  breakdown  into  the  different  con- 
tributions is  given.  If  the  model  is  correct,  the  observed 
X 

trend  in  (AP)z  is  due  to  the  trend  for  the  (AP)  and/or  the 
(AP)^  terms.  At  this  point,  we  cannot  further  separate  the 
( Ap)  2 (AP)pd  terms,  for  it  is  not  clear  whether  one 

dominates,  or  what  is  the  percentage  that  each  term  contri- 
butes to  the  total  { (AP)  X + (AP)X  )}  . 

2 pd 

These  results  can  be  contrasted  with  those  obtained  from 
the  HY  series  (37) . In  that  series,  the  (AP)^  vectors  varied 
as  Y is  changed  across  the  periodic  table  with  increasing 
atomic  number  (and  presumably  EN) , while  the  (AP)  vectors 
were  constant.  For  the  CH^X  series,  variation  of  X down  the 
column  of  the  periodic  table  (Z  constant  but  EN  changes) 


changes  (AP)  but  not  (AP)^ 
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TABLE  4-5 


MAGNITUDE  OF  THE  DIFFERENT  CONTRIBUTIONS 

TO  THE  HALOGEN  (AP)X  VECTORS 

z 


(Units  are  e-A) 


Molecule 


(AP) 


X 


X 

(AP)  + (AP) 

2 hyb 


X xc 

. (AP)  + (AP) 

2 pd 


CH3F 


-0.93 


-0.70 


-0.90 


ch3ci 


-0.45 


-0.22 


-0.42 


CH  Br 
3 


-0.30 


-0.07 


-0.27 


CH  I 
3 


-0.12 


+ 0.11 


-0.09 


a.  Taken  from  Table  4-1. 

X 

b.  The  values  of  (AP)^  for  the  CH3X  series  were  averaged  to 
give^(AP)^  = -0.23  e-%;  this  number  is  subtracted  from  the 
(AP) z values  to  obtain  the  number  in  this  column. 

c.  These  numbers  were  obtained  by  subtracting  the  constant 
value  of  (AP)  = + 0.20  e-R,  given  by  Bruns  and  Person  (37), 
from  numbers  in  the  second  column. 


d.  The  value  of 


(AP)^  is  expected  to  be  negligible. 
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H 

The  magnitudes  of  the  (AP)  vectors  are  plotted 

z ' 

against  EN  in  Figure  4-8.  Again  we  find  a fairly  good 

linear  relationship  between  (AP)  , (and  also  0 ) and  EN. 

z z ' 

. . H 

A similar  plot  for  the  (AP)  vectors  is  shown  in  Figure 

4-9.  Here  the  errors  are  quite  large,  but  the  same  trend  is 

H 

discernible.  The  angles  0 for  the  (AP)  . vectors  are  al- 

x'  x' 


most  independent  of  EN  within  the  uncertainties  in  0 


It 


x1 


H H 

is  interesting  to  note  that  (AP)  increases,  while  (AP) 

z x' 

decreases  by  almost  the  same  amount,  with  increasing  EN. 

The  magnitude  |(AP)yI|  is  independent  of  EN  within  ex- 
perimental error,  as  shown  in  Figure  4-10.  This  "out-of- 
plane"  motion  is  quite  insensitive  to  a change  in  the  halogen 

atom  and  appears  to  be  quite  transferrable  in  the  series. 

H 

Assuming  that  the  (AP)  vectors  measure  the  contribu- 

y' 

tion  from  motion  of  the  static  charge,  as  we  did  in  discuss- 
X 

mg  the  (AP)  vectors,  we  calculated  the  contributions  to  the 
z 

H H 

hydrogen  vectors  (AP)  , and  (AP)  . that  come  from  (AP)  and 

z x 2 

(AP) 


We  did  not  attempt  an  additional  breakdown  of 

(AP)  into  (AP)  and  (AP)  since  the  motion  of  the  hydro- 
hyb  sp  pd  J 

gen  atoms  was  not  expected  to  result  in  an  appreciable  contri- 
bution from  (AP)  . We  note  that  (AP)  is  a vector  parallel 

hyb  1 

to  the  displacement  6x'  or  5z'.  The  data  for  (AP)  and  for 

H z' 

(AP)  f are  given  in  Tables  4-6  and  4-7,  respectively. 

As  Tables  4-6  and  4-7  indicate,  the  observed  trends 

(Figures  4-8,  4-9)  are  due  to  the  change  in  the  (AP)  + 

(AP)  ^ vectors  as  the  electronegativity  changes.  If  the 

(AP)  ^ term  is  nearly  constant  down  the  series  of  halogens. 


H 

Figure  4-8:  A plot  of  the  magnitude  of  the  (AP)  , 

vectors  against  the  electronegativity 
of  the  halogen  atoms. 


See  the  reference  in  Figure  4-7  for 
the  electronegativity  data. 


Electronegativity 


Figure  4-9:  A plot  of  the  magnitudes  of  the 

(AP)^,  vectors  against  the  electro- 
negativity of  the  halogen  atoms. 


See  the  reference  in  Figure  4-7 
for  the  electronegativity  data. 


Figure  4-10:  A plot  of  the  magnitude  of  the  (AP) 
vectors  against  the  electronega- 
tivity of  the  halogen  atoms . 


H 

y' 


See  reference  in  Figure  4-7  for  the 
electronegativity  data. 
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(AP)2  must  be  causing  the  trend.  We  see  that  (AP)  for  the 
6x'  displacement  becomes  more  negative  by  nearly  a factor 
of  two  down  the  series  from  CH3F  to  CH3I,  while  (AP) 2 for  the 
displacement  becomes  less  negative  by  about  the  same 
factor. 

This  discussion  has  suggested  one  possible  approach  to 
rationalizing  the  trends  found  empirically  in  examining  the 
vectors  giving  the  dipole  moment  changes.  In-  Chapter  5,  we 
shall  return  to  this  simple  approach  and  compare  the  results 
with  contributions  calculated  by  the  CNDO/2  procedure,  and 
show  that  our  model  is,  perhaps,  too  simple. 

There  are  other  effects  that  could  contribute  to  the 

a 

observed  trends  in  the  (AP)  vectors  which  we  have  not  con- 
sidered explicitly.  For  example,  we  did  not  specifically 
include  in  our  model  the  effect  of  dipole-dipole,  charge- 
dipole,  and  charge-charge  induction  effects  on  the  dipole 
moment  changes,  and  these  effects  should  be  considered  in 
any  quantitative  explanation  of  the  observed  trends.  There 
may  or  may  not  be  a cause-and-ef feet  relationship  between 

Oi 

EN  and  the  observed  trends  in  the  (AP)  vectors,  although  our 
results  suggest  that  the  EN  of  the  halogen  is  an  important 
factor.  Keep  in  mind,  however,  that  EN  correlated  with  many 
fundamental  quantities  such  as  atomic  number,  or  electron 
density , so  that  a correlation  with  EN  may  reflect  a correla- 
tion with  a more  fundamental  and  theoretically  sound  principle 
that  is  hidden  in  the  qualitative  EN  numbers  (38) . 

c 1 

The  significance  of  the  correlation  of  the  (AP)  vectors 

T 
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with  EN  is  that  there  is  a smooth  trend  that  may  be  useful 
for  predicting  polar  tensors  of  these  atoms  in  other  molecules. 
The  insensitivity  of  the  (AP)^,  and  (AP)*,  vectors  to  EN 
certainly  suggests  that  these  quantities  may  be  quite 
transferable . 

Correlations  with  Atomic  Polarizability 

The  atomic  polarizability,  a , of  the  halogen  atoms  is 

another  parameter  that  might  be  expected  to  correlate  with 
H 

the  (AP)^  vectors.  The  atomic  polarizability  of  X may  affect 
(Ap)  because  the  moment  induced  in  X by  the  oscillating  di- 
pole from  the  motion  of  the  hydrogen  atom  may  add  a term  to 
the  intrinsic  values  of  (AP  that  is  proportional  to  a. 

This  induced  moment  is  expected  to  contribute  to  the  deviation 
of  the  direction  of  (AP)^,  or  (AP)H,  from  the  direction  of 
motion  since  the  induced  moment  (AP)^n^  is  not  parallel  to 

the  direction  of  motion,  and  (AP)H  = (AP°)H  + (AP) 

x x ind 

H 

A plot  of  (AP)  vs.  the  atomic  polarizabilities  of  the 

A — — 

halogen  atoms  is  shown  in  Figure  4-11.  The  magnitude  (AP)H, 

z 

decreases  smoothly  as  the  polarizability  of  X increases,  im- 
plying an  induced  moment  that  opposes  the  direction  of  the 
intrinsic  moment.  A plot  of  9 vs.  the  atomic  polarizabili- 
ties  is  shown  in  Figure  4-12;  as  the  polarizability  of  X in- 
creases, 0 is  found  to  increase  smoothly  in  qualitative 

agreement  with  the  preceding. 

H 

A plot  of  (AP)  , vs.  the  atomic  polarizability  is  given 
in  Figure  4-13.  Within  experimental  error,  there  is  a linear 


H 

Figure  4-11:  A plot  of  the  magnitude  of  the  (AP) 
vectors  against  the  atomic  z 

polarizabilities  of  the  halogen  atoms 


The  atomic  polarizabilities  of  the 
halogen  atoms  were  taken  from 
Reference  39. 
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0 


Figure  4-12:  A plot  of  the  angle  0Z  against  atomic 
polarizability  of  the  halogen  atoms 
for  displacement  of  H along  z ' . 


See  Reference  39  for  atomic  polarizability 
data . 
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Oy  (degrees) 


a x SO"  (ern°) 


H 

Figure  4-13:  A plot  of  the  magnitude  of  the  (AP) 
vectors  against  the  atomic  x 

polarizability  of  the  halogen  atoms. 


See  Reference  39  for  atomic 
polarizability  data. 
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H 

correlation  between  (AP)  , and  the  polarizability  of  the 
halogen  atoms,  in  contrast  with  the  curve  found  for  (AP) 


H 

z ' 


There  is  no  change  in  0 with  a change  in  the  X atom  within 

X 

experimental  error,  except  possibly  for  CH^F. 

pp 

The  (AP)  t vectors  are  parallel  to  y ' . A plot  of  the 
H ^ 

(AP)  , vectors  vs.  the  atomic  polarizabilities  of  the  X atoms 

pp 

is  given  in  Figure  4-14.  The  (AP)  , vectors  are  constant 

y 

within  experimental  error  and  independent  of  changes  in  the 
polarizability  of  the  X atom. 

In  order  to  determine  whether  the  contribution  from  the 
induced  moments  in  the  halogen  atoms  are  expected  to  be  as 
large  as  required  to  explain  these  results  (see  discussion 
above) , we  carried  out  the  following  simple  calculation.  We 
assumed  that  an  oscillating  point  dipole  with  a magnitude 

pp 

and  direction  equal  to  that  of  the  z'  component  of  (AP)  , 

z 

pp 

(and  to  that  of  the  x'  component  for  (AP)  ,)  was  located  at 

X 

the  hydrogen  atom.  Using  the  atomic  polarizability,  a , of 
the  halogen  atoms  (39) , we  calculated  the  polar  coordinate 
components  of  the  dipole  induced  at  the  halogen  atom,  y^_  and 
y^  (40)  using 


and 


+ H 3 

y-  = 2a  P (cos$/r  ) 


-*  H 3 

y = aP  (sin  $/r  ) 


(4-5) 


(4-6) 


H 

Here  P^_  , the  components  of  the  (AP)  vectors,  are  co-linear 
with  the  direction  of  motion  of  the  H atom  along  the  x axis; 
(x  = z ' or  x 1 ) ; r is  the  distance  from  the  point  dipole  to 


Figure  4-14:  A plot  of  the  magnitudes  of  the  (AP)  ( 
vectors  against  the  atomic  Y 

polarizability  of  the  halogen  atoms. 


See  Reference  39  for  the  atomic 
polarizability  data. 
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the  polarizable  halogen  atom  and  $ is  the  angle  between  the 
axis  of  the  dipole  and  the  r vector,  as  defined  in  Figure 
4-15. 

The  component  of  the  induced  moment  is  positive  in 

the  same  sense  as  r,  and  the  y^  component  is  perpendicular 

to  r and  has  the  same  directional  sense  as  $ , defined  in 

Figure  4-15.  The  resulting  induced  moment  y^  is  just  the 

vector  sum  of  the  two  components.  With  a knowledge  of  the 

molecular  geometry  and  of  these  calculated  values  of  y and 

- 

y$  , yi  can  be  transformed  back  to  the  rotated  Cartesian 

coordinate  axes  (x',  y'  z')  to  find  the  contribution  from 

->  h -> 

y.  to  (Ap)  . The  data  required  to  calculate  y for  each 
1 z ' 1 

methyl  halide  are  given  in  Table  4-8,  as  are  the  calculated 

• “V 

induced  components  y and  y,.  These  values  are  converted 

r - 

to  the  resultant  induced  moment  y and  the  calculated  con- 

H 1 

tributions  to  (Ap)  , are  compared  with  experimental  values. 

z 

We  see  in  Table  4-8  that  the  dipole  induced  on  the 

halogen  atom  due  to  the  dipole  in  the  C-H  bond  generated 

when  the  H atom  is  displaced  by  6z'  is  calculated  to  have  a 

qualitative  resemblance  to  the  experimental  behavior  of 

(Ap  ( ) . This  can  be  seen  by  comparing  the  calculated  x' 
z 

components  of  y . with  the  experimental  values  of  (Ap  , ) , 

1 Z X 

shown  in  Table  4-8  in  parentheses.  We  see  that  this  esti- 
mated induced  moment,  if  added  vectorially  to  an  intrinsic 
-*0  a 

moment  (Ap  )T directed  along  the  C-H  bond,  gives  a resultant 
H - ->q  oi 

value  of  (Ap)  > {=  (Ap  )^_  + y^}  that  behaves  qualitatively 
like  the  observed  value,  in  terms  of  its  trend  with  changing 


Figure  4-15:  A schematic  diagram  of  the  model  used 

for  calculating  the  induced  moments  in 
the  halogen  atoms  assuming  an  oscillating 
unit  dipole  at  the  hydrogen  atom  parallel 
to  the  C-H  bond. 


150 


Oscillating  Unit  Dipole 


151 


— 

as 

— 

• 

^ — % 

. 

x — *. 

N 

p 

00  CN 

r-  <j\ 

oo  in 

in  i" 

- P 

CD 

OS 

• • 

Sh  3 

33 

ro  oo 

00  CN 

«o  ro 

O CN 

O -rH 



** 

i— l vr 

i — 1 LO 

CN  VO 

p p 

" — 

' — " 

" — 

V 0 

as  o 

> H 

• 

ro 

r~- 

rH 

CN 

o 

p 

• 

• 

• 

• 

•H  3 

>- 

(1) 

VO 

VO 

00 

cn 

+ 3-  as 

33 

CN 

i — 1 

rH 

p 

as  £ 

p 

- 

P 33 

N 

in 

LO 

ro 

CN 

3 

o 

o 

o 

o 

33  0 

- 

o 

o 

o 

o 

3 p 

E N 

• 

• 

• 

• 

• 

3 

• 

P 

as 

o 

o 

o 

o 

E 

00 

<3 



1 

1 

1 

1 

r->  | 

1 

' ' 

- U 

• N 

cn  ' — ’ as 

as 

p 

- 

33 

i P 

rH 

w 

X 

H ID  4J 

• p 

u - 

O H 

CD  CO 

r- 

rH  LO 

cd  as 

E N 

- 

<2 

rH  CN 

pH  LO 

rH  LD 

i — 1 U0 

x g 

1 Eh 

P 

E N 

• 

O O 

o o 

o o 

o o 

as  3 o 

S3  O 

P 

OS 

rH  H 

g 

H E-i 

< 

" 

o o 

o o 

o o 

o o 

P 33  P 

as  o 

— 

+ 

+ 

+ 

+ 

ITS  3 

i — 1 ^H 

W P 

— 

■ — 

Eh  0 3 

p p 

E H 

P O 

as 

Eh  P 

3 -H 

Eh  - 

— 

^ — x 

1 — 1 

t" 

1 — 1 

•HEP 

N 

E 2 

•H 

«2 

rH 

rH 

i — 1 

1 — 1 

1 ITS 

g ® 

O « 

+ 3- 

• 

o 

o 

o 

o 

CUP 

0 

2 

0) 

• 

• 

• 

• 

as  o 

Sh  p 

E P 

— 

o 

o 

o 

o 

> as  sh 

p o 

o 

•h  p 

H E 

pp  4H 

us  us 

EH 

0 

as  as 

2 P 

US  3 

3 3 

CO 

P O 

— 



CO 

r- 

VO 

its  as  3 

rH  rH 

1 

D 

•e 

*2 

o 

o 

o 

o 

as  o 

as  aS 

U S3 

+ 3- 

• 

o 

o 

o 

o 

as  £ -h 

• > > 

P O 

d) 

• 

• 

• 

• 

p p p 

g 

w 

2 H 

" 

o 

o 

o 

o 

3 as  o 

O » — 1 rH 

p 

U EH 

u p as 

P as  as 

m 

O 

as  Sh 

as  P P 

2 

w g 

h as  -H 

3 3 

&H 

E 

O iH  35 

E as  as 

Eh  w 

— 

o 

in 

CN 

o 

g p 

g g 

E 

Sh 

e2 

rH 

l — 1 

1 — 1 

rH 

3 as 

as  -H  -H 

E EH 

t3- 

• 

o 

o 

o 

o 

as  3 us  p Sh  3 

O , 

as 

• 

• 

• 

• 

P -h 

P as  as 

P E 

• — ' 

o 

o 

o 

o 

p as  £ 

0.  a. 

o 

p p 

3 X X 

CO  P 

P P o 

o as  as 

Eh 

o o 

P X 

us  h 35  as  as 

P 

• 

ID 

CTi 

cr» 

o 

• >i  m o 

as  u h 

CO  E 

p 

p 

• 

• 

• 

• 

cn  Sh 

3 03  3 

M H 

o 

as 

o 

in 

f" 

o 

oo  p 33  3 

as 

P 

33 

in 

as  as 

P us  us 

EH 

as  g c p 

3 as  as 

P S3 

[SOHO 

as  us  us 

3 as  p 

o as  as 

<[ J ^ 

as  p as  as 

p p 

O 

sh  ts  us 

US  P P 

2 g 

p 

. 

CN 

in 

p 

CD 

as  g 3 

•H  3 3 

Eh 

Sh 

C2 

O 

CO 

CD 

p o us  as 

x as  as 

2 

• 

t 

• 

• 

as  Sh  -H  us 

as  3 H 

P 

CN 

CN 

CN 

CN 

E p 

3 3 

>-  as 

- CU  ft 

co 

g 3 3 P 

g 

o as  as  p 

3 3 

o 

00 

CO 

rH 

Si  PH 

3 -H  -H 

co 

CN 

00 

rH 

p nS  P us  p 

3 

• 

• 

• 

• 

H C ITS  P US  US 

a 

CN 

o 

CN 

oo 

in 

a 3 as  p 

OS  3 

o 

as  o 

3 3 3 

1 — 1 

ASH  013 

^ I — 1 rH 

X 

nS  as  p a 

0 3 3 

rH 

U 

Eh  U Eh  as  rP  > > 

Em 

u 

a 

H 

oo 

ro 

ro 

00 

E 

E 

E 

E 

• • • 

• • 

U 

U 

U 

U 

as  p a 

a3  as 

152 


X.  This  trend  is  shown  in  Figure  4—16,  which  pictures  the 

vector  addition,  showing  the  apparent  increase  in  0 with 

z ' 

decreasing  EN  as  the  magnitude  (and  direction)  of  y.  changes 

relative  to  the  intrinsic  moment  (AP^) • 

We  note  that  the  magnitude  of  the  z ' component  of  the 

induced  moment  is  estimated  to  be  so  small  as  to  be  negligible, 

compared  to  the  z'  component  of  the  intrinsic  moments. 

In  conclusion,  the  induced  dipole  may  contribute  some- 

what  to  the  apparent  "rotation"  of  (AP)  , with  changing  X, 

z 

since  the  effect  of  the  induced  moment  is  qualitatively  in 
agreement  with  the  observed  behavior.  However,  the  magni- 
tude calculated  for  this  effect  is  much  smaller  than  observed, 
so  that  we  must  conclude  that  the  induced  moment  makes  only 
a very  small  contribution  to  the  observed  effect.  (See  the 

comparison  between  calculated  and  experimental  0 ■ values  ) 

z 

The  calculation  of  the  induced  moment  due  to  the  oscillat- 
ing dipole  parallel  to  the  x'  axis  and  the  associated  angles 
is  presented  in  Table  4-9.  The  coordinates  are  defined  in 
Figure  4-17.  From  the  table  we  see  that  generally  the  dipole 

induced  moment  is  smaller  for  the  (AP)  , vectors  than  for  the 

H x 
(AP)2,  vectors  in  Table  4-8.  In  the  case  of  the  x'  components 

of  1-b  ’ the  errors  in  the  experimental  values  are  about  the 

same  as  the  calculated  induced  moment  components . These 

calculated  z ' components  of  y^  are  a little  larger  than  the 

x components,  but  are  still  quite  small.  One  must  conclude 

again  that  the  effect  due  to  the  induced  moment  is  small, 

and  that  it  does  not  explain  the  observed  trends. 
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DATA  AND  RESULTS  FOR  THE  CALCULATION  OF  THE  INDUCED  MOMENT 
IN  X FOR  THE  MOTION  OF  H PARALLEL  TO  X' 
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Figure  4-17:  A schematic  diagram  of  the  model  used 

for  calculating  the  induced  moments  in 
the  halogen  atoms  assuming  an  oscillating 
unit  dipole  at  the  hydrogen  atom  perpen- 
dicular to  the  C-H  bond. 
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The  dipole  induced  moment  for  the  motion  of  H parallel 
to  y'  can  be  calculated  using  the  model  discussed  above.  In 
this  instance,  however,  there  is  only  one  non-zero  component 
of  y ; namely,  the  y$  component.  Since  $ = 90  , Equation 

4-5  becomes  zero  and  Equation  4-6  reduces  to 

= (APy-/r3)  • a * (4-7) 

The  dipole  induced  moments  are  parallel  to  the  y'  axis 

H 

and  in  the  opposite  direction  to  the  (AP)  vectors.  The 

y' 

induced  moments  are  given  in  Table  4-10.  For  CH3C1,  CH^Br , 
and  CH^I,  the  calculated  induced  moment  is  essentially 
constant,  as  Equation  4-7  suggests.  The  magnitudes  of  y ■ 
are  comparable  to  the  errors  in  (AP)  , , which  seems  to 

support  our  earlier  hypothesis  that  the  dipole  moment  change 
for  this  "out-of-plane"  motion  is  due  only  to  the  motion  of 
the  fixed  charge  on  the  hydrogen  atom. 

The  simple  calculation  of  the  induced  moment  in  the 
halogen  atom  has  qualitatively  reproduced  some  of  the  ob- 
served trends  in  the  experimental  results;  however,  the 
model  is  by  no  means  quantitative.  We  have  neglected  some 
effects  which  may  be  just  as  important  as  the  induced 
moments.  For  example,  we  completely  neglected  the  charge 
effect  of  the  hydrogen,  halogen,  and  the  methyl  carbon  in 
our  simple  approach. 

As  was  the  case  in  the  correlations  with  EN,  the  trends 
exhibited  by  the  experimental  vectors  are  the  important  fea- 
tures. The  simple  model  of  the  induced  dipole  effect  was 
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TABLE  4-10 

THE  DIPOLE  INDUCED  MOMENT  FOR  MOTION 
OF  H PARALLEL  TO  Y' 


CH  F 
3 


CH  Cl 
3 


CH^Br 


ch3i 


a 

(electrons) 

-0.003 

(0.073) 


-0.013 

(0.073) 


-0.015 

(0.068) 


-0.016 

(0.060) 


a.  The  values  in  parentheses  are  the  experimental 
(AP)^,  vectors. 
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suggested  as  one  possible  effect  that  should  contribute  to 
the  observed  trends.  The  calculations  suggest  that  it  is 
actually  of  minor  importance. 

Effective  Charges 

King,  Mast  and  Blanchette  (KMB)  (10)  first  introduced 
the  concept  of  effective  charge  of  the  ctth  atom,  , relat- 

ing it  to  the  intensity  sum  rule  of  Crawford  (11) . They 
found  for  hydrocarbons  with  symmetrically  equivalent  hydrogen 
atoms  that: 

(nH/mH^H  + (nc/mCKC  = K^Ai  + ' (4-8) 

Here  na  and  ma  are  the  number  of  atoms  of  type  a and  their 
masses,  ft  is  the  rotational  correction.  If  data  for  the 
completely  deuterated  isotopes  are  available,  a correspond- 
ing equation  can  be  written  for  the  isotopically  substituted 
species  (with  = Cq  » nH  = nj-,)  and  the  two  equations  solved 
for  the  resulting  £H  value.  By  this  procedure,  KMB  obtained 
a set  of  £h  values  from  the  various  hydrocarbons  whose  in- 
tensities had  been  studied;  the  values  of  were  surprisingly 
constant  for  the  entire  series  and  appeared  to  be  quite  trans- 
ferable for  the  molecules  for  which  data  were  available. 

For  the  CH3X  compounds,  KMB  calculated  the  effective 
charge  on  the  hydrogen  atoms  from  the  intensity  sums  for 
CH3X  and  CD^X,  as  described  for  the  hydrocarbons,  but  they 
were  not  able  to  resolve  the  effective  charge  for  the  CX  group 
into  the  atomic  effective  charges,  and  £x-  Instead,  they 
determined  an  apparent  effective  charge  for  the  CX  group: 
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2 2 2 

^cx  = (vcZc  + yxCx)/(yc  + • (4-9) 

From  the  point  of  view  of  the  atomic  polar  tensors,  the 
effective  charge  is  given  by  Equation  2—31.  Using  that  equa- 
tion, we  see  that  £^  is  an  important  invariant  of  the  atomic 

polar  tensors.  It  is  clear  that  £ = | 9P/9$  I,  so  that  £ is 

a a a 

the  magnitude  of  the  gradient  of  the  dipole  moment.  Further- 
more, this  gradient  is  given  by  the  vector  sum  of  the  APa 
vectors  defined  above: 


A?/9R 

a 


a ™ a 

AP  + AP  + AP  . 
x y z 


(4-10) 


The  effective  charges,  £^  / ?x  > an<3  for  the  atoms  in 

the  methyl  halide  molecules  are  given  in  Table  4-11,  along 

with  the  values  computed  by  KMB  from  the  intensity  sums  of 

CH^X  and  CD3X.  The  agreement  between  the  values  for  £ and 

H 

?CX  by  KMB  and  the  values  we  found  is  probably  about  as 

good  as  can  be  expected,  considering  the  rather  large  errors 
in  the  infrared  intensity  sums.  KMB  estimated  that  the  pro- 
pagated error  in  £ is  nearly  as  large  as  £ itself;  the 
propagated  errors  in  £^  are  smaller,  however,  and  nearly  the 

same  as  the  differences  here  between  the  £ values. 

H 

Examining  Tabls  4 — 11  r ws  sbb  that  thB  valuBS  of  £ foir 

H 

the  methyl  halides  from  our  analysis  are  nearly  the  same  (with- 
in experimental  error)  as  the  average  value  for  £^  in  hydro- 
carbons (£h  = 0.15  e) , thus  supporting  the  suggestion  by  KMB 
that  these  quantities  may  be  transferable.  in  order  to  ex- 
amine this  question  further,  we  include  in  Table  4-11  values 
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of  effective  charges  calculated  from  the  intensity  sums  for 
CF4  and  CCl^ . We  use  equation  4-11: 

kK  = (4.0  yx  + yc)  , (4-11) 


since  ft  = 0 for  these  non-polar  molecules.  Here  K is 
-2  2 

(1.537  x 10  /4.8)  (41),  y^  and  y^  are  the  inverse  masses 

(amu)  of  the  X and  C atoms,  and  ^A . is  the  intensity  sum 

-1  1 

(cm  mmole  ) of  the  infrared-active  fundamental  modes  of 

vibration. 

In  order  to  separate  £ and  £ in  CF  , we  recall  that 

F c 4 

for  tetrahedral  symmetry  with  z along  the  CF  bond: 


F 

P = 
X 


3p  /3x 


x 


3p  /3y 
Y F 


3p  /3z 
2 F 


Here  (3p  /3x  ) - (3p  /3y  ) , and  it  is  reasonable  to  assume 

r y r 

that  it  is  the  same  value  ( - 0.25  e)  as  found  for  the  CF  bond 


in  CH  F . We  can  estimate  (3p  /3z  ) in  several  ways  — for 
•3  z F 

example,  assume  that  5 is  the  same  in  CF  as  in  CH  F and  then 

F 4 3 

obtain  3p^/3z^  from  Equation  4-10.  Alternatively,  we  can 

obtain  3pz/3zF  from  the  value  given  by  Levin  and  Lewis  (4  2)  for 

, realizing  that  3^/3R  = 3p  /3Z  . Since  the  latter 

CF  2 F 

is  readily  available,  we  obtain  3p/3z  from  3p  /3R  = -0.984  e 

F z z — ' 

From  this  value,  and  the  value  above  for  3p  /3x  ( - 0.25 e ) 

x F — ' 
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we  calculate  £ = 1.043.  Then,  using  Equation  4-11,  we  find 

= 2*98  £.  The  values  for  CC14  are  obtained  in  a similar 

fashion,  using  the  value  of  -0.544  e (or  -0.628  e)  for 

ap/3Rccl  from  Lindsay  and  Schatz  (43) . 

From  the  results  given  in  Table  4-11  for  E from  CH  F 

F 3 

and  CF4  , and  for  ^ from  CI^Cl  and  CC14  , we  conclude  that 
?F  and  €cl  are  nearly  transferable,  as  is  E.  The  average 

effective  charge  found  by  KMB  for  the  hydrocarbons  (6  = 0 15  e) 

* - 

was  somewhat  larger  than  the  average  given  in  Table  4-11  for 
the  four  methyl  halides.  This  difference  in  values  is  only 
slightly  larger  than  experimental  error;  in  fact,  we  believe 
the  values  are  probably  the  same  within  experimental  error, 
although  it  is  possible  that  £ may  depend  somewhat  on  the 

ri 

substituent  X groups  in  x-C-H  molecules. 

Notice  that  E is  nearly  exactly  equal  to  E,  for  the 

X 

CH3X  compounds.  {We  remember  that  9p/9Rc  must  have  a sign 
opposite  (positive  here)  that  for  9p/9R  (negative  here) . } 

This  equality  implies  that  the  vector  sum  of  the  gradients 

(9p/9R^)  for  the  three  H atoms  is  zero  for  each  one  of  the 
methyl  halides. 

The  value  of  E for  the  C atom  in  CF4  , for  example,  is 
related  to  Ep  , since  (see  Equation  2-33) : 


For  the  Cartesian  coordinate  systems  we  have  been  considering 

(CF-l  bond  defines  the  z axis),  (AP)°  is  diagonal  and  the  tetra- 

x 

hedral  symmetry  requires  that  9p  /9x  = 9p  /9y  = 3p  /gz 

x C v C tz/  r ’ 
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Furthermore , 

9px/9xC  = ~ I 3px/9xF  ' (4-12) 

i i 

from  Equation  2-33.  The  polar  tensor  elements  for  the 
different  F atoms  are  related  to  each  other  by  Equation  4-1. 
For  these  tetrahedral  molecules  (44): 

'9px/3xC  = 9py/9yC  = 9pz/9zC  = 4/3{2Tll  + T33)- 

(4-13) 

Here  T^  = 3px/3xF  = T22  = 3py/3yp  and  = 3pz/9zp.  Hence, 

£2  = 3{16/9(2T11  + T33)2}  . (4-14) 

For  CF4  , Equation  4-14  yields  3.43  e and  2.34  e (or  2.56  e, 

depending  on  the  9p/9R  chosen)  for  CC14 . The  discrepancies 

between  these  numbers  and  those  in  Table  4-11  (2.98  e for  CF4 

and  2.06  e or  2.02  e for  CC14)  are  probably  within  the  errors 

due  to  the  intensity  measurements. 

Finally  it  is  of  interest  to  plot  the  effective  charges 

(Cx)  as  a function  of  the  electronegativity  of  x in  CH3X,  in 

Figure  4-18.  As  we  saw  in  the  case  of  the  (AP)X  vectors 

z 

(Figure  4-1)  , so  also  does  £ increase  linearly  with  the 
increase  in  electronegativity.  In  fact,  the  linearity  of 
this  plot  is  remarkable. 

In  order  to  test  the  possible  predictive  power  of  this 
correlation,  we  have  included  in  Figure  4-18  the  effective 
charge  for  oxygen  (?q) , calculated  from  formaldehyde  by 
Person  and  Newton  (12) , and  £ for  hydrogen  from  the  hydro- 

ri 
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carbon  series  of  KMB.  The  agreement  of  these  two  extra  points 
with  the  values  predicted  by  the  linear  correlation  is  sur- 
prisingly good,  especially  when  we  consider  that  oxygen  in 

2 

H^CO  is  double-bonded  to  an  sp_  carbon  and  is  expected  to 

have  a different  effective  electronegativity  relative  to  oxy- 

3 

gen  bonded  to  an  sp  carbon  (4  5)  . 

It  is  apparent  that  we  could  have  predicted  quite 

accurately  the  values  of  £ for  CH  Cl  and  CH  Br  if  we  had 

X d 3 

known  the  values  for  CH^F  and  CH^I.  This  ability  to  predict 
infrared  intensities  of  molecules  from  values  obtained  for 
related  molecules  is  expected  to  be  of  considerable  value  in 
future  work. 

Predicted  Intensities 

A more  general  form  of  Equation  4-8  is  (5,  12) : 

K I A + n = l £2 3 (1/m  ) . (4-15) 

i  *  1 a a a 

Here,  the  index  "i"  runs  over  all  infrared  active  fundamentals 
and  "a"  runs  over  all  atoms  in  the  molecule.  The  quantities 
K,  A, , Q,  £ , and  m have  been  previously  defined  (see  Equa- 

l a a ^ 

tion  4-8  and  Equation  4-11) . 

2 

The  intensity  sum  and  £ sum  for  the  CH  X molecules  are 

a 3 

compared  in  Table  4-12.  The  small  discrepancy  between  the 

intensity  sum  and  the  squared  effective  charge  sum  is  due  to 

the  fact  that  the  £^'s  used  for  the  computations  are  average 

values  computed  for  the  CH  X and  CD  X molecules  (taken  from 

3 3 

Table  4-11) . 

In  the  previous  section  we  predicted  the  £ value  of  CF 

C 4 


COMPARISON  OF  THE  INTENSITY  SUMS  TO  THE  SUM  OF  THE 
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taken  from  Reference  18. 
were  computed  using  the  data  in  Table 
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to  be  3.43  e (by  Equation  4-12)  from  the  values  we  found  for 
the  CH^F  fluorine  tensor  and  the  experimental  bond  moment 
of  the  C-F  bond  in  CF^. 

In  the  following  table  we  have  predicted  the  values 

of  CH.,  CPIF  and  CH  F based  on  the  values  for  CH_F  and  CF. 

4 3 2 2 3 4 

in  Table  4-11  by  plotting  B,^  as  a function  of  the  number  of 

fluorine  atoms  using  CH^F  and  CF^  as  reference  points.  The 

intercept  (to  give  the  value  of  B,  for  CH  ) is  ^ = 

C 4 C 

CHoF  CF/i 

4/3  B,  J 1/3  B,  4 = 0.274  e.  This  value  is  very  close  to 
C c 

the  experimental  value  (0.293  e)  reported  by  KMB  (see  Table 

4-11) , increasing  our  confidence  in  the  procedure. 

Using  the  £ values  given  in  Table  4-13  and  the  B and 
C H 

B,  values  calculated  for  the  CH_F  molecule  in  Table  4-11  (first 
F 3 

and  third  column,  respectively) , we  attempted  to  predict  the 

total  intensity  of  the  CH^ , CH2F2  and  CHF^  molecules  using 

Equation  4-15.  The  results  are  given  in  Table  4-14. 

The  predicted  intensity  sums  for  CHI  and  for  CH~F~  are 

4 2 2 

in  excellent  agreement  with  the  experimental  values.  The 

predicted  value  of  Ta.  for  CHF  is  somewhat  smaller  than  the 

L 1 3 

experimental  value,  but  there  is  reason  to  think  that  the 
experimental  value  reported  in  Table  4-14  is  somewhat  too 
high  (46)  . The  predicted  value  for  CF4  is  somewhat  higher 
than  the  most  recent  experimental  value  (0.971  e)  (42).  However, 
the  average  of  the  four  values  reported  in  Table  4-14  is 
1.09  + 0.14  e,  which  agrees  with  the  predicted  value  (1.19  e) 
within  the  experimental  scatter. 
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TABLE  4-13 


PREDICTED  £ VALUES  FOR  CH  , CH  F , AND  CH  F 


(Units  are  e) 


2 2 


CH 

4 

0.27  4 


b 

CH  F CH  F 

3 2 2 

1.063  1.85 


CHF 


2.64 


CF 


3.43 


a.  These  values  are  predicted  by  assuming  a linear 
relationship  between  and  the  number  of  fluorine 
atoms  in  the  molecule,  based  on  the  CH3F  and  CF. 
values  given  in  Table  4-11. 

b.  CH^F  and  CF.  are  taken  as  reference  points.  The  £ 
value  for  CH^F  was  taken  from  Table  4-11. 

The  £c  value  was  calculated  by  Equation  4-12  (see  text) . 


c . 
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TABLE  4-14 

A COMPARISON  OF  THE  PREDICTED  INTENSITY  SUM  (£a.) 

of  ch4,  ch2f2,  chf3  and  cf4 

AND  THE  EXPERIMENTAL  TOTAL  INTENSITIES 

2 

(Units  are  e ) 


*e«  (1/v 

a 

kJa  + Q 

i 

ft 

(predicted) a 

(experimental) ° 

ch4 

0.113 

0.110C 

0.0 

ch2f2 

0.444 

0 . 450d 
(0 . 405) e 

0.0066 

chf3 

0.765 

1.03f 

0.0068 

CF4 

1.19 

0.97lJ 
(1.39) J 
(1.05)  . 
(.956)  11 

o 

• 

o 

a.  The  predicted  values  are  calculated  from  Table  4-11  and 
4-13,  using  Equation  4-14. 

b.  The  data  required  for  fi  (Ix  , I ) are  based  on  the  geometry 
reported  in  "Tables  of  Inter-Atomic  Distances  and  Con- 
figurations of  Molecules  and  Ions,"  Vol.  1,  Special 
Publication  No.  11,  The  Chemical  Society,  London,  1958, 

and  the  values  reported  in  "Selected  Values  of  Electric 
Dipole  Moments  for  Molecules  in  the  Gas  Phase,"  R.D.  Nel- 
son, Jr.,  D.L.  Lide,  Jr.,  and  Moryott,  eds . , NSRDS-NBS  10, 
1967. 

c.  Experimental  intensities  taken  from  E.T.  Ruf,  M.S.  Thesis, 
University  of  Minnesota,  1959. 

d.  Experimental  data  taken  from  Reference  7. 

e.  Experimental  data  from  private  communication  (S.  Rondo  and 
S . Saeki) . 


(Continued  to  next  page) 
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(Table  4-14,  continued) 

f.  Experimental  data  from  J.  Morcillo,  J.  Herranz,  J.F. 
Biarge,  Spectrochim.  Acta,  15 , 110  (1959). 

g.  Experimental  data  taken  from  I.W.  Levin  and  T.P.  Lewis 
(Reference  42) . 

h.  B.  Schurin,  J.  Chem.  Phys . , 30_,  1 (1959). 

i.  A. A.  Chalmers  and  D.C.  McKean,  Spectrochim.  Acta,  22 , 
251  (1966)  . 

j.  R.W.  Hannah,  Ph.D.  Thesis,  Purdue  University,  1957,  as 
reported  in  Reference  47. 


175 


The  preceding  exercise  in  predicting  intensity  sums 
demonstrates  quite  well  the  transferability  of  £ values 
through  this  series  of  compounds,  and  suggests  the  usefulness 
of  the  concept. 


CHAPTER  5 


DISCUSSION  OF  THE  CNDO  RESULTS 

In  Chapter  3 we  used  the  CNDO/2  calculated  tensors  to 

help  select  the  signs  of  the  3p/3Q's.  In  this  chapter,  we 

shall  compare  in  some  detail  the  tensors,  the  associated 
a 

(AP) T vectors,  and  effective  charges  obtained  from  CNDO  with 
the  experimental  values  for  CH3F  and  CH3C1.  For  easy  reference 
we  summarize  in  Table  5-1  the  experimental  and  CNDO  tensor 
elements  for  CH3F  and  CH^Cl  that  had  been  given  previously  in 
Table  4-1. 

As  we  have  seen  previously,  the  CNDO-calculated  tensor 

elements  are  in  fairly  good  agreement  with  the  experimental 

values,  especially  for  CH^F.  There  are,  however,  notable 

F 

exceptions.  In  the  case  of  the  P tensor,  the  CNDO  values 

X 

for  3p  /3x  and  3p  /3y  agree  with  the  observed  values  almost 

x y 

within  experimental  error,  but  the  calculated  value  of  3p  /g 

c z z 

is  far  outside  the  observed  error  limit.  It  appears  that  for 
this  motion  of  the  F atom  (C-F  bond  stretching) , the  CNDO 
technique  has  underestimated  the  charge  on  the  F atom.  The 
CNDO  calculation  breaks  down  the  total  calculated  dipole 
moment  into  contributions  from  two  terms:  the  density  term 
and  the  sp  (and  £d  for  second  row  elements)  rehybridization 
(hybridization-polarization)  term  (14) . It  is  of  some  interest 
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to  examine  these  different  contributions  to  the  polar  tensors, 
in  somewhat  more  detail  than  given  previously  in  Chapter  4, 
where  we  examined  the  experimental  data.  The'  contributions  to 
the  tensor  elements  from  the  electron  density  term  and  from  the 

F 

sp  rehybridization  term  for  are  given  in  Table  5-2,  along 

ct 

with  the  contribution  from  the  static  charge  (AP^)  . Notice 

that  the  contribution  to  9pz/9z  from  the  £p  term  is  quite  small 

and  that  the  density  term  predominates.  For  9p  /3x  and  9p  /3y 

x y 

the  contribution  from  the  ££  terms  is  still  small,  but  relative- 
ly larger  than  for  3p  /9z.  It  is  apparent  from  Tables  5-1  and 
5-2  that  the  CNDO  procedure  has  underestimated  either  the  den- 
sity contribution  or  the  rehybridization  contribution  (or  both) 
for  the  3pz/3z  element,  but  has  only  very  slightly  overestimated 

one  or  the  other  (or  both)  contribution  to  the  9p  /9x  and  3p  /3y 

x y 

elements.  However,  since  the  sjo  rehybridization  term  is  so  small 

it  can  apparently  be  neglected  in  considering  the  contribution 

to  the  polar  tensors  for  the  F atom. 

A breakdown  of  the  CNDO  P , tensor  into  contributions  from 

—x1 

the  density  and  sjs  rehybridization  terms  is  given  in  Table  5-3. 

In  the  first  column  of  tensors,  the  sp_  and  density  terms  oppose 
each  other,  but  the  resultants  (Table  5-1)  are  in  fairly  good 
agreement  with  the  experimental  values  (within  the  fairly  large 
experimental  errors) . In  the  second  column,  the  s£  and  density 
terms  also  oppose  each  other;  the  resultant  CNDO-calculated 
value  is  far  outside  experimental  error.  Notice  that  in  the 
last  column  the  sp_  and  density  terms  both  have  the  same  sign, 
and  the  agreement  between  the  resultant  CNDO  tensor  elements 
and  the  experimental  values  is  fairly  good. 
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TABLE  5-2 
F 

BREAKDOWN  OF  THE  CNDO  P TENSOR  INTO 
ELECTRON  DENSITY  AND  SP_  REHYBRIDIZATION  TERMS 
(Units  Are  Electrons) 


P (density) 
—x 


. 


0.220 


0.0 


0.0 


0.0 

-0.220 

0.0 


0.0 

0.0 

-0.517 


F 

P (sp) 
“X  


-0.072 

0.0 

0.0 


0.0 

-0.072 

0.0 


0.0 

0.0 

-0.009 


(APX) 


-o. 


194 
0.0 
0.0 


0.0 

-0.194 

0.0 


0.0 

0.0 


7 


•0.194 


a.  Calculated  from  CNDO  equilibrium  charge  on  F. 
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TABLE  5-3 


BREAKDOWN  OF  THE  CNDO  P x , TENSOR  FOR  CH  F 
INTO  ELECTRON  DENSITY  AND  SP  REHYBRIDIZATION  TERMS 
(All  Units  Are  Electrons) 


1 -0.053 

O 

• 

o 

0.005 

H , 

P , (density) 

o 

• 

o 

-0.066 

o 

• 

o 

— X 

. 0.019 

o 

• 

o 

-0.076 

/ 


0.033 

o 

• 

o 

0.0077  \ 

H 

P , (sp) 
-x ' — 

0.0 

0.075 

0.0 

^ -0.045 

o 

• 

o 

-0.048 

1 0.005 

0.0 

0.0  \ 

(api  x * 

0.0 

0.005 

0.0 

o 

• 

O 

0.0 

0.005  j 

a.  Calculated  from  CNDO  equilibrium  charge  on  H 
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CH3CI 

The  contributions  from  the  density  and  from  the  sp  and 

Cl 

pd  rehybridization  terms  calculated  (CNDO)  for  the  P ten- 

X 

Cl 

sors  and  (AP)  are  given  in  Table  5-4.  In  the  3p  /3x  and 
x x 

3p^/3y  elements,  the  density  and  £p  terms  both  have  the  same 

sign,  but  the  pd  term  has  the  opposite  sign.  The  density 

contribution  is  essentially  cancelled  by  the  pd  term.  For 

the  3p  /3z  element,  the  sp  and  pd  terms  are  both  large  and 
z 

positive,  as  opposed  to  the  negative  density  term  so  that 
the  net  result  is  calculated  to  be  positive.  This  calculated 
positive  sign  is  contrary  to  our  expectations  and  is  the 
major  reason  for  our  suspicion  concerning  the  reliability 
of  the  calculation  for  this  tensor  alluded  to  earlier.  It 
appears  that  the  CNDO  procedure  grossly  overestimates  the 
contribution  for  the  sp  and  pd  terms  and/or  underestimates 
the  contribution  from  the  density  term.  We  believe  the  latter 
term  is  underestimated  because  even  if  we  neglect  the  positive 
sp  and  pd  terms  completely,  we  still  calculate  only  a relative- 
ly small  negative  density  term  contribution  to  3p_/3z. 

<u 

The  breakdown  of  the  P CNDO  calculated  tensor  for  CH_C1 

—x  3 

into  its  respective  terms  is  given  in  Table  5-5.  There  does 
not  seem  to  be  any  systematic  explanation  in  terms  of  the 
density,  S£  or  pd  terms  that  correlates  with  the  agreement 
(or  lack  of  agreement)  between  the  CNDO-calculated  tensor 
elements  and  the  experimental  tensor  elements. 

It  is  generally  conceded  that  if  the  CNDO-calculated 
values  for  the  contributions  to  the  dipole  moments  are  rela- 
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TABLE  5-4 


Cl 

BREAKDOWN  OF  THE  CNDO  P , TENSOR  INTO  ELECTRON  DENSITY 
AND  SP  AND  PD  REHYBRIDIZATION  TERMS 
(All  Units  Are  Electrons) 


/ -0.072 

o 

• 

o 

o 

• 

o 

Cl 

*x' 

(density) 

o 

• 

o 

-0.072 

o 

• 

o 

\ 0.0 

o 

• 

o 

-0.193 

-0.190 

o 

• 

o 

o 

• 

o 

PC! 

-x ' 

(sp) 

o 

• 

o 

-0.190 

o 

• 

o 

o 

• 

o 

o 

• 

o 

0.142 

0.067 

0.0 

o 

• 

o 

Cl 

P 

~x ' 

(pd) 

0.0 

0.067 

0.0 

o 

• 

o 

0.0 

0.153 

\ 


(Continued  on  next  page) 
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(Table  5-4  continued) 


a.  Calculated  from  the  CNDO  equilibrium  charge  on  Cl 
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TABLE  5-5 
H 

BREAKDOWN  OF  THE  CNDO  P_  , TENSORS  FOR  CH3C1 

INTO  CONTRIBUTIONS  FROM  THE  ELECTRON  DENSITY  AND 
SP  AND  PD  REHYBRIDIZATION  TERMS 

(Units  Are  Electrons) 


/ 

-0.089 

0.0 

0.021  \ 

H 

P (density) 

o 

• 

o 

-0.047 

0.0 

0.027 

0.0 

-0.087 

/ 0.070 

o 

• 

o 

0.001  \ 

H 

P , (sp) 

0.0 

0.078 

o 

• 

o 

\ -0.047 

0.0 

-0.030 

-0.015 

0.0 

0.012  \ 

H 

E*.  (Pd) 

o 

* 

o 

0.006 

o 

• 

o 

^ 0.038 

0.0 

-0.006  / 

1 0.012 

o 

• 

o 

o 

• 

o 

<*V*' 

0.0 

0.012 

0.0 

o 

• 

o 

0.0 

0.012J 

Calculated  from 

the  CNDO 

equilibrium  charge 

on  H. 
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tively  large,  with  opposing  signs,  any  conclusions  about  the 
sign  of  that  dipole  moment  derivative  are  even  more  suspect 
than  they  would  be  otherwise  (47) . Our  results  for  the  F 
and  CL  tensors  support  this  idea,  but  in  the  case  of  the  H 
tensor  (Tables  5-3  and  5-5)  it  is  not  clear  what  effect  these 
opposing  terms  may  actually  have. 

Recall  in  Chapter  4 that  we  broke  down  the  polar  tensors 
into  contributions  from  the  equilibrium  static  charge,  (AP)^  , 
the  charge  reorientation,  (AP)^  / and  the  term  due  to  re- 

c i 

hybridization,  (AP)^^*  In  addition,  we  assumed  that  the 

"out-of-plane"  motion  of  the  halogen  and  the  hydrogen  atoms 

reflected  the  static  charge  contribution,  (AP)^. 

Restricting  our  attention  to  the  halogen  tensors,  we 

F 

see  from  Table  5-2  that  the  calculated  (AP)^  value  is  in  good 

agreement  with  the  -0.23  e value  assumed  in  Chapter  4.  However, 

Cl 

in  Table  5-4  we  see  that  for  (AP)^  the  calculated  value  is 

-0.05  e,  whereas  in  Chapter  4 we  assumed  it  to  be  -0.23  e 

as  for  the  (AP)^  value.  Since  the  value  of  (AP)^  appears  to 

be  fairly  constant  (=0.23  e)  for  all  the  halogens,  it  appears 

Cl 

reasonable  to  assume  that  the  calculated  (AP)^  value  in 
Table  5-4  is  in  error  and  the  CNDO  procedure  has  underestimated 
the  equilibrium  static  charge  contribution  to  the  dipole 
moment  change  for  motion  of  the  chlorine  atom  in  CH^Cl. 

X 

In  Chapter  4 we  also  assumed  that  the  (AP)  term  for 

sp 

the  halogen  atoms  was  constant  at  0.20  e as  predicted  by 
Bruns  and  Person  (37)  based  on  the  analysis  of  the  diatomic 
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hydrides.  In  Tables  5-2  and  5-4,  we  see  that  the  (AP)  terms 

sp 

are  -0.009  e and  +0.142  e for  fluorine  and  chlorine,  respec- 
tively. This  discrepancy  is  not  understood,  but  may  be  due 

to  failure  of  the  simple  model  used  in  Chapter  4. 

x 

Notice  also  that  the  (AP) ^ terms  can  be  calculated  from 

x 

the  P (density)  terms  and  the  equilibrium  static  charge  term 

(AP)^  since  the  sum,  (AP) ^ and  (AP)^  / must  equal  the  total 

density  contribution.  Notice  that  in  the  case  of  the  diatomic 

hydrides,  (AP)^  was  constant,  and  (AP) 2 varied  as  Z . In  our 

case,  however,  the  calculations  seem  to  indicate  that  neither 
X X 

(AP)^  nor  (AP)2  for  the  halogens  is  constant. 

Turning  our  attention  to  the  hydrogen  tensors,  we  recall 

TJ 

that  in  Chapter  4 we  assumed  a value  of  0.070  e for  the  (AP)^ 

static  charge  contribution  for  the  CH^X  series.  In  Table  5-3 

we  see  that  the  CNDO/2  procedure  predicts  a value  of  0.005  e 

for  (AP)^  for  CH^F,  and  in  Table  5-5,  a value  of  0.012  e for 

the  same  quantity  for  CH^Cl.  The  fact  that  the  calculated 

values  of  (AP)^  are  small  is  encouraging,  even  though  the 

actual  magnitude  is  not  constant  for  CH^F  and  CH^Cl.  (Notice 

that  the  average  value  of  the  calculated  (AP)^  for  CH^F  and 

CH^Cl  is  about  +0.085  e . ) 

H H 

The  (AP)  . and  (AP)  . vectors  and  associated  angles,  0 

X Z X 

and  0z  , for  the  CNDO-calculated  tensors  are  compared  with 
the  experimental  values  in  Table  5-6.  The  comparison  is 
given  graphically  in  Figures  5-1  and  5-2  respectively.  The 
agreement  is  pretty  good,  except  for  the  obvious  discrepancy 
between  the  calculated  and  experimental  (AP)  , 
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CH^F  (see  Figure  5-1) . As  we  have  mentioned,  this  particular 

vector  is  essentially  undetermined  experimentally,  because  of 

its  large  error.  It  is  interesting  to  note,  however,  that 
, H 

the  CNDO-calculated  (AP)  , vector  is  in  the  fourth  quadrant 

X 

with  the  experimental  (AP)  , vectors  of  CH  Cl,  CH  Br  and  CH  I, 

X o j o 

and  not  in  the  third  quadrant,  as  indicated  experimentally. 

The  qualitative  agreement  between  the  calculated  and 
experimental  (AP)  , vectors  is  very  encouraging.  We  believe 

X 

this  result  supports  a conclusion  that  the  "error"  alluded  to 

earlier  in  the  CNDO  calculation  for  CH^Cl  is  confined  largely 

Cl 

to  the  chlorine  tension,  P , and  that  the  CNDO-calculated 

—x 

H 

value  for  the  hydrogen  tensor  P may  be  nearly  correct. 

—x 

Calculated  Band  Intensities 

We  shall  now  turn  our  attention  to  the  actual  band  in- 
tensities predicted  by  CNDO  for  CH  F and  CH  Cl.  By  using 

J 

a 

the  complete  CNDO  calculated  P matrix  containing  all  P 

X X 

tensors  (Equation  2-25)  and  Equation  2-27  to  2-29,  we  obtain 
the  P^  tensor  and  from  it  we  can  calculate  the  band  intensi- 
ties using  Equation  1-2  (Chapter  1) . In  our  calculation, 
we  did  not  actually  compute  the  P tensors  for  the  two  hy- 

~ X 

drogen  atoms  that  are  out  of  the  xz  plane  by  the  CNDO  pro- 
gram. Instead,  since  the  hydrogens  are  symmetrically  equiva- 
lent, we  simply  rotated  the  H(3)  tensor  to  coincide  with  the 
two  remaining  hydrogen  atoms.  This  procedure  is,  of  course, 
general  so  that  for  a molecule  containing  two  or  more  identi- 
cal atoms  that  are  symmetrically  equivalent,  one  need  only 
compute  the  polar  tensor  for  one  of  the  atoms  and  then  rotate 
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that  tensor  to  obtain  the  tensor  for  the  remaining  symmetri- 
cally equivalent  atoms.  For  fairly  symmetric  molecules,  this 
feature  may  save  a considerable  amount  of  computer  time  in  the 
quantum  mechanical  calculation.  The  tensor  rotations  were 
carried  out  by  program  TRANSC  (Appendix  III)  and  the  computa- 
tion of  Pq  was  done  by  program  PVDTEN  (Appendix  III) . 

The  results  of  the  intensity  calculations  are  given  in 
Table  5-7,  along  with  the  experimentally  observed  intensities. 
For  CH^F , the  agreement  is,  perhaps,  better  than  we  should 
expect,  considering  the  approximate  nature  of  the  CNDO  model. 
The  v^,  and  bands  are  essentially  in  quantitative 
agreement  with  the  experimental  results,  and  is  the  right 

order  of  magnitude.  The  band  is  essentially  a C-F  stretch- 
ing motion  and  as  we  have  shown,  the  8p  /3z  element  has  been 

z 

underestimated  and  that  element  probably  accounts  for  most  of 

the  intensity  of  this  band.  The  v,  and  v bands  are  well  out- 

6 

side  the  observed  error  and  have  opposite  relative  values 
from  the  experimental  results.  In  spite  of  the  obvious  dis- 
crepancies, the  agreement  is  better  than  we  had  hoped  for. 

With  the  exception  of  the  v and  v9  bands  in  CH  Cl,  CNDO 
greatly  overestimated  or  underestimated  the  intensities  in  a 

non-systematic  way.  The  v band,  as  in  the  case  of  CH  F,  is 

3 3 

the  halogen  stretching  mode  and  CNDO  has  grossly  underesti- 
mated the  intensity.  This  is  a reflection  of  the  error  in 
the  3p  /Sz  tensor  element.  We  should  point  out  that  we  did 
not  take  into  account  the  experimental  errors  in  L when 
computing  the  intensities  given  in  Table  5-7.  The  inclusion 
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of  the  errors  in  L probably  would  not  change  our  results  by 
a significant  amount.  That  is,  we  do  not  believe  that  re- 
sulting errors  in  the  calculated  intensities  would  be  out- 
side the  observed  errors  of  the  experimental  values. 

Cl 

Since  we  believe  that  the  Op  /3z)  tensor  element  cal- 

z 

culated  by  CNDO  is  wrong,  we  replaced  that  element  by  the  ex- 
perimental value  of  the  ( 3pz/9z) ^ tensor  element  from  our  pre- 
ferred sign  set  to  see  how  much  improvement  there  was  in  the 
total  intensities  of  the  A block  (it  will  affect  only  the 
z symmetry  block) . By  adjusting  the  carbon  tensor  accord- 
ingly so  that  = 0,  we  computed  the  semi-empirical  inte- 

a 

grated  intensities  of  the  A block  of  CH  Cl  reported  in 

A.  f Z 3 

Table  5-7  in  the  last  column.  This  simple  change  in  the  PC1 
C 

and  P tensors  greatly  improved  the  v,  and  intensities. 
x 13 

There  is  no  significant  improvement  for  the  symmetric  C-H  bend- 
ing mode,  , however;  in  fact,  the  resulting  calculated  in- 
tensity of  this  mode  has  a larger  deviation  from  the  experi- 
mental value  in  the  opposite  sense. 

The  significance  of  this  semi-empirical  approach  is  not 
that  we  were  able  to  improve  the  CH3C1  intensities,  but  that 
such  an  approach,  the  combination  of  the  experimental  and  theo- 
retical tensor  elements,  may  be  very  useful  in  the  quantita- 
tive prediction  of  the  infrared  spectra  of  complex  systems . 

We  believe  that  the  preceding  calculations  of  the  inte- 
grated band  intensities  readily  illustrated  the  ease  and  use- 
fulness of  the  polar  tensor  technique  for  quantum  mechanical 
calculations  of  intensities,  even  though  the  results  from 
this  preliminary  study  are  not  completely  satisfying. 
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We  carried  through  the  analysis  of  the  different  contri- 
butions to  the  CNDO-calculated  tensors  as  an  illustration  of 
its  application  and  a possible  indication  of  where  CNDO  pro- 
cedure is  overestimating  or  underestimating  effects.  Others 
have  reparameterized  the  CNDO  procedure  (47)  for  better 
prediction  of  specific  effects.  Although  this  is  an  attrac- 
tive alternative  to  our  procedure  (using  the  standard  parame- 
ter values  in  our  calculations) , we  feel  that  at  present  it 
is  not  worthwhile.  It  defeats  the  purpose  of  the  CNDO  pro- 
cedure to  reparameterize  for  every  problem.  The  CNDO  cal- 
culation provides  an  inexpensive  general  purpose  quantum- 
mechanical  calculation  of  molecular  properties  that  generally 
leads  to  fairly  reliable  predictions.  We  believe  that  the 
CNDO  calculation  should  be  calibrated  (not.  reparameterized 
indiscriminately)  by  comparison  with  more  sophisticated  quan- 
tum mechanical  calculations. 

CNDO  Calculated  Effective  Charge 

Comparison  of  the  CNDO/2  calculated  effective  charge  to 

the  experimental  effective  charge  for  CH  F and  CH  Cl  is  given 

4 3 

in  Table  5-8.  The  CNDO  £ values  given  in  the  table  were 

a 

calculated  from  the  CNDO  P tensors  given  in  Table  4-1  by 

—x 

use  of  Equation  2-31.  The  qualitative  agreement  between  the 

CNDO/2  computed  effective  charge  and  the  experimental  values 

is  fairly  good  for  both  CH3F  and  CH  Cl,  in  spite  of  the 

somewhat  low  value  calculated  (for  £cl  and  £c  for  CH^Cl)  due 

to  the  anomalous  CNDO  8p  /2z  value  for  CHLC1.  For  the  E and 

z 3 c 

? values  in  both  compounds,  the  CNDO  values  are  somewhat 
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TABLE  5-8 

COMPARISON  OF  THE  CNDO  CALCULATED  EFFECTIVE  CHARGES  TO  THE 
EXPERIMENTAL  EFFECTIVE  CHARGES  FOR  CH  F AND  CH  Cl 

(Units  are  electrons) 


| 


5 

x 

5 


c 


— 3- 

a b 

CNDO  Experimental 

0.130  (0.164  + .007) 

0.670  (1.000  + .004) 

0.884  (1.063  + .417) 


' ch3ci 

a . t 

CNDO  Experimental 

0.138  (0.119  + .003) 

0.291  (0.557  + .002) 

0.415  (0.588  + .060) 


a.  Calculated  from  the  CNDO  polar  tensor  (Table  5-1)  using 
Equation  2-31. 

From  Table  4-11. 


b . 
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underestimated,  presumably  because  of  an  underestimation 
of  the  charge  on  the  halogens  as  discussed  earlier.  It  is 
quite  possible  that  this  underestimation  may  be  a character- 
istic fault  of  the  CNDO  calculation,  but  notice  that  the  trends 
in  the  calculated  and  experimental  results  for  £ are  parallel. 
One  should  also  keep  in  mind  that  the  values  are  determined 
by  the  and  £TT  values  and  if  "errors"  are  present  in  these 
values,  £c  will  reflect  them. 

In  the  case  of  £ for  CH  F,  the  calculated  value  is 

H J 

slightly  smaller  than  the  experimental  value  for  that  com- 
pound, and  for  CH-C1  the  calculated  £ value  is  slightly 

H 

larger  than  the  experimental  value.  Notice,  however,  that 
in  both  cases,  CH  F and  CH-C1,  the  calculated  £ values  are 

J J ri 

almost  constant  and  very  close  to  the  average  experimental 
value  of  0.130  e for  the  CH^X  series.  Although  CNDO  predicts 
an  exactly  constant  value  of  £„  for  CH_F  and  CH_C1,  the 

hi  J J 

experimental  values  do  not  appear  to  be  quite  the  same  for 
the  two  compounds,  since  the  two  values  differ  by  more  than 
2.5  (the  experimental  error)  . 


APPENDIX  I 


NORMAL  COORDINATES  OF  CH  X AND  CD  X 

3 3 

As  we  indicated  in  Chapter  3,  we  have  recalculated  the 
normal  coordinates  of  the  methyl  halides  in  order  to  insure 
consistency  in  coordinate  definitions.  We  used  the  normal 
coordinate  programs  WMAT  and  CHARLY  from  the  Molecular 
Spectroscopy  Laboratory  of  the  University  of  Minnesota, 
made  available  to  us  by  Dr.  J.  Overend.  The  programs  were 
originally  written  for  the  CDC  6400  computer  and  were  re- 
vised by  Dr.  G.  Sanchez  for  use  on  the  IBM  370/165  computer 
at  the  University  of  Florida.  These  programs  have  been 
described  in  detail  elsewhere  (19)  and  will  not  be  discussed 
further  here. 

The  CH3X  molecules  belong  to  the  C symmetry  point 
group  and  have  three  infrared  active  vibrations  belonging 
to  the  totally  symmetric  A^  representation  and  three  be- 
longing to  the  doubly  degenerate  E representation.  (The 
values  of  vVs  are  from  References  20  and  21.) 

The  molecular  geometries  of  the  methyl  halides  and  the 
trideuterated  isotopes  used  in  this  work  are  given  in  Table 
AI-1.  The  standard  equilibrium  Cartesian  coordinates,  cal- 
culated by  program  CART  (see  reference  at  bottom  of  Table 
3-2)  are  given  in  Table  AI-2  and  are  based  on  the  Cartesian 


199 


200 


axis  defined  in  Figure  AI-1.  The  internal  displacement 

coordinates  are  defined  in  Figure  AI-2,  and  the  symmetry- 

displacement  coordinates  are  defined  in  Table  AI-3  in  terms 

of  the  un-normalized  U matrix.  The  symmetrized  G elements 

are  given  in  Table  AI-4  and  AI-5  for  CH^X  and  CD  X. 

3 3 

We  carried  through  the  normal  coordinate  calculations 
using  both  Overend's  (18)  and  Duncan's  (21)  force  constants 
and  we  found  that  there  was  no  significant  difference  in 
the  resulting  polar  tensors.  Since  Duncan's  force  field 
is  better  determined  (see  Chapter  1) , we  have  reported  only 
the  results  obtained  from  his  force  field.  Duncan's 
symmetrized  force  constants  are  given  in  Table  AI-6,  and 
the  resulting  normal  coordinate  transformation  matrices 
are  given  in  Tables  AI-7  and  AI-8  for  both'  the  CH^X  and 
CD^X  series.  The  L matrices  are  given  in  Tables  AI-9 
and  AI-10.  The  calculated  harmonic  frequencies  are  given 
in  Tables  AI-11  and  AI-12.  A description  of  the  infrared 
band  assignments  is  given  in  Table  AI-13. 
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Figure  AI-1 


Cartesian  Coordinate  Axes  and  Atom  Sequence  for  CH  X 
Used  in  the  Normal  Coordinate  Analysis  3 


Cartesian  Axis  Atom  Sequence 


Atoms  1,  2,  and  3 are  in  the  xz  plane. 


Atom  # 

Atom 

a 

Atomic  Mass  (AMU) 

1 

C 

12.000000 

2 

X (X=F ,C1 ,Br , I ) 

18,99840, 

78.9183, 

34,96885, 

126.9004 

3 

H(D) 

1.007825 

(2.014102) 

4 

H(D) 

1.007825 

(2.014102) 

5 

H (D) 

1.007825 

(2.014102) 

Handbook  of  Chemistry  and  Physics, 
Rubber  Company,  Cleveland,  1969. 


a. 


50th  Edition,  Chemical 
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Figure  AI-2 

Definition  of  Internal  Displacement 
Coordinates  for  Methyl  Halides 


R = 6a  (3,  1,  4) 
7 


(continued  to  next  page) 
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(Figure  AI-2,  continued) 


R = 
8 


R 


10 


68 

(3, 

1, 

2) 

68 

(4, 

1, 

2) 

66 

(5, 

1, 

2) 

a.  Taken  from  Reference  34. 

b.  6r^  ■ means  bond  stretching  coordinate  between  atoms 
i and  j . 

c.  a (i,  j,  k)  means  angle  bending  coordinate  for  angle 
formed  by  atoms  i,  j,  k.  a is  an  HCH  angle . 

d.  8 is  an  XCH  angle. 


TABLE  AI-1 


EQUILIBRIUM  MOLECULAR  GEOMETRIES  FOR 
CH  X USED  IN  THE  a 

NORMAL  COORDINATE  CALCULATION 


ch3f 

CH3C1 

CH3Br 

ch3i 


rCH(X) 


1.095 


1.084 


1.084 


1.083 


rCX  ^ 
1.382 

1.778 

1.935 

2.136 


110°30 

110°50 

111°20 

111°30 


a 


Data  taken  from  Reference  21 
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TABLE  AI-2 


STANDARD  EQUILIBRIUM  CARTESIAN  COORDINATES 
FOR  THE  METHYL  HALIDES,  CH3X 
(Units  are  in  Angstroms) 


CH  F 
3 

Atom  - 
1 
2 

3 

4 

5 


X 

y 

z 

0 

• 

0 

0.0 

0 

• 

0 

0.0 

0.0 

1.382 

-1.038888 

0 

• 

0 

-0.346030 

0.519444 

-0.899703 

-0.346030 

0.519444 

0.899703 

-0.346030 

1 

0 

• 

0 

0.0 

0 

• 

O 

2 

0 

• 

O 

0.0 

1.778 

3 

-1.030523 

0.0 

-0.336272 

4 

0.515261 

-0.892459 

-0.336272 

5 

0.515261 

0.892459 

-0.336272 

CH3Br 


1 

0 

• 

0 

0.0 

0 

• 

0 

2 

0 

• 

0 

0.0 

1.935 

3 

-1.033613 

0.0 

-0.326650 

4 

0.516806 

-0.895135 

-0.326650 

5 

0.516806 

0.895135 

-0.326650 

Continued  to  next  page 
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(TABLE  AI-2  continued) 


CH  I 


Atom 

X 

y 

z 

1 

o 

• 

o 

0.0 

o 

• 

o 

2 

o 

• 

o 

0.0 

2.136 

3 

-1.034707 

o 

• 

o 

-0.319799 

4 

0.517353 

-0.896082 

-0.319799 

5 

0.517353 

0.896082 

-0.319799 
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TABLE  AI-4 


UPPER  TRIANGLE  OF  THE  SYMMETRIZED  G MATRIX 

FOR  CH-X 

-1  a 

(all  units  are  amu  ) 


Symmetry 

Block 


ch3f 

ch3ci 

CH3Br 

ch3i 

gll 

1.0172005 

1.0162935 

1.0149364 

1.0140343 

gl2 

-0.0949003 

-0.0936552 

-0.0902867 

-0.0880854 

g!3 

-0.0456120 

-0.0447755 

-0.0434943 

-0.0426214 

g22 

1.9513454 

1.9656754 

1.9266357 

1.9035587 

g23 

0.1733833 

0.1743044 

0.1729860 

0.1722240 

g33 

0.1359693 

0.1119302 

0.0960046 

0.0912135 

g44 

1.1047525 

1.1052065 

1.1058836 

1.1063356 

945 

0.1481217 

0.1511650 

0.1535017 

0.1552210 

g46 

-0.1200393 

-0.1008428 

-0.0947300 

-0.0884736 

g55 

2.2960701 

2.3571835 

2.3790112 

2.3982687 

g56 

0.2393007 

0.2647356 

0.2627645 

0.2648623 

966 

0.9969396 

0.9480025 

0.9284539 

0.9171713 

. The  valence  bending-stretching  interaction  elements 

have  been  weighted  by  1A  and  the  valence  bending-bending 
elements  have  been  weighted  by  (1A)  . 


a 
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TABLE  AI-5 

UPPER  TRIANGLE  OF  THE  SYMMETRIZED  G MATRIX 

FOR  CD3X  ~ 

(all  units  are  amu  )a 


Symmetry 

Block 

cd3f 

cd3ci 

CD  Br 
3 

cd3i 

gll 

0.5214643 

0.5205573 

0.5192000 

0.5182982 

g!2 

-0.0949002 

-0.0936524 

-0.0902716 

-0.0880544 

A, 

gi3 

-0.0456120 

-0.0447755 

-0.0434943 

-0.0426214 

I 

g22 

1.1566544 

1.1656857 

1.1430759 

1.1295519 

g23 

0.1733833 

0.1742998 

0.1729565 

0.1721638 

g33 

0.1359693 

0.1119302 

0.0960046 

0.0912135 

g44 

0.6090165 

0.6094702 

0.6101485 

0.6105996 

g45 

0.1481215 

0.1511648 

0.1535016 

0.1552206 

E ,E 

g46 

-0.1200391 

-0.1008427 

-0.0947299 

-0.0884736 

x y 

g55 

1.2463379 

1.2805548 

1.2940044 

1.3055553 

g56 

0.0407912 

0.0650535 

0.0675585 

0.0724000 

g66 

0.5834892 

0.5261183 

0.5065695 

0.4945073 

a.  The  valence  bending-stretching  interaction  elements 

have  been  weighted  by  1A  and  the  valence  bending-bending 
elements  have  been  weighed  by  (1A)^. 
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TABLE  AI-6 

UPPER  TRIANGLE  OF  THE  SYMMETRIZED 
F MATRIX  FOR  CH  X 

3. 

(all  units  are  mdynes/X) 


Symmetry 

' CH  F 

CH  Cl 

CH  Br 

CH  I 

Block 

3 

3 

3 

3 

Hi 
M 
i— 1 

5.248 

5.493 

5.536 

5.533 

f12 

-0.222 

0.028 

0.076 

0.126 

A 

1 

0.385 

0.173 

0.121 

0.072 

f 22 

0.757 

0.636 

0.600 

0.559 

f 23 

-0.690 

-0.526 

-0.469 

-0.411 

f 33 

5.692 

3.500 

2.940 

2.394 

f 44 

5.276 

5.364 

5.436 

5.450 

f 45 

-0.132 

-0.132 

-0.146 

-0.142 

E f« 

0.176 

0.061 

0.081 

0.128 

f 55 

0.578 

0.542 

0.534 

0.524 

f 56 

-0.066 

-0.019 

-0.019 

-0.011 

f 66 

0.906 

0.718 

0.641  ' 

0.557 

a.  Taken 

from  Reference 

21. 
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TABLE  AI-7 

CALCULATED  SYMMETRIZED  L MATRIX  FOR  CH  X 
(all  units  are  amu  2) 


Symmetry 

Block 

ch3f 

CH3C1 

CH  Br 
3 

ch3i 

1.004814 

1.008100 

1.007339 

1.006582 

*12 

0.086789 

0.001071 

-0.013873 

-0.028588 

cn 

rH 

<=>? 

-0.004100 

-0.003938 

-0.003150 

-0.002689 

£21 

-0.211729 

-0.093734 

-0.070199 

-0.048175 

Ai 

Z 

22 

1.367030 

1.388852 

1.379216 

1.372646 

1 

*23 

0.194278 

0.167260 

0.139536 

0.130703 

*31 

-0.049811 

-0.043251 

-0.040972 

-0.038827 

*32 

0.068103 

0.083930 

0.093749 

0.097128 

*33 

0.358956 

0.320959 

0.292467 

0.283323 

*44 

1.050864 

1.051079 

1.051319 

1.051511 

*45 

0.012734 

0.020922 

0.023556 

0.019036 

*46 

-0.016576 

-0.001320 

-0.006915 

-0.017190 

E ,E 

*54 

0.119930 

0.113393 

0.110444 

0.117192 

x y 

*55 

1.499687 

1.514077 

1.528255 

1.535962 

*56 

-0.180620 

-0.227802 

-0.176767 

-0.159220 

*64 

-0.102633 

-0.101162 

-0.090480 

-0.074124 

*65 

0.282452 

0.319943 

0.284446 

0.272934 

*66 

0.952171 

0.914005 

0.916166 

0.914974 
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TABLE  AI-8 

CALCULATED  SYMMETRIZED  L MATRIX  FOR  CD0X 

— -U  3 

(all  units  are  amu  ) 


Symmetry 

Block 


A 


1 


cd3f 

cd3ci 

CD3Br 

CD  I 
3 

*11 

0.717525 

0.721401 

0.720552 

0.719798 

£12 

0.077600 

0.011547 

-0.001174 

-0.013565 

*13 

-0.024478 

0.002148 

0.002200 

0.002223 

*21 

-0.250852 

-0.146994 

-0.123636 

-0.102532 

*22 

0.994273 

1.069451 

1.061608 

1.057103 

0 

*23 

-0.324261 

0.018732 

0.027807 

0.039545 

*31 

-0.080539 

-0.065321 

-0.060928 

-0.057224 

*32 

0.241159 

0.148880 

0.148887 

0.147693 

*33 

0.267067 

0.292400 

0.264811 

0.257149 

*44 

0.780369 

0.780630 

0.781108 

0.781384 

*45 

-0.005721 

0.000589 

0.002940 

0.000588 

*., 

46 

0.002858 

0.009331 

0.003177 

-0.006225 

*54 

0.197614 

0.193876 

0.192573 

0.197260 

*55 

1.096870 

1.111297 

1.119186 

1.123386 

*56 

0.064468 

-0.089384 

-0.065887 

-0.068152 

*64 

-0.156403 

-0.137638 

-0.124543 

-0.107874 

*65 

0.021440 

0.138734 

0.122412 

0.124860 

*66 

0.747373 

0.698517 

0.689981 

0.683579 
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TABLE  AI-9 

CALCULATED  SYMMETRIZED  L INVERSE  MATRIX  FOR  CH  X 

' __  1^*  ^ 

(all  units  are  amu  2) 


Symmetry 

CH  F 

CH-C1 

CH-.Br 

CH  I 

Block 

L-1 

£11 

3 

0.983857 

J 

0.992372 

j 

0.993631 

3 

0.994450 

*12 

-0.064769 

-0.001550 

0.009578 

0.020720 

*13 

0.046292 

0.012984 

0.006132 

-0.000121 

*21 

0.136663 

0.052519 

0.037714 

0.022665 

A1 

%22 

0.742783 

0.743353 

0.749720 

0.753579 

*23 

-0.400459 

-0.386735 

-0.357284 

-0.347427 

*31 

0.110594 

0.120001 

0.127109 

0.128501 

*32 

-0.149914 

-0.194592 

-0.238972 

-0.255497 

*33 

2.868244 

3.218544 

3.534535 

3.648592 

£.  . 
44 

0.954261 

0.952611 

-0.953182 

0.953720 

45 

-0.010844 

-0.012781 

-0.015541 

-0.014560 

l.c 

46 

0.014555 

-0.001810 

0.004197 

0.015384 

E ,E 

*54 

-0.061719 

-0.052706 

-0.055985 

-0.062818 

X y 

*55 

0.644506 

0.628135 

0.632574 

0.632490 

*56 

0.121184 

0.156477 

0.121628 

0.108883 

*64 

0.121167 

0.123881 

0.111529 

0.095997 

*65 

-0.192356 

-0.221290 

-0.197932 

-0.189855 

*66 

1.015856 

1.039117 

1.054165 

1.061696 
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TABLE  AI-10 

CALCULATED  SYMMETRIZED  L INVERSE  MATRIX  FOR  CD  X 

— — i-  3 

(all  units  are  amu  z) 


Symmetry 

Block 

L"1. 

£11 

cd3f 

1. 356573 

1 

£12 

-0.105088 

-0 

£13 

-0.003254 

-0 

*21 

0.367464 

0 

£22 

0.748494 

0 

A1 

£23 

0.942466 

-0 

£31 

0.077287 

0 

£32 

-0.707577 

-0 

£33 

2.892391 

3 

CD  Cl 

CD  Br 

CD  I 

3 

3 

3 

382575 

1.387374 

1.390926 

013636 

0.003198 

0.019957 

009283 

-0.011860 

-0.015092 

186280 

0.155505 

0.126037 

941629 

0.956402 

0.968564 

061694 

-0.101719 

-0.150036 

214015 

0.231785 

0.237144 

482500 

-0.536983 

-0.551849 

449344 

3.830689 

3.971611 

44 

1.278629 

1.277682 

1.279922 

1.281858 

lAC 

45 

0.006776 

0.001431 

-0.002689 

-0.001947 

£46 

-0.005475 

-0.016886 

-0.006151 

0.011479 

£54 

-0.246505 

-0.199466 

-0.204490 

-0.210488 

£55 

0.911924 

0.885487 

0.884696 . 

0.880739 

£56 

-0.077719 

0.115974 

0.085422 

0.085892 

£64 

0.274651 

0.291368 

0.267305 

0.240734 

£65 

-0.024742 

-0.175585 

-0.157445 

-0.161177 

£66 

1.339117 

1.405252 

1.433064 

1.449007 
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TABLE  AI-11 

CALCULATED  HARMONIC 

WAVENUMBERS 

(cm  ■*" ) FOR 

ch3x 

Symmetry 

CH  F 

CH  Cl 

CH  Br 

CH  I 

Block 

3 

3 

3 

3 

W1 

3024.92 

3075.53 

3084.39 

3081.52 

Ai 

W2 

1488.11 

1384.62 

1331.35 

1277.14 

W3 

1062.74 

738.88 

617.66 

539.08 

to 

4 

3131.53 

3164.36 

3184.37 

3187.13 

E 

to 

5 

1493.47 

1481.47 

1472.11 

1464.81 

to 

6 

1207.14 

1038.97 

974.97 

902.48 
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TABLE  AI-12 


CALCULATED  HARMONIC  WAVENUMBERS  (cm”1)  for  CD3X 


Symmetry 

CD  F 

CD  Cl 

CD  Br 

3 

CD_  I 

Block 

3 

3 

3 

3 

03 

1 

2178.26 

2201.91 

2205.04 

2200.39 

fli 

U2 

1136.25 

1042.13 

1003.43 

962.63 

W3 

1009.14 

706.53 

582.65 

506.31 

03 

4 

2323.17 

2349.54 

2362.86 

2363.55 

E 

03 

5 

1086.01 

1072.36 

1068.61 

1064.06 

924.09 

780.29 

725.19 

667. 97 
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TABLE  AI-13 

DESCRIPTION  OF  THE  INFRARED  BAND  ASSIGNMENTS 
OF  THE  METHYL  HALIDES 

"Description  of  Motion" 
Symmetric  CH  stretch 
Symmetric  XCH  bend  (umbrel la- motion) 
Symmetric  CX  stretch 
Asymmetric  CH  stretch 
Asymmetric  HCH  bend 
Asymmetric  XCH  bend 


W1 

w2 


APPENDIX  II 


COMPUTER  PROGRAMS 


Program  PVDTEN 

Program  PVDTEN  calculated  the  P , V,  and  P tensors  for 

— s — —x 

a molecular  system  containing  up  to  sixteen  atoms.  The  pro- 
gram is  designed  to  accept  the  symmetrized  L ^ matrix  and 
symmetrized  B matrix  used  in  a normal  coordinate  analysis 
to  calculate  P^  and  V.  At  the  user's  option,  rotational 
corrections  are  computed  and  P is  calculated.  In  addition, 

X 

the  program  provides  options  for  calculating  L 1 from  I,  and 


symmetrizing  B.  An  option  for  calculating  G is  available 

as  an  internal  check  on  L ( L L'  = G)  and  B (B  M_1B ' = G) . 

An  option  for  calculating  A and  G 1 is  also  available.  The 

a 

program  will  punch  the  V and  P tensors  on  cards  in  the 

X 

correct  format  for  use  in  program  TRANSC  if  desired. 

Program  PVDTEN  consists  of  a main  program  and  four  sub- 
routines (ROTCOR,  LINV,  SDYM,  AMAT) . The  main  program  calcu- 
lates Pg  and  V directly  from  L ^ and  P^  and  the  symmetrized  B. 
The  subroutine  ROTCOR  calculates  the  rotational  corrections 
and  transfers  them  to  the  main  program  to  be  used  with  V 
for  the  P calculation.  Subroutine  LINV  is  called  if  L ^ is 

X 

to  be  calculated  from  L,  F,  and  A;  BSYM  is  called  if  B is  to 

be  symmetrized.  Subroutine  AMAT  calculated  A and  G_1  {G_1= 

(LL 1 ) ^ and  A=M  ^B'  G and  contains  an  option  for  calculating  It 
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and  from  a given  P^  matrix.  A description  of  the  in-put 
data  is  given  below: 

ITEM  (each  item  begins  a new  card) : 

1.  NOPROB  — The  number  of  the  problem.  If  the  number  of  the 
problem  is  7777,  the  program  terminates;  otherwise, 
problems  may  be  stacked  for  consecutive  runs.  FORMAT  (14). 

2.  RECORD (I)  --  Alphanumeric  data  for  identification  of  the 
problem.  The  first  column  must  be  blank.  FORMAT  (IX, 

12A6) . 

3.  NS,  NA,  NI , NRC  --  NS  is  the  number  of  independent 
symmetry  coordinates  exclusive  of  any  redundant  coordi- 
nates; NA  is  the  number  of  atoms;  NI  is  the  number  of 
internal  coordinates;  and  NRC  is  the  number  of  redundant 
coordinates.  FORMAT  (413)  . 

The  inclusion  of  NRC  is  an  artificial  device  so  that  the 
punched  B from  program  WMAT  (19)  can  be  utilized  in  PVDTEN 
if  redundant  coordinates  are  used  in  WMAT. 

4.  IFLI , IFSB , IFCHL , IFCHLI , IFCHB , IFROC,  IFPV,  IFPD,  IFAM  -- 

This  is  the  option  control  card.  IFLI  is  the  option  flag 

-1  -1 
for  calculating  L from  L,  F,  and  A.  If  IFLI  = 1,  L 

is  calculated  from  L by  subroutine  LINV.  If  IFLI  = 0, 

-1 

L is  read  directly  into  the  program.  IFSB  is  the 
option  flag  for  symmetrizing  B.  If  IFSB  = 1,  the  B 
matrix  is  symmetrized  by  subroutine  BSYM;  if  IFSB  = 0, 
the  symmetrized  B is  read  directly  into  the  program. 

IFCHL  is  the  option  flag  for  checking  L by  calculating 
G.  If  IFCHL  = 1,  G is  computed  from  L;  if  IFCHL  = 0, 
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G is  not  computed  from  L.  IFCHLI  is  the  option  flag  for 
checking  L 1 as  calculated  by  LINV.  If  IFCHLI  = 1,  L L_1 
is  computed  as  a check  on  L ^ ; if  IFCHLI  = 0 , no  check 
is  made.  The  options  IFCHL  and  IFCHLI  are  in  effect 
only  if  LINV  has  been  called  (IFLI  = 1) . IFCHB  is  the 
option  flag  for  checking  the  symmetrized  B matrix.  If 
IFCHB  = 1,  G is  computed,  and  if  IFCHB  =0,  it  is  not. 
This  option  is  in  effect  only  if  BSYM  has  been  called 
(IFBS  = 1) . IFROC  is  the  option  flag  for  calculating 
the  rotational  corrections.  If  IFROC  = 1,  rotational 
corrections  are  calculated  by  subroutine  ROTCOR,  and  P 
may  then  be  computed.  IFPV  and  IFPD  are  the  option  flags 

CL 

for  punching  the  V and  P tensors  on  cards  for  use  in 
program  TRANSC.  If  these  two  flags  are  equal  to  "1", 
the  tensors  are  punched;  if  the  flags  are  equal  to  "0", 
the  tensors  are  not  punched.  IFAM  is  the  option  flag 
for  calculating  A and  G ^ by  subroutine  AMAT.  In  order 
to  utilize  this  subroutine,  however,  IFSB  must  be  equal 
to  "1".  This  is  due  to  the  fact  that  G used  in  AMAT  is 
calculated  by  BSYM.  Hence,  if  IFSB  = 1 and  IFAM  = 1, 
then  A and  G ^ are  computed.  If  IFSB  = 0 and/or  IFAM 
= 0,  A and  G ^ are  not  computed.  FORMAT  (913). 

5.  XM(I)  — The  atomic  masses  are  entered  in  the  order  used 

in  the  normal  coordinate  analysis.  FORMAT  (6F12.6). 

^ -1  . 

If  L is  to  be  calculated  from  L,  IFLI  = 1,  the  following 

sets  of  data  are  required: 

6.  F ( I ) — The  frequencies  (cm  ) as  computed  by  the  normal 
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coordinate  analysis.  FORMAT  (6F12.6). 

7.  XL ( I , J ) — The  symmetrized  L matrix  is  entered  by 
reading  the  row  number,  column  number,  and  the  L 
element.  The  row  number  following  the  last  L element 
is  set  to  - 06  (end  code) . Only  non-zero  elements  need 
to  be  entered.  FORMAT  {4  (213 ,F12 . 6) } . 

8.  FC (I , J)  — The  upper  triangle  of  the  symmetrized  F matrix 
is  entered  using  the  same  format  as  item  7.  Only  non- 
zero elements  need  to  be  entered.  The  end  code  is  - 07. 

If  L ^ is  read  directly,  IFLI  = 0,  then  items  6,  7,  and  8 
are  omitted  and  only  the  symmetrized  L is  read  using  the 
same  format  as  item  7,  except  the  end  code  is  "0". 

If  B is  to  be  symmetrized,  IFSB  = 1,  the  following  sets 
of  data  are  required: 

9.  B(I,J)  --  The  non-symmetrized  B matrix  (as  punched  by 
WMAT)  is  read  using  the  same  format  as  item  7,  except 
the  end  code  is  - 05.  Only  non-zero  elements  need  to 
be  entered. 

10.  U(I,J)  --  The  non-normalized  U matrix  used  in  the  normal 
coordinate  analysis  is  entered  in  the  same  format  as 
item  7,  except  the  end  code  is  - 03.  Only  non-zero 
elements  need  to  be  entered. 

If  B is  already  symmetrized,  IFSB  = 0,  items  9 and  10  are 
omitted,  and  B is  read  directly  using  the  format  of  item  7, 
except  the  end  code  is  "0". 

11.  A(I,J)  — The  3p/0£b  matrix  is  entered  by  giving  the  row 
and  column  number  followed  by  the  absolute  value  of 
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3p/9Q.«  The  end  code  is  "0",  and  only  non-zero  elements 
need  to  be  entered.  FORMAT  {4  (213 ,F12 . 6) } . 

If  rotational  corrections  are  required,  the  following  sets  of 
data  are  read: 


12.  CC (I , J)  --  The  center  of  mass  Cartesian  coordinates  are 

read  giving  the  coordinate  number  (x  = 1,  y = 2,  z = 3), 
the  atom  number,  and  the  value  of  the  coordinate,  respec- 
tively. The  end  code  is  "0".  Only  non-zero  elements 
need  to  be  entered.  FORMAT  {4 (2I3,F12.6) }. 

0 0 o 

13.  u , u , u --  The  x,  y,  and  z components  of  the 

X y Z 

permanent  dipole  moment  are  entered  in  the  order  x,  y,  z. 


14. 


FORMAT  (3F12.6). 

XI,  YI,  ZI  — The  principal  moments  of  inertia  are  entered 
in  the  order  I , I , and  I with  the  same  format  as 

XX  yy  Z Z 


13. 


If  no  rotational  corrections  are  required,  IFROC  = 0,  items 
12,  13,  and  14  are  omitted. 

15.  SIGN (I, J)  — The  signs  (a.)  of  the  9p/3Q  's  are  entered 

1 i 

with  the  same  format  as  item  11. 

16.  IFADSC  — A flag  to  signify  whether  additional  sign  com- 
binations for  the  9p/3Q.'s  are  to  follow.  If  IFADSC  = 1, 

l 

an  additional  set  of  signs  are  read  using  the  same  for- 
mat as  item  15.  Thus,  any  number  of  sets  of  signs  may 
be  entered  after  the  initial  set  (item  15)  if  each 
different  set  is  preceeded  by  a card  with  IFADSC  = 1. 

If  IFADSC  = 0,  the  program  continues  to  the  calculation 
of  A and  G if  IFAM  = 1 or  proceeds  to  the  next  problem 
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if  IFAM  = 0. 

17.  IFCAL  --  An  option  flag  in  subroutine  AMAT  to  indicate 

that  P and  P are  to  be  calculated  from  a P matrix. 

~S  “Q  -X 

If  IFCAL  = 1,  a P matrix  is  to  follow  (item  18) . If 

X 

IFCAL  = 0,  the  program  proceeds  to  the  next  problem. 
FORMAT  (13). 

18.  CAL(I,J)  — -The  calculated  P matrix  for  the  subsequent 

X 

calculation  of  P and  P is  read  in  using  the  format 

S -Q 

of  item  11. 

Items  17  and  18  can  be  used  only  if  AMAT  has  been  called. 

To  terminate  the  program,  a problem  number  card  (item  1) 
is  placed  at  the  end  of  the  data  with  NOPROB  = 7777. 

PVDTEN  was  initially  written  to  use  symmetry  coordinates, 
but  internal  coordinates  could  be  used  if  certain  variables 
are  assigned  to  compensate  for  the  differences  between  the 
coordinates.  In  the  case  of  internal  coordinates,  NS  would 
be  equal  to  NI  and  a unit  matrix  could  be  entered  for  the 
required  U matrix  if  the  options  in  subroutines  BSYM  and 
AMAT  were  required. 

The  program  length  is  about  181  K bytes  and  could  easily 
be  expanded  to  accept  molecules  larger  than  sixteen  atoms  by 
increasing  the  dimensions  of  the  arrays. 

A computer  print-out  of  program  PVDTEN  is  given  below. 
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Program  TRANSC 

Program  TRANSC  rotates  a Cartesian  3x3  tensor  T about 
the  origin  of  the  Cartesian  coordinate  system.  The  trans- 
formation is  given  by 


3 


I 

j,l=l 


a . 

l 


. a t 
3 l,k  j,l 


(AII-1) 


where  t.  is  the  magnitude  of  the  i,  kth  element  of  the  new 

1 r K 

I _ I ->  I 

tensor  T . The  unit  vectors  e_^  and  e^  are  parallel  to  the 

new  i and  k coordinate  axes,  respectively.  The  coefficients 

a.  . and  a are  just  the  cosines  of  the  angles  between  the 
1 / J 1 / K 

original  coordinate  axes  and  the  new  coordinate  axes. 

Program  TRANSC  consists  of  a main  program  which  trans- 

I 

forms  the  tensor  T to  the  new  tensor,  T , and  one  subroutine, 

1 

ERROR, which  propagates  the  dispersions  in  T through  to  T . 

A description  of  the  in-put- data  is  given  below: 

ITEM  (each  item  begins  on  a new  card) : 

1.  NOPROB  --  The  number  of  the  problem.  If  NOPROB  = 7777, 
the  program  terminates;  otherwise,  problems  may  be 
stacked  for  consecutive  runs.  FORMAT  (14). 

2.  RECORD  (I)  — Alphanumeric  data  for  identification  of  the 
problem.  The  first  column  must  be  blank.  FORMAT  (IX, 


12A6) . 

3.  X(I,J)  --  The  nine  angles  (in  degrees)  between  the  old 

i i i 

and  new  axes  are  entered  in  the  order  x x,  x y,  x z, 

iii  i 

y x,  yy,  y z...z  z,  with  three  angles  per  card.  FORMAT 
(3F12.8) . 

4.  T(I,J)  — The  original  tensor  is  entered  in  row  order  with 


247 


three  elements  per  card.  FORMAT  (3F12.8). 

5.  IFERR  — If  dispersions  in  the  old  tensor  are  known  and 
the  resulting  dispersions  in  the  new  tensor  are  required, 
then  set  IFERR  = 1.  This  activates  subroutine  ERROR  for 
the  calculations  of  the  new  dispersions.  If  IFERR  = 0, 
no  dispersions  are  computed..  The  IFERR  card  is  included 
for  each  transformation.  FORMAT  (14) . 

6.  ERR ( I , J)  — If  IFERR  = 1,  the  dispersions  in  the  old 
tensor  are  read  using  the  same  format  as  item  4.  If 
IFERR  = 0,  the  dispersions  are  omitted. 

7.  IFADT  — If  additional  tensors  are  to  be  transformed 
using  the  original  angles  of  item  3,  then  the  IFADT  flag 
is  set  equal  to  "1".  Any  number  of  tensors  may  be  trans- 
formed using  the  original  angles  if  each  tensor  is  pro- 
ceeded by  a card  with  IFADT  = 1.  (Note,  however,  that 
the  IFERR  card,  item  5,  must  also  be  present  for  each 
additional  tensor.)  If  IFADT  = 0,  the  program  proceeds 
to  a new  problem.  FORMAT  (14) . 

TRANSC  is  about  4.5  K bytes  in  length.  A listing  of  the 


program  is  given  below. 


//TRANSC  JOB  (2003,  1078,020,05,0200),'  JAMES  NEWTON  ' ,CLASS  = 

//  EXEC  F4GCXS  FORT  G COMPILE  (NODECK),  EXECUTE,  CLASS  S 

//FORT.SYSIN  DC  * 
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55  CONTINUE 
WR1TE(6,58) 

58  FORMAT! IX, //IX, 'PRODUCT  OF  THE  NEW  TENSOR  WITH  ITS  TRANSPOSE') 
CO  59  1=1,3 

WRITE(6, 18)  ( I , J , TTT ( I , J ) , J = 1 , 3 ) 


CALCULATE  THE  SC.  ROOT  OF  THE  SUM  OF  THE  SQ.'  OF  THE  DIAGONAL 
ELEMENTS  OF  THE  PROCUCT  OF  THE  NEW  TENSOR  AND  ITS  TRANSPOSE 
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END 


APPENDIX  III 


CORIOLIS  CONSTANTS  FOR  THE  CH  X SERIES 

z y 

The  Coriolis  constants  z,  and  z,  for  the  methyl  halide 

.1  i,j 

series  were  calculated  by  the  procedure  outlined  by  Meal  and 
Polo  {J.H.  Meal  and  S.R.  Polo,  J.  Chem.  Phys.,  24_,  1119  (1956)}. 

T 

As  described  there,  the  £ (x  = x,  y,  or  z)  matrix  is  given  by 

x -1  x -1  i 

£ = L C (L  ) , (AIII-1) 

-1 

where  L is  the  symmetrized  inverse  normal  coordinate  trans- 
formation matrix  given  in  Appendix  I,  and  C is  defined  as 


T i 

C = D M D 


(AIII-2 ) 


3N  x 3N  matrix  composed  of  N 3 x 3 sub- 
diagonal as  shown  below. 
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The  submatrices  MX 


have  the  form 


0 

0 

0 

0 

0 

-i  \ 

l3x 

II 

0 

0 

1 

MY  = 

0 

0 

0 

. 0 

-1 

0 1 

1 

0 

0 / 

1 

0 

0 


0 1 
0 
0 


The  matrix  D is  the  mass  weighted  symmetrized  B (D 


B M 2) 


matrix  used  in  the  normal  coordinate  analysis.  (Do  not 

— ^ T 

confuse  M 2 with  M . ) 

Our  calculated  £ constants  are  given  in  Table  AIII-1. 

The  experimental  values  are  given  in  parentheses  where 

y 

available.  The  £ values  are  the  coupling  constants  between 

i 

z 

the  degenerate  E modes  and  the  £ values  are  the  coupling 

i/ j 

constants  between  the  ith  A mode  and  the  jth  E mode. 

— l,z  — 

One  should  be  aware  of  the  fact  that  the  zeta  constants 
are  not  free  of  sign  ambiguities,  for  as  shown  in  Equation 
AIII-1  and  AIII-2,  £ is  a function  of  both  L ^ and  B.  Thus, 
to  pick  a sign  of  the  zeta  constant,  one  must  know  the  L ^ 
matrix  (that  is,  the  sense  of  Q)  and  the  original  internal 
coordinate  definitions  (B) , since  the  sign  of  zeta  is  not 
derived  directly  from  experiment  (27,  34). 

We  are  primarily  interested  in  this  work  in  the  trends 

z z z 

in  £ and  £ . We  see  from  the  table  that  £ is  ex- 

2,5  2,6  2,5 

pected  to  be  almost  constant  through  the  CH^X  series,  and 
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values  of  £ are  expected  to  change  smoothly.  The  fact 

2,6 

that  there  does  not  appear  to  be  any  anomalous  behavior  in 


these  calculated  values  of  £ through  the  series  CH^X  prompts 

us  to  predict  that  the  sign  relationship  between  3p/3Q  and 

->  -*  -> 
dp/dQ^  , as  well  as  that  between  Sp/SQ^  and  3p/3Q^  , is  the 

same  for  the  entire  methyl  halide  series  as  found  for  CH^F 

and  CD^Cl  by  DiLauro  and  Mills  (15)  . 
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